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fT^ Abstract 

^ ■ This review is focused on tests of Einstein's theory of General Relativity with gravitational 

2i waves that are detectable by ground-based interferometers and pulsar timing experiments. 

. Einstein's theory has been greatly constrained in the quasi-linear, quasi-stationary regime, 

\l where gravity is weak and velocities are small. Gravitational waves will allow us to probe 

'>— ' a complimentary, yet previously unexplored regime: the non-linear and dynamical strong- 

_ ' field regime. Such a regime is, for example, applicable to compact binaries coalescing, where 

. . characteristic velocities can reach fifty percent the speed of light and compactnesses can reach 

^ a half. This review begins with the theoretical basis and the predicted gravitational wave 

observables of modified gravity theories. The review continues with a brief description of the 

^ . detectors, including both gravitational wave interferometers and pulsar timing arrays, leading 

J3 to a discussion of the data analysis formalism that is applicable for such tests. The review 

ends with a discussion of gravitational wave tests for compact binary systems. 

Contents 

1 Introduction [3] 

1.1 The Importance of Testing [3] 

1.2 Testing General Relativity versus Testing Alternative Theories |3] 

1.3 Gravitational Wave Tests versus Other Tests of General Relativity [S] 

1.4 Ground-Based vs Space-Based and Interferometers vs Pulsar Timing [5] 

1.5 Notation and Conventions [7] 



2 Alternative Theories of Gravity [8] 

2.1 Desirable Theoretical Properties [5] 

2.2 Well-Posedness and Effective Theories [S] 

2.3 Explored Theories [11] 

2.3.1 Scalar-Tensor Theories E] 

2.3.2 Massive Graviton Theories and Lorentz Violation [T3] 

2.3.3 Modified Quadratic Gravity [TH] 

2.3.4 Variable G Theories and Large Extra Dimensions [201 

2.3.5 Non- Commutative Geometry [531 

2.3.6 Gravitational Parity Violation [25| 

2.4 Currently Unexplored Theories in the Gravitational Wave Sector [27| 



3 Detectors 

3.1 Gravitational Wave Interferometers [28| 

3.2 Pulsar Timing Arrays [21] 



4 Testing Techniques 

4.1 Coalescence Analysis [36] 

4.1.1 Matched Filtering and Fisher Analysis [M] 

4.1.2 Bayesian Theory and Model Testing [37] 

4.1.3 Systematics in Model Selection [3S] 

4.2 Burst Analyses [10] 

4.3 Stochastic Background Searches [43] 

5 Compact Binary Tests [47] 

5.1 Direct and Generic Tests [17] 

5.2 Direct Tests [48] 

5.2.1 Scalar- Tensor Theories [15] 

5.2.2 Modified Quadratic Gravity [H 

5.2.3 Non-Commutative Geometry [57] 

5.3 Generic Tests [10] 

5.3.1 Massive Graviton Theories and Lorentz Violation [60] 

5.3.2 Variable G Theories and Large Extra-Dimensions [64] 

5.3.3 Parity Violation [HT] 

5.3.4 Parameterized Post-Einstcinian Framework [70] 

5.3.5 Searching for Non-Tcnsorial Gravitational Wave Polarizations [7S] 

5.3.6 LLove-Q Tests [78] 

5.4 Tests of the No-Hair Theorems [80] 

5.4.1 The No-Hah Theorems [50] 

5.4.2 Extreme Mass-Ratio Tests of The No-Hah Theorem [H] 

5.4.3 Ringdown Tests of The No-Hah Theorem [55] 

5.4.4 The Hairy Search for Exotica [57] 



6 Musings About the Future 

7 Acknowledgements 1911 



1 Introduction 

1.1 The Importance of Testing 

The era of precision gravitational wave astrophysics is at our doorstep. With it, a plethora of 
previously unavailable information will flood in, allowing for unprecedented astrophysical mea- 
surements and tests of fundamental theories. Nobody would question the importance of more 
precise astrophysical measurements, but one may wonder whether fundamental tests are truly nec- 
essary, considering the many successes of Einstein's theory of General Relativity. Indeed, General 
Relativity has passed many tests, including Solar System ones, binary pulsar ones and cosmological 
ones (for a recent review, see [399] ). 

What all of these tests have in common is that they sample the quasi-stationary, quasi-linear 
weak field regime of General Relativity. That is, they sample the regime of spacetime where the 
gravitational field is weak relative to the mass-energy of the system, the characteristic velocities 
are small relative to the speed of light, and the gravitational field is stationary or quasi-stationary. 
A direct consequence of this is that gravitational waves emitted by weakly-gravitating, quasi- 
stationary sources are necessarily extremely weak. To make this more concrete, let us define the 
gravitational compactness as a measure of the strength of the gravitational field: 

where M is the characteristic mass of the system, Jl is the characteristic length scale associated 
with gravitational radiation, and henceforth we set G = c = 1. For binary systems, the orbital 
separation serves as this characteristic length scale. The strength of gravitational waves and the 
mutual gravitational interaction between bodies scale linearly with this compactness measure. Let 
us also define the characteristic velocities of such a system V as a measure of the rate of change of 
the gravitational field. We can then more formally define the weak field as the region of spacetime 
where the following two conditions are simultaneously satisfied: 

WeakField: e<l, V<1. (2) 

By similarity, the strong field is defined as the region of spacetime where both conditions in Eq. (pi) 
are not valid simultaneousljF] 

Let us provide some examples. For the Earth-Sun system, M is essentially the mass of the sun, 
while 51 is the orbital separation, which leads to C w 9.8 x 10~^ and V sa 9.9 x 10^^. Even if an object 
were in a circular orbit at the surface of the Sun, its gravitational compactness would be O(10~^) 
and its characteristic velocity O(10~'^). We thus conclude that all Solar System experiments are 
necessarily sampling the weak field regime of gravity. Similarly, for the double binary pulsar J0737- 
3039 [213 1212], e w 6 X lO-*^ and V « 2 x IQ-^, where we have set the characteristic length 3? 
to the orbital separation via 3? ~ [MP^/(47r^)]^/'^ « 10^ km, where P = 0.1 days is the orbital 
period and M sa iM^ is the total mass. Although neutron stars are sources of strong gravity 
(the ratio of their mass to their radius is of order one tenth), binary pulsars are most sensitive 
to the quasi-static part of the post-Newtonian effective potential or to the leading-order (so-called 
Newtonian piece) of the radiation-reaction force. On the other hand, in compact binary coalescence 
the gravitational compactness and the characteristic speed can reach values much closer to unity; 
these sources would serve as strong-field probes of gravity. 

Even though current data does not give us access to the full non-linear and dynamical regime 
of General Relativity, Solar System tests and binary pulsar observations have served (and will 

^Notice that "strong-field" is not synonymous with Planck scale physics in this context. In fact, a stationary 
black hole would not serve as a probe of the strong-field, even if one were to somehow acquire information about the 
gravitational potential close to the singularity. This is because any such observation would necessarily be lacking 
information about the dynamical sector of the gravitational interaction. 



continue to serve) an invaluable role in testing Einstein's theory. Solar system tests effectively 
cured an outbreak of modified gravity theories in the 1970s and 1980s, as summarized for example 
in |399) . Binary pulsars were crucial as the first indirect detectors of gravitational waves, and later 
to heavily constrain theories that predict dipolar energy loss, as we will see in Sees. [2] and [Sj Binary 
pulsars are probes of General Relativity in a certain sector of the strong-field: in the dynamical 
but quasi-linear sector, verifying that compact objects move as described by a perturbative, post- 
Newtonian analysis to leading-order. Similarly, future electromagnetic observations of black hole 
accretion disks may probe General Relativity in another strong-field sector: the non-linear but 
fully stationary regime, verifying that black holes are described by the Kerr metric. As of the 
writing of this review article, only gravitational waves will allow for tests of General Relativity in 
the full strong-field regime, where gravity is both heavily non-linear and inherently dynamical. 

No experiments exist to date that validate Einstein's theory of General Relativity in the strong 
field region. Due to previous successes of General Relativity, one might consider such validation 
unnecessary. However, as most scientists would agree, the role of science is to predict and verify and 
not to assume without proof. Moreover, the incompleteness of General Relativity in the quantum 
regime, together with the somewhat unsatisfactory requirement of the dark sector of cosmology 
(including dark energy and dark matter), have prompted more than one physicist to consider 
deviations from General Relativity more seriously. Gravitational waves will soon allow us to verify 
Einstein's theory in a regime previously inaccessible to us, and as such, these tests are invaluable. 

In many areas of physics, however. General Relativity is so ingrained that questioning its va- 
lidity (even in a regime where Einstein's theory has not yet been validated) is synonymous with 
heresy. Dimensional arguments are usually employed to argue that any quantum gravitational 
correction will necessarily and unavoidably be unobservable with gravitational waves, as the for- 
mer are expected at a (Planck) scale that is inaccessible to gravitational wave detectors. This 
rationalization is dangerous, as it introduces a theoretical bias in the analysis of new observations 
of the Universe, thus quenching the potential for new discoveries. For example, if astrophysicists 
had followed such a rationalization when studying supernova data, they would not have discovered 
that the Universe is expanding. Dimensional arguments suggest that the cosmological constant 
is over 100 orders of magnitude larger than the value required to agree with observations. When 
observing the Universe for the first time in a completely new way, it seems more conservative to 
remain agnostic about what is expected and what is not, thus allowing the data itself to guide our 
efforts toward theoretically understanding the gravitational interaction. 

1.2 Testing General Relativity versus Testing Alternative Theories 

When testing General Relativity, one considers Einstein's theory as a null hypothesis and searches 
for generic deviations. On the other hand, when testing alternative theories one starts from a 
particular modified gravity model, develops its equations and solutions and then predicts certain 
observables that then might or might not agree with experiment. Similarly, one may define a 
bottom-up approach versus a top-down approach. In the former, one starts from some observables 
in an attempt to discover fundamental symmetries that may lead to a more complete theory, as 
was done when constructing the standard model of elementary particles. On the other hand, a 
top-down approach starts from some fundamental theory and then derives its consequence. 

Both approaches possess strengths and weaknesses. In the top-down approach one has complete 
control over the theory under study, being able to write down the full equations of motion, answer 
questions about well-posedness and stability of solutions, and predict observables. But as we shall 
see in Sec. [2] carrying out such an approach can be quixotic within any one model. What is worse, 
the lack of a complete and compelling alternative to General Relativity makes choosing a particular 
modified theory difficult. 

Given this, one might wish to attempt a bottom-up approach, where one considers a set of 



principles one wishes to test without explicit mention to any particular theory. One usually starts 
by assuming General Relativity as a null-hypothesis and then considers deformations away from 
General Relativity. The hope is that experiments will be sensitive to such deformations, thus 
either constraining the size of the deformations or pointing toward a possible inconsistency. But 
if experiments do confirm a General Relativity deviation, a bottom-up approach fails at providing 
a given particular action from which to derive such a deformation. In fact, there can be several 
actions that lead to similar deformations, all of which can be consistent with the data within its 
experimental uncertainties. 

Nonetheless, both approaches are complementary. The bottom-up approach draws inspiration 
from particular examples carried out in the top-down approach. Given a plausible measured 
deviation from General Relativity within a bottom-up approach, one will still need to understand 
what plausible top-down theories can lead to such deviations. From this standpoint, then, both 
approaches are intrinsically intertwined and worth pursuing. 

1.3 Gravitational Wave Tests versus Other Tests of General Relativity 

Gravitational wave tests differ from other tests of General Relativity in many ways. Perhaps one 
of the most important differences is the spacetime regime gravitational waves sample. Indeed, as 
already mentioned, gravitational waves have access to the most extreme gravitational environments 
in Nature. Moreover, gravitational waves travel essentially unimpeded from their source to Earth, 
and thus, they do not suffer from issues associated with obscuration. Gravitational waves also exist 
in the absence of luminous matter, thus allowing us to observe electromagnetically dark objects, 
such as black hole inspirals. 

This last point is particularly important as gravitational waves from inspiral black hole binaries 
are one of the cleanest astrophysical systems in Nature. In the last stages of inspiral, when such 
gravitational waves would be detectable by ground-based interferometers, the evolution of binary 
black holes is essentially unaffected by any other matter or electromagnetic fields present in the 
system. As such, one does not need to deal with uncertainties associated with astrophysical matter. 
Unlike other tests of General Relativity, such as those attempted with accretion disk observations, 
binary black hole gravitational wave tests may well be the cleanest probes of Einstein's theory. 

Of course, what is an advantage here, can be also a huge disadvantage in another context. 
Gravitational waves from compact binaries are intrinsically transient (they turn on for a certain 
amount of time and then shut off). This is unlike binary pulsar systems, for which astrophysicists 
have already collected tens of years of data. Moreover, gravitational wave tests rely on specific 
detections that cannot be anticipated before hand. This is in contrast to Earth-based laboratory 
experiments, where one has complete control over the experimental setup. Finally, the intrinsic 
feebleness of gravitational waves makes detection a very difficult task that requires complex data 
analysis algorithms to extract signals from the noise. As such, gravitational wave tests are limited 
by the signal-to-noise ratio and affected by systematics associated with the modeling of the waves, 
issues that are not as important in other loud astrophysical systems. 

1.4 Ground-Based vs Space-Based and Interferometers vs Pulsar Tim- 
ing 

This review article focuses only on ground-based detectors, by which we mean both gravitational 
wave interferometers, such as LIGO, Virgo and ET, as well as pulsar timing arrays (for a recent 
review of gravitational wave tests of General Relativity with space-based detectors, see |1641I407| V 
Ground-based detectors have the limitation of being contaminated by anthropomorphic noise, such 
as ground and air traffic, logging, earthquakes, ocean tides and waves, which arc clearly absent in 
space-based detectors. Ground-based detectors, however, have the clear benefit that they can be 



continuously upgraded and repaired in case of malfunction, which is obviously not possible with 
space-based detectors. 

As far as tests of General Relativity are concerned, there is a drastic difference in space-based 
and ground-based detectors: the gravitational wave frequencies these detectors are sensitive to. 
Space-based interferometers will be sensitive in the milli-Hz band, while ground-based interfer- 
ometers are sensitive in the deca- and hecta-Hz band. Different types of interferometers are then 
sensitive to different types of gravitational wave sources. For example, when considering binary 
coalescences, ground-based interferometers are sensitive to late inspirals and mergers of neutron 
stars and stellar-mass black holes, while space-based detectors will be sensitive to supermassive 
black hole binaries with masses around lO^Ai0. 

The impact of a different population of sources in tests of General Relativity depends on the 
particular modified gravity theory considered. When studying quadratic gravity theories, as we 
will see in Sec. [2] the Einstein-Hilbert action is modified by introducing higher order curvature 
operators, which are naturally suppressed by powers of the inverse of the radius of curvature. Thus, 
space-based detectors will not be ideal at constraining these theories, as the radius of curvature of 
supermassive black holes is much larger than that of stellar-mass black holes at merger. Moreover, 
space-based detectors will not be sensitive to neutron star coalescences, and thus, they cannot 
constrain theories that introduce modifications to General Relativity in non-vacuum spacetimes. 

On the other hand, space-based detectors are unique in their potential to probe the spacetime 
geometry of supermassive black holes through gravitational waves emitted during extreme mass- 
ratio inspirals. These inspirals consist of a stellar-mass compact object in a generic decaying orbit 
around a supermassive black hole. Such inspirals produce millions of cycles of gravitational waves in 
the sensitivity band of space-based detectors (in fact, they can easily out-live the observation time!). 
Therefore, even small changes to the radiation-reaction force, or to the background geometry, can 
lead to noticeable effects in the waveform observable and thus strong tests of General Relativity, 
albeit constrained to the radius of curvature of the supermassive black hole. For recent work on 
such systems and tests, see p2l [8M1 l335l [2M1 ITTHl l45l [2^^ 11681 ISMl l482l [28l [2M1 llB^ 

Space-based detectors also have the advantage of range, which is particularly important when 
considering theories where gravitons do not travel at light speed |284) . Space-based detectors have 
a horizon distance much larger than ground-based detectors; the former can see black hole mergers 
to redshifts of order 10 if there are any at such early times in the Universe, while the latter are 
confined to events within redshift 1. Gravitational waves emitted from distant regions in spacetime 
need a longer time to propagate from the source to the detectors. Thus, theories that modify the 
propagation of gravitational waves will be best constrained by space-based type systems. Of course, 
such theories are also likely to modify the generation of gravitational waves, which ground-based 
detectors should also be sensitive to. 

Another important difference between detectors is in their response to an impinging gravita- 
tional wave. Ground-based detectors, as we will see in Sec. [3j cannot separate between the two 
possible scalar modes (the longitudinal and the breathing modes) of metric theories of gravity, due 
to an intrinsic degeneracy in the response functions. Space-based detectors in principle also possess 
this degeneracy, but they may be able to break it through Doppler modulation if the interferometer 
orbits the Sun. Pulsar timing arrays, on the other hand, lack this degeneracy altogether, and thus, 
they can in principle constrain the existence of both modes independently. 

Pulsar timing arrays differ from interferometers in their potential to test General Relativity 
mostly by the frequency space they are most sensitive to. The latter can observe the late inspiral 
and merger of compact binaries, while the former is restricted to the very early inspiral. This 
is why pulsar timing arrays do not need very accurate waveform templates that account for the 
highly-dynamical and non-linear nature of gravity to detect gravitational waves; leading-order 
quadrupole waveforms are sufficient jll3| . In turn, this implies that pulsar timing arrays cannot 
constrain theories that only deviate significantly from General Relativity in the late inspiral, while 



they are exceptionally well-suited for constraining low-frequency deviations. 

We therefore see a complementarity emerging: different detectors can test General Relativity 
in different complementary regimes: 

• Ground-based detectors are best at constraining higher-curvature type modified theories that 
deviate from General Relativity the most in the late inspiral and merger phase. 

• Space-based detectors are best at constraining modified graviton dispersion relations and the 
geometry of supermassive compact objects. 



• 



Pulsar-timing arrays are best at independently constraining the existence of both scalar 
modes and any deviation from General Relativity that dominates at low orbital frequencies. 



Through the simultaneous implementation of all these tests, General Relativity can be put on a 
much firmer footing in all phases of the strong-field regime. 

1.5 Notation and Conventions 

We follow mainly the notation of [286) . where Greek indices stand for spacetime coordinates and 
spatial indices in the middle of the alphabet (i, j, fc, . . .) for spatial indices. Parenthesis and square 
brackets in index lists stand for symmetrization and anti-symmetrization respectively, e.g. ^(^^i/) = 
{Afj_^ + A^^)/2 and ^[,^1/] = {A^^i, — A^^)/2. Partial derivatives with respect to spacetime and 
spatial coordinates are denoted d^A ~ A^j^ and diA — Ai respectively. Covariant differentiation is 
denoted V^A = A.^^, multiple covariant derivatives V^'^" = V^V . . ., and the curved spacetime 
D'Alembertian DA — V^V^^. The determinant of the metric g^i, is g, R^ivSa is the Riemann 
tensor, R^i, is the Ricci tensor, and R is the Ricci scalar. The Levi-Civita tensor and symbol are 
g/ii/(5cr g^j^j ^^ll^5<J respectively, with e^^^s _ _|_2 y^ g^j^ orthonormal, positively oriented frame. We 
use geometric units {G — c — \) and the Einstein summation convention is implied. 

We will be mostly concerned with metric theories, where gravitational radiation is only defined 
much farther than a gravitational wave wavelength from the source. In this far- or radiation-zone, 
the metric tensor can be decomposed as 

with 77^1, the Minkowski metric and h^i, the metric perturbation. If the theory considered has 
additional fields 0, these can also be decomposed in the far-zone as 

^^^Q+i^, (4) 

with (f)Q the background value of the field and V' a perturbation. With such a decomposition, the 
field equations for the metric will usually be wave equations for the metric perturbation and for 
the field perturbation, in a suitable gauge. 



2 Alternative Theories of Gravity 

In this chapter, we discuss the many possible alternative theories that have been studied so far in 
the context of gravitational wave tests. We begin with a description of the theoretically desirable 
properties that such theories must have. We then proceed with a review of the theories so far 
explored as far as gravitational waves are concerned. We will leave out the description of many 
theories in this chapter, especially those which currently lack a gravitational wave analysis. We 
will conclude with a brief description of unexplored theories as possible avenues for future research. 

2.1 Desirable Theoretical Properties 

The space of possible theories is infinite, and thus, one is tempted to reduce it by considering 
a subspace that satisfies a certain number of properties. Although the number and details of 
such properties depend on the theorist's taste, there is at least one fundamental property that all 
scientists would agree on: 

1. Precision Tests. The theory must produce predictions that pass all Solar System, binary 
pulsar and experimental tests that have been carried out so far. 

This requirement can be further divided into the following: 

l.a General Relativity Limit. There must exist some limit, continuous or discontinuous, 
such as the weak- field one, in which the predictions of the theory are consistent with 
those of General Relativity within experimental precision. 

l.b Existence of Knovi^n Solutions |386j . The theory must admit solutions that cor- 
respond to observed phenomena, including but not limited to (nearly) flat spacetime, 
(nearly) Newtonian stars, and cosmological solutions. 

l.c Stability of Solutions |386j . The special solutions described in property (l.b) must 
be stable to small perturbations. For example, perturbations to (nearly) Newtonian 
stars, such as impact by asteroids, should not render such solutions unstable. 

Of course, these properties are not all necessarily independent, as the existence of a weak-field 
limit usually also implies the existence of known solutions. On the other hand, the mere existence 
of solutions does not necessarily imply that these are stable. 

In addition to these fundamental requirements, one might also wish to require that any new 
modified gravity theory possesses certain theoretical properties. These properties will vary depend- 
ing on the theorist, but the two most common ones are listed below: 

2. Well-Motivated from Fundamental Physics. There must be some fundamental theory 
or principle from which the modified theory (effective or not) derives. This fundamental 
theory would solve some fundamental problem in physics, such as late time acceleration or 
the incompatibility between quantum mechanics and General Relativity. 

3. Well-posed Initial Value Formulation [386| . A wide class of freely specifiable initial data 
must exist, such that there is a uniquely determined solution to the modified field equations 
that depends continuously on this data. 

The second property goes without saying at some level, as one expects modified gravity theory 
constructions to be motivated from some (perhaps yet incomplete) quantum gravitational descrip- 
tion of nature. As for the third property, the continuity requirement is necessary because otherwise 
the theory would lose predictive power, given that initial conditions can only be measured to a 
finite accuracy. Moreover, small changes in the initial data should not lead to solutions outside 



the causal future of the data; that is, causality must be preserved. Section |2.2| expands on this 
well-posedness property further. 

One might be concerned that property (2) automatically implies that any predicted deviation 
to astrophysical observables will be too small to be detectable. This argument usually goes as 
follows. Any quantum gravitational correction to the action will "naturally" introduce at least 
one new scale, and this, by dimensional analysis, must be the Planck scale. Since this scale is 
usually assumed to be larger than 1 TeV in natural units (or 10"'^^ meters in geometric units), 
gravitational wave observations will never be able to observe quantum gravitational modifications 
(see e.g. |140| for a similar argument). Although this might be true, in our view such arguments 
can be extremely dangerous, since they induce a certain theoretical bias in the search for new 
phenomena. For example, let us consider the supernova observations of the late time expansion of 
the universe that led to the discovery of the cosmological constant. The above argument certainly 
fails for the cosmological constant, which on dimensional arguments is over 100 orders of magnitude 
too small. If the supernova teams had respected this argument, they would not have searched for 
a cosmological constant in their data. Today, we try to explain our way out of the failure of such 
dimensional arguments by claiming that there must be some exquisite cancelation that renders the 
cosmological constant small; but this, of course, came only after the constant had been measured. 
One is not trying to argue here that cancelations of this type are common and that quantum 
gravitational modifications are necessarily expected in gravitational wave observations. Rather, 
we are arguing that one should remain agnostic about what is expected and what is not, and allow 
oneself to be surprised without quenching the potential for new discoveries that will accompany 
the new era of gravitational wave astrophysics. 

One last property that we will wish to consider for the purposes of this review is the following: 

4. Strong Field Inconsistency. The theory must lead to observable deviations from General 
Relativity in the strong-field regime. 

Many modified gravity models have been proposed that pose infrared or cosmological modifications 
to General Relativity, aimed at explaining certain astrophysical or cosmological observables, like 
the late expansion of the Universe. Such modified models usually reduce to General Relativity in 
the strong-field regime, for example via a Vainshtein like mechanism |373l I125J . Gravitational wave 
tests of General Relativity, however, are concerned with modified theories that predict deviations 
in the strong-field, precisely where cosmological modified models do not. Clearly, property (4) is 
not necessary for a theory to be a valid description of nature, but it is a necessary feature if one 
is to test this theory with gravitational wave observations. 

2.2 Well-Posedness and Effective Theories 

Property (3) does not only require the existence of an initial value formulation, but also that it 
be well-posed, which is not necessarily guaranteed. For example, the Cauchy-Kowalewski theorem 
states that a system of n partial differential equations for n unknown functions (f>i of the form 
0i,tt = Fiix^] (^j.fi'i 4'j,ti'T 4'j,ik)j with Fi analytic functions has an initial value formulation (see 
e.g. [387] V This theorem, however, does not guarantee continuity or the causal conditions described 
above. For this, one has to rely on more general energy arguments, for example constructing a 
suitable energy measure that obeys the dominant energy condition and using it to show well- 
posedness (see e.g. |197[ I387J ). One can show that second-order, hyperbohc partial differential 
equations, i.e. equations of the form 

V^V^0 + A^V^0 + B0 + C = O, (5) 

where A^ is an arbitrary vector field and {B, C) are smooth functions, have a well-posed initial 
value formulation. Moreover, the Leray theorem proves that any quasilinear, diagonal, second- 
order hyperbolic system also has a well-posed initial value formulation 



Proving the well-posedness of an initial value formulation for systems of higher-than-second- 
order, partial differential equations is much more difficult. In fact, to our knowledge, no general 
theorems exist of the type described above that apply to third, fourth or higher-order, partial, 
non-linear and coupled differential equations. Usually, one resorts to the Ostrogradski theorem 
to rule out (or at the very least cast serious doubt on) theories that lead to such higher-order 
field equations. Ostrogradski's theorem states that Lagrangians that contain terms with higher 
than first time-derivatives possess a linear instability in the Hamiltonian (see e.g. 1404 1 for a nice 
reviewr] As an example, consider the Lagrangian density 



whose equations of motion 



m T muj o qm n , , 

q + UJ^q^^-^"q , (7) 

obviously contains higher derivatives. The exact solution to this differential equation is 

q~Ai cos kit + Bi sin kit + A2 cos k2t -f B2 sin /c2i , (8) 

where {Ai, Bi) are constants and k\ 2/^^ = (1 T \/r^^4^)/(2g). The on-shell Hamiltonian is then 

H = "^VT^Igkl {Al + Bl) "^V^^Tgki {Aj + Sf) , (9) 

from which it is clear that the (1) mode carries positive energy, while the (2) mode carries negative 
energy and forces the Hamiltonian to be unbounded from below. The latter implies that dynamical 
degrees of freedom can reach arbitrarily negative energy states. If interactions are present, then 
an "empty" state would instantaneously decay into a collection of positive and negative energy 
particles, which cannot describe the Universe we live in [404] . 

The Ostrogradski theorem, however, can be evaded if the Lagrangian in Eq. ^ describes an 
effective theory, i.e. a theory that is a truncation of a more general or complete theory. Let us 
reconsider the particular example above, assuming now that the coupling constant g is an effective 
theory parameter and Eq. m is only valid to linear order in g. One approach is to search for 
perturbative solutions of the form gpert — xq + gxi + . . ., which leads to the system of differential 
equations 

Xn+Uj'^Xn^ ^'i'ri-l, (10) 

with x_i =0. Solving these set of n differential equations and resumming, one finds 

^port = ^1 cosfcii + i?i sinfcii . (11) 

Notice that gpdt contains only the positive (well-behaved) energy solution of Eq. (Is]), ie. pertur- 
bation theory acts to retain only the well-behaved, stable solution of the full theory in the 5^-0 
limit. One can also think of the perturbative theory as the full theory with additional constraints, 
ie. the removal of unstable modes, which is why such an analysis is sometimes called perturbative 
constraints [TTUl [TTTl li?7] . 

Another way to approach effective field theories that lead to equations of motion with higher- 
order derivatives is to apply the method of order-reduction. In this method, one substitutes the 
low-order derivatives of the field equations into the high-order derivative part, thus rendering the 



■^Stability and well-posedness are not the same concepts and they do not necessarily imply each other. For 
example, a well-posed theory might have stable and unstable solutions. For ill-posed theories, it does not make 
sense to talk about stability of solutions. 
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resulting new theory usually well-posed. One can think of this as a series resummation, where 
one changes the non-linear behavior of a function by adding uncontrolled, higher-order terms. Let 
us provide an explicit example by reconsidering the theory in Eq. ([6]). To lowest order in g, the 
equation of motion is that of a simple harmonic oscillator, 

q + iu^q = 0{g), (12) 

which is obviously well-posed. One can then order-reduce the full equation of motion, Eq. (l7| , by 
substituting Eq. (12) into the right hand side of Eq. (l7|. Doing so, one obtains the order-reduced 



equation of motion 

q + u'q = gq + Oig^), (13) 

which clearly now has no high-order derivatives and is well posed provided g -^ 1. The solution 
to this order-reduced differential equation is qpcrt once more, but with ki linearized in g <C 1. 
Therefore, the solutions obtained with a perturbative decomposition and with the order-reduced 
equation of motion are the same to linear order in g. Of course, since an effective field theory is 
only defined to a certain order in its perturbative parameter, both treatments are equally valid, 
with the unstable mode effectively removed in both cases. 

Such a perturbative analysis, however, can say nothing about the well-posedness of the full 
theory from which the effective theory derives, or of the effective theory if treated as an exact one 
(ie. not as a perturbative expansion). In fact, a well-posed full theory may have both stable and 
unstable solutions. The arguments presented above only discuss the stability of solutions in an 
effective theory, and thus, they are self-consistent only within their perturbative scheme. A full 
theory may have non-perturbative instabilities, but these can only be studied once one has a full 
(non-truncated in g) theory, from which Eq. ([6]) derives as a truncated expansion. Lacking a full 
quantum theory of Nature, quantum gravitational models are usually studied in a truncated low- 
energy expansion, where the leading-order piece is General Relativity and higher order pieces are 
multiplied by a small coupling constant. One can perturbatively explore the well-behaved sector 
of the truncated theory about solutions to the leading-order theory. Such an analysis, however, is 
incapable of answering questions about well-posedness or non-linear stability of the full theory. 

2.3 Explored Theories 

In this subsection we briefly describe the theories that have so far been studied in some depth 
as far as gravitational waves are concerned. In particular, we focus only on those theories that 
have been sufficiently studied so that predictions of the expected gravitational waveforms (the 
observables of gravitational wave detectors) have been obtained for at least a typical source, such 
as the quasi-circular inspiral of a compact binary. 

2.3.1 Scalar- Tensor Theories 

Scalar-tensor theories in the Einstein frame |150l 11491 I161J are defined by the action (where we 
will restore Newton's gravitational constant G in this section) 



c(E) _ 



^ J d^X^-g [R ~ 2 (5^<^) (d^cp) - Viifi)] + 5,„at[^mat, ^'('^).9m.], (14) 



where ip is a scalar field, A{(p) is a coupling function, V{ip) is a potential function and i/'mat 
represents matter degrees of freedom. For more details on this theory, we refer the interested reader 



to the reviews [3991 1396] . One can of course consider more complicated scalar-tensor theories, for 
example by multiplying the action by arbitrary functionals of ip |426j . but we will ignore such 
generalizations here. 



11 



The Einstein frame is not the frame where the metric governs clocks and rods, and thus, 
it is convenient to recast the theory in the Jordan frame through the conformal transformation 






i?_^(aM0)(a^^)_0V 



+ <S'mat[V'mat,5/ii/], (15) 



where g^i, is the physical metric, the new scalar field (/> is defined via cj) = A~ , the coupling field is 
uj{(j)) = (a~^ — 3)/2 and a = A^^/A. When cast in the Jordan frame, it is clear that scalar-tensor 
theories are metric theories (see |399| for a definition), since the matter sector depends only on 
matter degrees of freedom and the physical metric (without a direct coupling of the scalar field) . 
When the coupling w(0) = wbd is constant, then Eq. (Ils]) reduces to the massless version of 
Jordan-Fierz-Brans-Dicke theory |76j . 

The modified field equations in the Einstein frame are 

4 dip Sip 

G^.=87rGfc* + r^^)) , (16) 



where 

, - 1 r 1 .1 1 

(17) 



^^ 47r 



^,i,^^,v - i^gt,v^,&^'^ - -^9t^uV{ip) 



is a stress-energy tensor for the scalar field. The matter-stress energy tensor is not constructed 
from the Einstein-frame metric alone, but by the combination A{ipYg^^. In the Jordan frame and 
neglecting the potential, the modified field equations are [3961 [84] 



diiJ 



U<t) = — - SttT 



..M 



G^. = ^T;r + ^ (0.^0,. - \9,ur^,}j + ^ (0.M- - ~9^..n<P) , (18) 

where T™**' is the trace of the matter stress-energy tensor T^Jf* constructed from the physical 
metric g^^- The form of the modified field equations in Jordan frame suggest that in the weak-field 
limit one may consider scalar-tensor theories as modifying Newton's gravitational constant via 
G -^ G{(j)) = G/(j). 

Using the decompositions of Eqs. (plj-Q, the field equations of massless Jordan-Fierz-Brans- 
Dicke theory can be linearized in the Jordan frame to find (see e.g. |402| ) 

U^ef"" = -Wttt"-' , D^V = -167r5 , (19) 

where D^ is the D'Alembertian operator of flat spacetime, we have defined a new metric pertur- 
bation 

e^"" = h^"' - -ri^'^h - ^T^^"" , (20) 

2 00 

with h the trace of the metric perturbation and 

r^'' = 0olTCt+i^^ (21) 
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with cubic remainders in either the metric perturbation or the scalar perturbation. The quantity 
i^" is a function of quadratic or higher order in 9'^'^ or ip. These equations can now be solved given 
a particular physical system, as done for quasi-circular binaries in |402l 13381 1302) . Given the above 
evolution equations, Jordan-Fierz-Brans-Dicke theory possesses a scalar (spin-0) mode, in addition 
to the two transverse-traceless (spin-2) modes of General Relativity, ie. Jordan-Fierz-Brans-Dicke 
theory is of Type iVg in the E{2) classification P^[5M] . 

Let us now discuss whether scalar-tensor theories satisfy the properties discussed in Sec. |2.1| 
Massless Jordan-Fierz-Brans-Dicke theory agrees with all known experimental tests provided wbd > 
4 X 10''', a bound imposed by the tracking of the Cassini spacecraft [BS]. Massive Jordan-Fierz- 
Brans-Dicke theory has been recently constrained to wbd > 4 x 10^ and nis < 2.5 x 10~^° eV, with 
TOs the mass of the scalar field, through observations of the Shapiro time delay [18]. Observations 
of the Nordtvedt effect with Lunar Laser Ranging observations, as well as observations of the 
orbital period derivative of white-dwarf /neutron-star binaries, yield less stringent constraints |18) . 
Neglecting any homogeneous solutions to the scalar-field evolution equation, it is clear that in the 
limit w — )■ oo one recovers General Relativity, ie. scalar-tensor theories have a continuous limit to 
Einstein's theory. Moreover, Refs. |339|, 1248] have verified that scalar-tensor theories with minimal 
or non-minimal coupling in the Jordan frame can be cast in a strongly-hyperbolic form, and thus, 
they possess a well-posed initial value formulation. Therefore, scalar-tensor theories possess both 
properties (1) and (3). 

Scalar-tensor theories also possess property (2), since they can be derived from the low-energy 
limit of certain string theories. The integration of string quantum fluctuations leads to a higher- 
dimensional string theoretical action that reduces locally to a fleld theory similar to a scalar-tensor 
one |169l 1157] . the mapping being cj) — e^^'^', with tp one of the string moduli flelds. Moreover, 
scalar-tensor theories can be mapped to f{R) theories, where one replaces the Ricci scalar by some 
functional of R. In particular, one can show that f{R) theories are equivalent to Brans-Dicke 
theory with wbd = 0, via the mapping = df{R)/dR and ¥{(/)) = R df{R)/dR - f{R) [9711357]. 
For a recent Living Reviews article on this topic, see |151| . 

Black holes and stars continue to exist in scalar-tensor theories. Stellar configurations are 
modified from their General Relativity profile [4021 1206] . while black holes are not, provided one 
neglects homogeneous solutions to the scalar field evolution equation. Indeed, Hawking [196] has 
proved that Brans-Dicke black holes that are stationary and the cndpoint of gravitational collapse 
are identical to those of General Relativity. This proof has recently been extended to a general 
class of scalar-tensor models |359j . That is, stationary black holes radiate any excess "hair", 
ie. additional degrees of freedom, after gravitational collapse, a result sometimes referred to as the 
no-hair theorem for black holes in scalar-tensor theories. This result has recently been extended 
even further to allow for quasi-stationary scenarios in generic scalar-tensor theories through the 
study of extreme-mass ratio inspirals |426j (small black hole in orbit around a much larger one) , 
post-Newtonian comparable-mass inspirals [283j and numerical simulations of comparable-mass 
black hole mergers [201]. 

A small caveat must however, be discussed. If homogeneous solutions to the scalar field are 
allowed (e.g. (p = 4'^x'-''), then it is not completely clear that generic scalar-tensor theories retain 
property (1). That is, if the scalar field is allowed to vary at asymptotic spatial infinity, then, 
Jacobson's work [217j implies that spherically symmetric exterior, gravitational fields would be 
modified from the Schwarzschild solution, acquiring a certain hair. Indeed, Ref. |207j recently used 
this fact to argue that when such hair is present, black hole binaries would radiate in a dipolar 
fashion, just like neutron star binaries in scalar-tensor theories do, thus greatly altering their binary 
evolution. 

Another caveat, of course, is that one is working with Jordan-Fierz-Brans-Dicke theory, and 
not a more general class of scalar-tensor theories. Damour and Esposito-Farese |1221I123] proposed 



a different type of scalar-tensor theory, one that can be defined by the action in Eq. ( 15 1 but with 
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the conformal factor A{^) = e^'^ /^ or the couphng function uj{(t)) — —3/2 — 27rG/(/31og(/)), where 
f3 is another constant. When /3 < —4, non-perturbative effects that develop if the gravitational 
energy is large enough can force neutron stars to spontaneously acquire a non-trivial scalar field 
profile, to spontaneously scalarize. Through this process, binary neutron stars that initially had 
no scalar hair in their early inspiral would acquire it before they merge, when their binding energy 
exceeds some threshold |46) . 

As for property (4), scalar tensor theories are not built to introduce high-order curvature cor- 
rections to General Relativity, and thus, they are not truly strong-field tests of General Relativitjj^ 
Instead, they naturally lead to modifications of Einstein's theory in the weak-field, modifications 
that dominate in scenarios with sufficiently weak gravitational interactions. Although this might 
seem strange, it is natural if one considers, for example, one of the key modifications introduced by 
scalar-tensor theories: the emission of dipolar gravitational radiation. Such dipolar emission dom- 
inates over the General Relativistic quadrupolar emission for systems in the weak to intermediate 
field regime, such as in binary pulsars or in the very early inspiral of compact binaries. Therefore, 
one would expect scalar-tensor theories to be best constrained by experiments or observations of 
weakly-gravitating systems, as it has recently been explicitly shown in |426) . 

2.3.2 Massive Graviton Theories and Lorentz Violation 

Massive graviton theories are those in which the gravitational interaction is propagated by a massive 
gauge boson, ie. a graviton with mass rUg 7^ or Compton wavelength Xg = h/(rngc) < 00. 
Einstein's theory predicts massless gravitons and thus gravitational propagation at light speed, 
but if this were not the case, then a certain delay would develop between electromagnetic and 
gravitational signals emitted simultaneously at the source. Fierz and Pauli |154j were the first to 
write down an action for a free massive graviton, and ever since then, much work has gone into 
the construction of such models. For a detailed review, see e.g. |203j . 

Gravitational theories with massive gravitons are somewhat well-motivated from a fundamental 
physics perspective, and thus, one can say they possess property (2). Indeed, in loop quantum 
cosmology [351 [H] j the cosmological extension to loop quantum gravity, the graviton dispersion 
relation acquires holonomy corrections during loop quantization that endow the graviton with a 
mass [72] rUg = A^^/'^^^^{p/pc), with 7 the Barbero-Immirzi parameter, A the area operator, 
and p and pc the total and critical energy density respectively. In string theory-inspired effective 
theories, such as Dvali's compact, extra-dimensional theory |142j such massive modes also arise. 

Massive graviton modes also occur in many other modified gravity models. In Rosen's bimetric 
theory |329j . for example, photons and gravitons follow null geodesies of different metrics [39911396) . 
In Visser's massive graviton theory |385j . the graviton is given a mass at the level of the action 
through an effective perturbative description of gravity, at the cost of introducing a non-dynamical 
background metric, i.e. a prior geometry. A recent re-incarnation of this model goes by the name 
of Bigravity, where again two metric tensors are introduced [3131 I311J . In Bekenstein's Tensor- 
Vector-Scalar (TeVeS) theory [3^, the existence of a scalar and a vector field lead to subluminal 
GW propagation. 

Massive graviton theories have a theoretical issue, the van Dam-Veltman-Zakharov (vDVZ) 
discontinuity |3781I436] . which is associated with property l.a, ie. a General Relativity limit. The 
problem is that certain predictions of massive graviton theories do not reduce to those of General 
relativity in the nig — >■ limit. This can be understood qualitatively by studying how the 5 
spin states of the graviton behave in this limit. Two of them become the two General Relativity 
helicity states of the massless graviton. Another two become helicity states of a massless vector 
that decouples from the tensor perturbations in the rUg — >■ limit. The last state, the scalar mode. 



•^An exception to this is during spontaneous scalarization |122l 11231 through which non-pcrturbative corrections 
can affect the structure of neutron stars in certain classes of scalar-tensor theories. 
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however, retains a finite coupling to the trace of the stress-energy tensor in this hniit. Therefore, 
massive graviton theories in the uig — >■ hmit do not reduce to General Relativity, since the scalar 
mode does not decouple. 

The vDVZ discontinuity, however, can be evaded, for example by carefully including non- 
linearities. Vainshtein |373( I238[ I125J showed that around any source of mass M, there exists a 
certain radius r < ry ^ (rgA^)^'^, with rg the Schwarzschild radius, where linear theory cannot 
be trusted. Since ry — ;> oo as rrig — )■ 0, this implies that there is no radius at which the linear 
approximation (and thus vDVZ discontinuity) can be trusted. Of course, to determine then whether 
massive graviton theories have a continuos General Relativity limit, one must include non- linear 
corrections to the action (see also an argument by |31j). which are more difficult to uniquely 
predict from fundamental theory. Recently, there has been much activity in the development of 
new, non-linear massive gravity theories [SHI 11241 11891 [55] . 

Lacking a particular action for massive graviton theories that modifies the strong-field regime 
and is free of non-linear and radiatively-induced ghosts, it is difficult to ascertain many of its 
properties, but this does not prevent us from considering certain phenomenological effects. If the 
graviton is truly massive, whatever the action may be, two main modifications to Einstein's theory 
will be introduced: 

(i) Modification to Newton's Third law; 

(ii) Modification to gravitational wave propagation. 

Modifications of class (i) correspond to the replacement of the Newtonian potential by a Yukawa 
type potential (in the non-radiative, near-zone of any body of mass M): V — (M/r) — >■ [M /r) exp(— r/Ag), 
where r is the distance to the body |398| . Tests of such a Yukawa interaction have been proposed 
through observations of bound clusters, tidal interactions between galaxies jl77| and weak gravi- 
tational lensing [99], but such tests are model-dependent. 

Modifications of class (ii) are in the form of a non-zero graviton mass that induces a modified 
gravitational wave dispersion relation. Such a modification to the dispersion relation was originally 
parameterized via |398| 

2 2 4 

where Vg and mg are the speed and mass of the graviton, while E is its energy, usually associated 
to its frequency via the quantum mechanical relation E ^ hf. This modified dispersion relation is 
inspired by special relativity, a more general version of which, inspired by quantum gravitational 
theories, is |284j 

^ = 1 - A" , (24) 

where a is now a parameter that depends on the theory and A represents deviations from light 
speed propagation. For example, in Rosen's bimetric theory [329 , the graviton does not travel at 
the speed of light, but at some other speed partially determined by the prior geometry. In metric 
theories of gravity, A = Arnlc^ / E'^ , where A is some amplitude that depends on the metric theory 
(see discussion in |284p . Either modification to the dispersion relation has the net effect of slowing 
gravitons down, such that there is a difference in the time of arrival of photons and gravitons. 
Moreover, such an energy-dependent dispersion relation would also affect the accumulated gravi- 
tational wave phase observed at gravitational wave detectors, as we will discuss in Sec. [5] Given 
these modifications to the dispersion relation, one would expect the generation of gravitational 
waves to also be greatly affected in such theories, but again, lacking a particular healthy action to 
consider, this topic remains today mostly unexplored. 

From the structure of the above phenomenological modifications, it is clear that General Rela- 
tivity can be recovered in the nig — >■ limit, avoiding the vDVZ issue altogether by construction. 
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Such phenomenological modifications have been constrained by several types of experiments and 
observations. Using the modification to Newton's third law and precise observations of the motion 
of the inner planets of the Solar System together with Kepler's third law, |398| found a bound of 
Xg > 2.8 X 10^^ km. Such a constraint is purely static, as it does not sample the radiative sector 
of the theory. Dynamical constraints, however, do exist: through observations of the decay of the 
orbital period of binary pulsars, |156| found a bound of Xg > 1.6 x 10"'^° km. New constraints that 
use gravitational waves have been proposed, including measuring a difference in time of arrival of 
electromagnetic and gravitational waves |1191 1236) , as well as direct observation of gravitational 
waves emitted by binary pulsars (see Sec. pi) . 

Although massive gravity theories unavoidably lead to a modification to the graviton dispersion 
relation, the converse is not necessarily true. A modification of the dispersion relation is usually 
accompanied by a modification to either the Lorentz group or its action in real or momentum space. 
Such Lorentz-violating effects are commonly found in quantum gravitational theories, including 
loop quantum gravity 170] and string theory |1001 1364) , as well as other effective models [531 \EM ■ 
In Doubly Special Relativity [23l 12691 [24l [25] , the graviton dispersion relation is modified at high 
energies by modifying the law of transformation of inertial observers. Modified graviton dispersion 
relations have also been shown to arise in generic extra-dimensional models |344) , in Hof ava-Lifshitz 
theory J 2041 12051 13721 IB5] and in theories with non-commutative geometries |1661ll67lll68) . None of 
these theories necessarily requires a massive graviton, but rather the modification to the dispersion 
relation is introduced due to Lorentz-violating effects. 

One might be concerned that the mass of the graviton and subsequent modifications to the gravi- 
ton dispersion relation should be suppressed by the Planck scale. However, Collins, et al. |104[|103] 
have suggested that Lorentz violations in perturbative quantum field theories could be dramatically 
enhanced when one regularizes and renormalizes them. This is because terms that vanish upon 
renormalization due to Lorentz invariance do not vanish in Lorentz-violating theories, thus lead- 
ing to an enhancement |165j . Whether such an enhancement is truly present cannot be currently 
ascertained. 

2.3.3 Modified Quadratic Gravity 

Modified quadratic gravity is a family of models first discussed in the context of black holes and 
gravitational waves in |434L 1408] . The 4-dimensional action is given by 

+ a4f4WR^u5a*R^''''' - ^ [V^l?V^l9 + 2Vm + £,nat| • (25) 

The quantity *R^uS<j = {^l'^)(-s<j°''^ R'^yap is the dual to the Riemann tensor. The quantity ^mat 
is the external matter Lagrangian, while /i(-) are functionals of the field d, with (a^,/?) coupling 



constants and k = (167rG')~^. Clearly, the two terms second to last in Eq. ( 25 ) represent a canonical 
kinetic energy term and a potential. At this stage, one might be tempted to set /3 = 1 or the ai — 1 
via a rescaling of the scalar field functional, but we shall not do so here. 



The action in Eq. (25) is well-motivated from fundamental theories, as it contains all possible 
quadratic, algebraic curvature scalars with running (i.e. non-constant) couplings. The only restric- 
tion here is that all quadratic terms are assumed to couple to the same field, which need not be 
the case. For example, in string theory some terms might couple to the dilaton (a scalar field), 
while other couple to the axion (a pseudo scalar field). Nevertheless, one can recover well-known 
and motivated modified gravity theories in simple cases. For example. Dynamical Chern-Simons 
modified gravity [14| is recovered when 04 = — acs/4 and all other a.i = 0. Einstein-Dilaton- 
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Gauss-Bonnet gravity |308j is obtained when 014 = and (ai,Q;2,a3) = (1, — 4, l)aEDGEr| Both 
theories unavoidably arise as low-energy expansions of heterotic string theory |180l I181l1l31 [5^ . 
As such, modified quadratic gravity theories should be treated as a class of effective field theories. 
Moreover, Dynamical Chern-Simons Gravity also arises in loop quantum gravity |391 I330J when the 
Barbero-Immirzi parameter is promoted to a field in the presence of fermions [37l [131 1367| I279[ I172J 
One should make a clean and clear distinction between the theory defined by the action of 



Eq. (251 and that of f{R) theories. The latter are defined as functionals of the Ricci scalar only, 



while Eq. ( 25 1 contains terms proportional to the Ricci tensor and Riemann tensor squared. One 
could think of the subclass of f{R) theories with f{R) = R^ as the limit of modified quadratic 
gravity with only ai 7^ and /i(^) — 1. In that very special case, one can map quadratic gravity 
theories and f{R) gravity to a scalar-tensor theory. Another important distinction is that f{R) 
theories are usually treated as exact, while the action presented above is to be interpreted as an 
effective theory, ie. a truncation to quadratic order in the curvature of a more fundamental theory. 
This implies that there are cubic, quartic, etc. terms in the Riemann tensor that are not included 



in Eq. (25 1 and that presumably depend on higher powers of a^. Thus, when studying such an 
effective theory one should also order-reduce the field equations and treat all quantities that depend 
on ai perturbatively, the so-called small- coupling approximation. One can show that such an order 
reduction removes any additional polarization modes in propagating metric perturbations |3531I361"| 
that naturally arise in f{R) theories. In analogy to the treatment of the Ostrogradski instability 



in Sec. 2.1 one would also expect that order-reduction would lead to a theory with a well-posed 
initial value formulation. 

This family of theories is usually simplified by making the assumption that coupling functions 
fi{-) admit a Taylor expansion: /,(z9) = /,(0) + /;(0)i? + Oid'^) for smaU d, where /^(O) and //(O) 
are constants and •& is assumed to vanish at asymptotic spatial infinity. Reabsorbing /i(0) into the 
couphng constants a- = ai/i(0) and /-(O) into the constants al = aif^O), Eq. (25) becomes 



S = Sqr -I- 5*0 -f- 5*1 with 

Sgr= [ d^x^^{KR + L^,t} , (26a) 

5o = /" d^x^^ {af^R^ + 4°^i?^,i?^'' + af'R^.saR'"'^"} , (26b) 

Si= f d^xy^ ia[^'^^R^ + ai^'^^R^^R^'" + ai^^^'dR^.^s^R^"' 



ivba 



2 



otl^R^,i:R^-'" - ^ [V^W^i? + 2y(^)] \ . (26c) 



Here, iSgr is the Einstein-Hilbert plus matter action, while Sg and 5*1 are corrections. The 
former is decoupled from t9, where the omitted term proportional to <x^ does not affect the 
classical field equations since it is topological, i.e. it can be rewritten as the total 4-divergence 
of some 4-current. Similarly, if the Oi\ were chosen to reconstruct the Gauss-Bonnet invariant, 

(a^ , a2 J 03 ) = (1, —4, 1)q!gb, then this combination would also be topological and not affect the 
classical field equations. On the other hand, S\ is a modification to GR with a direct (non-minimal) 
couphng to 1?, such that as the field goes to zero, the modified theory reduces to GR. 

Another restriction one usually makes to simplify modified gravity theories is to neglect the a^- 
terms and only consider the S\ modification, the so-called restricted modified quadratic gravity. 
The c^l terms represent corrections that are non-dynamical and resemble that of a certain class of 



^Technically, Einstcin-Dilaton-Gauss-Bonnet gravity has a very particular set of coupling functions /i(i3) = 
/2('') = fai"^) oc e'''', where 7 is a constant. In most cases, however, one can expand about 71? <C 1, so that the 
functions become linear in the scalar field 
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/(i?) theories. As such, they can be mapped to scalar tensor theories with comphcated potentials, 
which can be heavily constrained by current experiments. The most popular such model is the one 
where only al is non-zero, but this has been greatly constrained by the torsion-balance Eot-Wash 
experiments to a]^ ' < 2 x 10~® m^ |208ll230ll57j . Moreover, these theories have a fixed coupling 
constant that does not run with energy or scale. In restricted modified gravity, the scalar field is 
effectively forcing the running of the coupling. 

Let us then concentrate on restricted modified quadratic gravity and drop the superscript in 
al . The modified field equations are 
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where we have defined 



Jjij =OR^i, + 2R^s„„R " 
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^dfifOR— y i^^R, 



SfiJ =8i? " Rf_iSua — '^QfiiyR "^ RSa + '^^^R^j.v ~ 2R Rf^^ + -gf^uR — 2V^^R , 

:kW ee - 4(V(^??)V,)i? - 2i?V^,t? + 5^, [2Rnd + 4(V^i?)V5i?] , 
JW EE - (V(^i?)V,)i? - 2V*i?V(^i?,)5 + 2V^^VsRt.u 
- 2Rs(^V'V,)d + g^, {V'^VsR + R^'^Vg.^) , 

a« EE - 8 (v^^) (V(,,i?,)5 - v^i?^.) + AR^su.y^^^ , 



-R^.Qi^ 



The t? stress-energy tensor is 



(V^^)(V,^) - -g^, (V,W^i9 - 2FW) 



The field equations for the scalar field are 
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aiR — a2R^„R^'' — a^R^^^saR^^" " — oiiRfj^i^sa* R'^'^ 






(27) 

(28a) 

(28b) 

(28c) 
(28d) 

(28e) 
(28f) 
(28g) 

(29) 
(30) 



Notice that unlike traditional scalar-tensor theories, the scalar field is here sourced by the geometry 
and not by the matter distribution. This directly implies that black holes in such theories are likely 
to be hairy. 

From the structure of the above equations, it should be clear that the dynamics of i? guarantee 
that the modified field equations are covariantly conserved exactly. That is, one can easily verify 



that the covariant divergence of Eq. (27) identically vanishes upon imposition of Eq. (30 1. Such 



a result had to be so, as the action is diffeomorphism invariant. If one neglected the kinetic and 
potential energies of ^ in the action, as was originally done in [216;, the theory would possess 
preferred-frame effects and will not be covariantly conserved. Moreover, such a theory requires an 



additional constraint, ie. the right-hand side of (30 1 would have to vanish, which is an unphysical 



consequence of treating § as prior structure ^431, 184 . 

One last simplification that is usually made when studying modified quadratic gravity theories 
is to ignore the potential V{-d), ie. set V{'d) = 0. This potential can in principle be non-zero, for 
example if one wishes to endow d with a mass or if one wishes to introduce a cosine driving term. 
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like that for axions in field theory. However, reasons exist to restrict the functional form of such a 
potential. First, a mass for t9 will modify the evolution of any gravitational degree of freedom only 
if this mass is comparable to the inverse length scale of the problem under consideration (such as 
a binary system). As one can see, this requires a significant amount of fine-tuning. Second, no 
mass term may appear in a theory with a shift symmetry, ie. invariance under the transformation 
'&—>'& + const. Such symmetries are common in four-dimensional, low-energy, effective string 
theories [751ll81[|180l[551[5^ . such as dynamical Chern-Simons and Einstein-Dilaton-Gauss-Bonnet 
theory. Similar considerations apply to other more complicated potentials, such as a cosine term. 
Given these field equations, one can linearize them about Minkowski space to find evolution 
equations for the perturbation in the small-coupling approximation. Doing so, one finds |408j 

n ,9 - - ^ {—\ t2 - ^ { — \ Tf"" T'"^' 

-'^-fr^^'Ks^.pK^''''^,, (31) 

where we have order-reduced the theory where possible and used the harmonic gauge condition 
(which is preserved in this class of theories |3531I361] ). The corresponding equation for the metric 
perturbation is rather lengthy and can be found in Eqs. (17)-(24) in |408j . Since these theories 
are to be considered effective, working always to leading order in a^, one can show that they are 
perturb at ively of type N2 in the E{2) classification |146j . ie. in the far zone, the only propagating 
modes that survive are the two transverse-traceless (spin-2) metric perturbations |353| . In the 
strong-field region, however, it is possible that additional modes are excited, although they decay 
rapidly as they propagate to future null infinity. 

Lastly, let us discuss what is known about whether modified quadratic gravities satisfy the 
requirements discussed in Sec. |2.1| As it should be clear from the action itself, this modified 
gravity theory satisfies the fundamental requirement, ie. passing all precision tests, provided the 
couplings ai are sufficiently small. This is because such theories have a continuous limit to General 
Relativity as a^ — >■ 0. Dynamical Chern-Simons gravity is constrained only weakly at the moment, 
^4 < 10^ km, where ^4 = a|/(/3K), only through observations of Lense-Thirring precession in the 
Solar System ^71. The Einstein-Dilaton-Gauss-Bonnet gravity coupling constant ^3 = a^/{/3K), 
on the other hand, has been constrained by several experiments: Solar System observations of the 
Shapiro time delay with the Cassini spacecraft placed the bound ^3 < 1.3 x 10^ km |66l[26]; the 
requirement that neutron stars still exist in this theory placed the constraint ^3 < 26km |307) . 
with the details depending somewhat on the central density of the neutron star; observations of 
the rate of change of the orbital period in the low-mass X-ray binary A0Q2Q — 00 |322l I226| has led 
to the constraint fg < 1.9 km |405j . 

Not all sub-properties of the fundamental requirement, however, are known to be satisfied. 
One can show that certain members of modified quadratic gravity possess known solutions and 
these are stable, at least in the small-coupling approximation. For example, in Dynamical Chern- 
Simons gravity, spherically symmetric vacuum solutions are given by the Schwarzschild metric 
with constant d to all orders in ai |216l 1431] . Moreover, one can show that such a solution is 
stable to small perturbations [2871 I170J , as well as non-spinning black holes and branes in anti De- 
Sitter space [129'. Spinning solutions, on the other hand, continue to be elusive, with approximate 
solutions in the slow-rotation/small-coupling limit known both for black holes |427l I241t 310, 4l5] 
and stars p!71l307j . On the other hand, in Einstein-Dilaton-Gauss-Bonnet theory even spherically 
symmetric solutions are modified |4341I310] . Recent work suggests, however, that these spherically 
symmetric solutions are stable to axial perturbations 
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The study of niodified quadratic gravity theories as effective theories is valid provided one 
is sufficiently far from its cut-off scale, ie. the scale beyond which higher-order curvature terms 
cannot be neglected anymore. One can estimate the magnitude of this scale by studying the size 
of loop corrections to the quadratic curvature terms in the action due to n-point interactions |416) . 
Simple counting requires that the number of scalar and graviton propagators, Pg and Pg , satisfy 
the following relation in terms of the number of vertices V: 



V V 



Ps^^, Pg = {n-1)-. (32) 



Loop corrections are thus suppressed by factors of aY M^^ " A"^, with Mpi the Planck mass 
and A the energy scale introduced by dimensional arguments. The cut-off scale above which the 
theory cannot be treated as an effective one can be approximated as the value of A at which the 
suppression factor becomes equal to unity: 

A,^A/pV'^"ay", (33) 

This cut-off scale automatically places a constraint on the magnitude of a^ above which higher- 
curvature corrections must be included. Setting the largest value of Ac to be equal to O(10/^m), 
thus saturating bounds from table-top experiments j230j . and solving for a^, we find 

a-/^ < O(lO^km). (34) 

Current Solar System bounds on a^ already require the coupling constant to be smaller than 
10^ km, thus justifying the treatment of these theories as effective models. 

As for the other requirements discussed in Sec. |2.1[ it is clear that modified quadratic gravity 
is well-motivated from fundamental theory, but it is not clear at all whether it has a well-posed 
initial value formulation. From an effective point of view, a perturbative treatment in a^ naturally 
leads to stable solutions and a well-posed initial value problem, but this is probably not the case 
when it is treated as an exact theory. In fact, if one were to treat such a theory as exact (to all 
orders in a^), then the evolution system would likely not be hyperbolic, as higher than second 
time derivatives now drive the evolution. Although no proof exists, it is likely that such an exact 
theory is not well-posed as an initial value problem. Notice, however, that this says nothing about 
the fundamental theories that modified quadratic gravity derives from. This is because even if 
the truncated theory were ill-posed, higher order corrections that are neglected in the truncated 
version could restore well-posedness. 

As for the last requirement (that the theory modifies the strong-field), modified quadratic 
theories are ideal in this respect. This is because they introduce corrections to the action that 
depend on higher powers of the curvature. In the strong-field, such higher powers could potentially 
become non-negligible relative to the Einstein-Hilbert action. Moreover, since the curvature scales 
inversely with the mass of the objects under consideration, one expects the largest deviations 
in systems with small total mass, such as stellar-mass black hole mergers. On the other hand, 
deviations from General Relativity should be small for small compact objects spiraling into a 
supermassive black hole, since here the spacetime curvature is dominated by the large object, and 
thus it is small, as discussed in |353) . 

2.3.4 Variable G Theories and Large Extra Dimensions 

Variable G theories are defined as those where Newton's gravitational constant is promoted to a 
spacetime function. Such a modification breaks the principle of equivalence (see |399j ) because 
the laws of physics now become local position dependent. In turn, this implies that experimental 
results now depend on the spacetime position of the laboratory frame at the time of the experiment. 
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Many known alternative theories that violate the principle of equivalence, and in particular, 
the strong equivalence principle, predict a varying gravitational constant. A classic example is 
scalar-tensor theory ' 396] . which, as explained in Sec. 2.3. 1[ modifies the gravitational sector of 



the action by multiplying the Ricci scalar by a scalar field (in the Jordan frame). In such theories, 
one can effectively think of the scalar as promoting the coupling between gravity and matter 
to a field-dependent quantity G — >■ G{(f>), thus violating local position invariance when varies. 
Another example are bimetric theories, such as that of Lightman-Lee [263j . where the gravitational 
constant becomes time-dependent even in the absence of matter, due to possibly time-dependent 
cosmological evolution of the prior geometry. A final example are higher-dimensional, brane- 
world scenarios, where enhanced Hawking radiation inexorably leads to a time-varying effective 
4D gravitational constant [126) . whose rate of change depends on the curvature radius of extra 
dimensions [226] . 

One can also construct /(i?)-type actions that introduce variability to Newton's constant. For 
example, consider the f{R) model [160] 



S — df^Xy/g kR 



1 + "0 In I ^- 



(35) 



where k — (IGttG) ^, ag is a coupling constant and Rq is a curvature scale. This action is motivated 
by certain renormalization group flow arguments |160| . The field equations are 



G.. - ^JT ~ f ^.^ - 2|||V[,i?V.ji? - ~9,.nR^ (36) 



where we have defined the new constant 






(37) 



Clearly, the new coupling constant k depends on the curvature scale involved in the problem, and 
thus, on the geometry, forcing G to run to zero in the ultraviolet limit. 

An important point to address is whether variable G theories can lead to modifications to a 
vacuum spacetime, such as a black hole binary inspiral. In Einstein's theory, G appears as the 
coupling constant between geometry, encoded by the Einstein tensor G^j/, and matter, encoded 
by the stress energy tensor T™**'. When considering vacuum spacetimes, T";f* = and one might 
naively conclude that a variable G would not introduce any modification to such spacetimes. In 
fact, this is the case in scalar-tensor theories (without homogeneous solutions to the scalar field 
equation), where the no- hair theorem establishes that black hole solutions are not modified. On the 
other hand, scalar-tensor theories with a homogenous scalar field solution can violate the no-hair 
theorem, endowing black holes with time-dependent hair, which in turn would introduce variability 
into G even in vacuum spacetimes |217l 1207) . 

In general, Newton's constant plays a much more fundamental role than merely a coupling 
constant: it defines the relationship between energy and length. For example, for the vacuum 
Schwarzschild solution, G establishes the relationship between the radius R of the black hole and 
the rest-mass energy E of the spacetime via R = 2GE/c'^. Similarly, in a binary black hole 
spacetime, each black hole introduces an energy scale into the problem that is quantified by a 
specification of Newton's constant. Therefore, one can treat variable G modifications as induced 
by some effective theory that modifies the mapping between the curvature scale and the energy 
scale of the problem, as is done for example in theories with extra dimensions. 

An explicit example of this idea is realized in braneworld models. Superstring theory suggests 
that physics should be described by 4 large dimensions, plus another 6 that are compactified and 
very small ^3181 I319J . The size of these extra dimensions is greatly constrained by particle theory 
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experiments. Braneworld models, where a certain higher-dimensional membrane is embedded in a 
higher dimensional bulk spacetime, can however evade this constraint as only gravitons can interact 
with the bulk. The ADD model [5^1 [50] is a particular example of such a braneworld, where the 
bulk is flat and compact and the brane is tensionless with ordinary fields localized on it. Since 
gravitational wave experiments have not yet constrained deviations from Einstein's theory in the 
strong-field, the size of these extra dimensions is constrained to micrometer scales only by table-top 
experiments [2301 \5\ . What is relevant to gravitational wave experiments is that in many of these 
braneworld models black holes may not remain static |1481 1366] . The argument goes roughly as 
follows: a five-dimensional black hole is dual to a four-dimensional one with conformal fields on 
it by the ADS/CFT conjecture [2701 17]. but since the latter must evolve via Hawking radiation, 
the black hole must be losing mass. The Hawking mass loss rate is here enhanced by the large 
number of degrees of freedom in the conformal field theory, leading to an effective modification 
to Newton's laws and to the emission of gravitational radiation. Effectively, one can think of the 
black hole mass loss as due to the black hole being stretched away from the brane into the bulk 
due to an universal acceleration, that essentially reduces the size of the brane-localized black hole. 
For black-hole binaries, one can then draw an analogy between this induced time-dependence in 
the black hole mass and a variable G theory, where Newton's constant decays due to the presence 
of black holes. Of course, this is only analogy, since large extra dimensions would not predict a 
time-evolving mass in neutron star binaries. 

Many experiments have been carried out to measure possible deviations from a constant G value, 
and they can broadly be classified into two groups: (a) those that search for the present or nearly 
present rate of variation (at redshifts close to zero); (b) those that search for secular variations over 
long time periods (at very large redshifts) . Examples of experiments or observations of the first class 
include planetary radar-ranging f314] , surface temperature observations of low-redshift millisecond 
pulsars . 2201 1326) . lunar ranging observations ^403) and pulsar timing observations |231l 1128) . 
the latter two being the most stringent. Examples of experiments of the second class include 
the evolution of the Sun [185] and Big Bang Nucleosynthesis (BBN) calculations jll2) |42] , again 
with the latter being more stringent. For either class, the strongest constraint are about G/G < 
10~^^ yr~^, varying somewhat from experiment to experiment. 

Lacking a particularly compelling action to describe variable G theories, one is usually left 
with a phenomenological model of how such a modification to Einstein's theory would impact 
gravitational waves. Given that the part of the waveform that detectors are most sensitive to is 
the gravitational wave phase, one can model the effect of variable G theories by studying how the 
rate of change of its frequency would be modified. Assuming a Taylor expansion for Newton's 
constant one can derive the modification to the evolution equation for the gravitational wave 
frequency, given whichever physical scenario one is considering. Solving such evolution equation 
then leads to a modification in the accumulated gravitational wave phase observed at detectors on 
Earth. In Sec. [5] we will provide an explicit example of this for a compact binary system. 

Let us now discuss whether such theories satisfy the criteria defined in Sec. |2.1[ The funda- 
mental property can be satisfied if the rate of change of Newton's constant is small enough, as 
variable G theories usually have a continuous limit to General Relativity (as all derivatives of G go 
to zero). Whether variable G theories are well- motivated from fundamental physics (property 2) 
depends somewhat on the particular effective model or action that one considers. But in general, 
property 2 is usually satisfied, considering that such variability naturally arises in theories with 
extra dimensions, and the latter are also natural in all string theories. Variable G theories, how- 
ever, usually fail at introducing modifications in the strong-field region. Usually, such variability 
is parameterized as a Taylor expansion about some initial point with constant coefficients. That 
is, the variability of G is not usually constructed so as to become stronger closer to merger. 

The well-posed property and the sub-properties of the fundamental property depend somewhat 
on the particular effective theory used to describe varying G modifications. In the f{R) case, one 
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can impose restrictions on the functional form /(•) such that no ghosts (/' > 0) or instabihties 
(/" > 0) arise |160j . This, of course, does not guarantee that this (or any other such) theory is 
weU-posed. A much more detailed analysis would be required to prove well-posedness of the class 
of theories that lead to a variable Newton's constant, but such is currently lacking. 

2.3.5 Non-Commutative Geometry 

Non-commutative geometry is a gravitational theory that generalizes the continuum Riemannian 
manifold of Einstein's theory with the product of it with a tiny, discrete, finite non-commutative 
space, composed of only two points. Although the non-commutative space has zero spacetime 
dimension, as the product manifold remains four dimensional, its internal dimensions are 6 to 
account for Weyl and chiral fermions. This space is discrete to avoid the infinite tower of massive 
particles that would otherwise be generated, as in string theory. Through this construction, one 
can recover the Standard Model of elementary particles, while accounting for all (elementary 
particle) experimental data to date. Of course, the simple non-commutative space described above 
is expected to be replaced by a more complex model at Planckian energies. Thus, one is expected to 
treat such non-commutative geometry models as effective theories. Essentially nothing is currently 
known about the full non-commutative theory of which the theories described in this section are 
an effective low-energy limit. 

Before proceeding with an action-principle description of non-commutative geometry theories, 
we must distinguish between the spectral geometry approach championed by Connes |107j . and 
Moyal-type non-commutative geometries |352(ll83]l289j . In the former, the manifold is promoted 
to a non- commutative object through the product of a Riemann manifold with a non-commutative 
space. In the latter, instead, a non-trivial set of commutation relations is imposed between oper- 
ators corresponding to position. These two theories are in principle unrelated. In this review, we 
will concentrate only in the former, as it is the only type of non-commutative General Relativity 
extension that has been studied in the context of gravitational wave theory. 

The effective action for spectral non-commutative geometry theories (henceforth, non-commutative 
geometries for short) is 



S= f d^x^^ {kR + aoC^.saC^'"^^ + roR*R* - ^qR \H\^) + 5„ 



(38) 



where H is related to the Higgs field, C^j^uSa is the Weyl tensor, (ao, tq, ^o) are couplings constants 
and we have defined the quantity 

R*R* :- ^e'^-^'ea^^^i^M-^^iip.^' • (39) 

Notice that this term integrates to the Euler characteristic, and since tq is a constant, it is topolog- 



ical and does not affect the classical field equations. The last term of Eq. (38 1 is usually ignored as 
H is assumed to be relevant only in the early universe. Finally, the second term can be rewritten 
in terms of the Riemann and Ricci tensors as 



C^.^SaC^"^" = Tji?' - '^Rt^.R'"' + Rf.uSaR'"'^^ ■ (40) 



Notice that this corresponds to the modified quadratic gravity action of Eq. ( 26 1 with all al — 



„(i) 



and {al ,03 ,0;)^ ) = (1/3, —2, 1) , which is not the Gauss-Bonnet invariant. Notice also that this 
model is not usually studied in modified quadratic gravity theory, as one usually concentrates on 
the terms that have an explicit scalar field coupling. 
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The field equations of this theory can be read directly from Eq. (27), but we repeat them here 
for completeness: 

G^, - ^ [2^^^ + R^'^] C^A.K = :^t;Z''' ■ (41) 

One could in principle rewrite this in terms of the Riemann and Ricci tensors, but the expressions 
become quite complicated, as calculated explicitly in Eqs. (2) and (3) of j434j . Due to the absence 
of a dynamical degree of freedom coupling to the modifications to the Einstein-Hilbert action, this 
theory is not covariantly conserved in vacuum. By this we mean that the covariant divergence of 



Eq. (41 1 does not vanish in vacuum, thus violating the weak-equivalence principle and leading to 
additional equations that might over-constrain the system. In the presence of matter, the equations 
of motion will not be given by the vanishing of the covariant divergence of the matter stress-energy 
alone, but now there will be additional geometric terms. 

Given these field equations, one can linearize them about a flat background to find the evolution 
equations for the metric perturbation |2931 1292] 

(i-rH)°'7V = -i6^r;i"'*' (42) 

where the term proportional to /3^ — (— 327rao)^^ acts like a mass term. Here, one has imposed the 
transverse-traceless gauge (a refinement of Lorenz gauge) , which can be shown to exist [2931 I292J . 
Clearly, even though the full non-linear equations are not covariantly conserved, its linearized 



version is, as one can easily show that the divergence of the left-hand side of Eq. (42) vanishes. 
Because of these features, if one works perturbative in I3~^ , then such a theory will only possess 
the two usual transverse-traceless (spin-2) polarization modes, ie. it is perturbatively of type N2 
in the E{2) classification 146 . 

Let us now discuss whether such a theory satisfies the properties discussed in Sec. |2.1| Non- 
commutative geometry theories clearly possess the fundamental property, as one can always take 
ao — >■ (or equivalently /3^^ -^ 0) to recover General Relativity. Therefore, there must exist 
a sufficiently small ao such that all precision tests carried out to date are satisfied. As for the 
existence and stability of known solutions, Refs. [2931 I292J have shown that Minkowski spacetime 
is stable only for a^ < Q, as otherwise a tachyonic term appears in the evolution of the metric 



perturbation, as can be seen from Eq. (42). This then automatically implies that j3 must be real. 



Current constraints on Weyl terms of this form come mostly from Solar System experiments. 



Ni [295 recently studied an action of the form of Eq. ( 38 ) minimally coupled to matter in light of 
Solar System experiments. He calculated the relativistic Shapiro time-delay and light deflection 
about a massive body in this theory and found that observations of the Cassini satellite place 
constraints on jaop ^ < 5.7 km [295) . This is currently the strongest bound we are aware of on uq. 
Many solutions of General Relativity are preserved in non-commutative geometries. Regarding 
black holes, all solutions that are Ricci flat (vacuum solutions of the Einstein equations) are also 



solutions to Eq. (41 ). This is because by the second Bianchi identity, one can show that 

V"^i?^A.« = V^i?,,« - UR^, , (43) 



and the right-hand side vanishes in vacuum, forcing the entire left-hand side of Eq. (41 ) to vanish. 
This is not so for neutron stars, however, where the equations of motion are likely to be modified, 
unless they are static |291j . Moreover, as of now there has been no stability analysis of black hole 
or stellar solutions and no study of whether the theory is well-posed as an initial value problem, 
even as an effective theory. Thus, except for the fundamental property, it is not clear that non- 
commutative geometries satisfy any of the other criteria listed in Sec. |2.1[ 
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2.3.6 Gravitational Parity Violation 

Parity, the symmetry transformation that flips the sign of the spatial triad, has been found to 
be broken in the Standard Model of elementary interactions. Only the combination of a parity 
transformation, time inversion and charge conjugation (CPT) remains still a true symmetry of the 
Standard Model. Experimentally, it is curious that the weak interaction exhibits maximal parity 
violation, while other fundamental forces seem to not exhibit any. Theoretically, parity violation 
unavoidingly arises in the Standard Model [SHI [H [E] , as there exists one-loop chiral anomalies 
that give rise to parity violating terms coupled to lepton number |389) . In certain sectors of 
string theory, such as in heterotic and in Type I superstring theories, parity violation terms are 
also generated through the Green-Schwarz gauge anomaly-canceling mechanism |18H 13191 [15]. 
Finally, in loop quantum gravity |37| , the scalarization of the Barbero-Immirzi parameter coupled 
to fermions leads to an effective action that contains parity- violating terms |367l [551 12791 1172) . 
Even without a particular theoretical model, one can generically show that effective field theories 
of inflation generically contain non-vanishing, second-order, parity violating curvature corrections 
to the Einstein-Hilbert action |390) . Alternatively, phenomenological parity- violating extensions 
of General Relativity have been proposed through a scalarization of the fundamental constants of 
Nature [TDEl . 

One is then naturally led to ask whether the gravitational interaction is parity invariant in 
the strong field. A violation of parity invariance would occur if the Einstein-Hilbert action were 
modified through a term that involved a Levi-Civita tensor and parity invariant tensors or scalars. 
Let us try to construct such terms with only single powers of the Riemann tensor and a single 
scalar field d: 

(ia) R^p^s e"'"'', (ib) R^f,^^ e^^^" V^.z?, 

(ic) R^p^^ t^^'- V^.^i? , (id) i?„^^^ t^f''- V'"^0,/z? . 

Option (ia) and (ib) vanish by the Bianchi identities. Options (ic) and (id) include the commutator 
of covariant derivatives, which can be rewritten in terms of a Riemann tensor, and thus it leads to 
terms that are at least quadratic in the Riemann tensor. Therefore, no scalar can be constructed 
that includes contractions with the Levi-Civita tensor from a single Riemann curvature tensor and 
a single field. One can try to construct a scalar from the Ricci tensor 

(iia) R^fi e"^^*V^5i?, (lib) i?„/3 e^^^^^^V ^s/d, 

but again (iia) vanishes by the symmetries of the Ricci tensor, while (lib) involves the commutator 
of covariant derivatives, which introduces another power of the curvature tensor. Obviously, the 
only term one can write with the Ricci scalar would lead to a double commutator of covariant 
derivatives, leading to extra factors of the curvature tensor. 

One is then forced to consider either theories with two mutually independent fields or theories 
with quadratic curvature tensors. Of the latter, the only combination that can be constructed and 
that does not vanish by the Bianchi identities is the so called Pontryagin density, ie. WR^ and 
therefore, the action |216l[n] 

" j4. ^ — / ^ a 



S= d^Xy/^UR+-dR*R] , (44) 

is the most general, quadratic action with a single scalar field that violates parity invariance, where 
we have rescaled the a prefactor to follow historical conventions. This action defines non-dynamical 
Chern-Simons modified gravity, initially proposed by Jackiw and Pi |216[ [14] . Notice that this is 



the same as the term proportional to a^ in the quadratic gravity action of Eq. ( 26 ) , except that here 



1? is prior geometry, ie. it does not possess self-consistent dynamics or an evolution equation. Such 
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a term violates parity invariance because the Pontryagin is a pseudo-scalar, while i? is assumed to 
be a scalar. 

The field equations for this theory arq^ 

G^. + £acW^i-T;r, (45) 



which is simply Eq. (27) with (ai, 02,013) set to zero and no stress-energy for •&. Clearly, these 
field equations are not covariantly conserved in vacuum, ie. taking the covariant divergence one 
finds the constraint 

aR*R^O. (46) 

This constraint restricts the space of allowed solutions, for example disallowing the Kerr met- 
ric [184i. Therefore, it might seem that the evolution equations for the metric are now overcon- 
strained, given that the field equations provide 10 differential conditions for the 10 independent 
components of the metric tensor, while the constraint adds one additional, independent differential 



condition. Moreover, unless the Pontryagin constraint, Eq. (46), is satisfied, matter fields will not 
evolve according to V^T^^* — 0, thus violating the equivalence principle. 

From the field equations, we can derive an evolution equation for the metric perturbation when 
linearizing about a flat background, namely 

hh Hi 

in a transverse-traceless gauge, which can be shown to exist in this theory [9l I422J . The constraint 



of Eq. ( 46 ) is here identically satisfied to second order in the metric perturbation. Without further 



information about d, however, one cannot proceed any further, except for a few general observa- 



tions. As is clear from Eq. (47), the evolution equation for the metric perturbation can contain 
third time derivatives, which generically will lead to instabilities. In fact, as shown in |10j the 
general solution to these equations will contain exponentially growing and decaying modes. The 



theory defined by Eq. (44), however, is an effective theory, and thus, there can be higher order 
operators not included in this action that may stabilize the solution. Regardless, when studying 
this theory order-reduction is necessary if one is to consider it an effective model. 

Let us now discuss the properties of such an effective theory. Because of the structure of the 
modification to the field equations, one can always choose a sufficiently small value for a such 
that all Solar System tests are satisfied. In fact, one can see from the equations in this section 
that in the limit a — >■ 0, one recovers General Relativity. Non-dynamical Chern-Simons gravity 
leads to modifications to the non-radiative (near-zone) metric in the so-called gravitomagnetic 
sector, leading to corrections to Lense-Thirring precession [11] [T^]. This fact has been used to 
constrain the theory through observations of the orbital motion of the LAGEOS satellites |351j to 
(«/«)■(? < 2 X 10^ km, or equivalently {K/a)'d~^ > 10^^'* eV. Much better constraints, however, 
can be placed through observations of the double binary pulsar [4331 [IB] : {a4/K)'d < 0.8 km. 

Some of the sub-properties of the fundamental requirement are satisfied in non-Dynamical 
Chern-Simons Gravity. On the one hand, all spherically symmetric metrics that are solutions to 
the Einstein equations are also solutions in this theory for a "canonical" scalar field {9 (x t) jl84j . 
On the other hand, axisymmetric solutions to the Einstein equations are generically not solutions 
in this theory. Moreover, although spherically symmetric solutions are preserved, perturbations of 
such spacetimes that are solutions to the Einstein equations are not generically solutions to the 
modified theory |431| . What is perhaps worse, the evolution of perturbations to non-spinning black 
holes have been found to be generically overconstrained [431' . This is a consequence of the lack of 



^The tensor 3CpJ is sometimes written as C^u and referred to as the C-tensor. 
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scalar field dynamics in the modified theory, which via Eq. (461 tends to overconstrain it. Such a 
conclusion also suggests that this theory does not posses a well-posed initial value problem. 

One can argue that Non-Dynamical Chern-Simons Gravity is well-motivated from fundamental 
theories [l3j, except that in the latter, the scalar field is always dynamical, instead of having to be 
prescribed a priori. Thus, perhaps the strongest motivation for such a model is as a phenomeno- 
logical proxy to test whether the gravitational interaction remains parity invariant in the strong 
field, a test that is uniquely suited to this modified model. 

2.4 Currently Unexplored Theories in the Gravitational Wave Sector 

The list of theories we have here described is by no means exhaustive. In fact, there are many 
fascinating theories that we have chosen to leave out because they have not yet been analyzed in 
the gravitational wave context in detail. Examples of these include the following: 

• Einstein- Aether Theory [218j and Horava-Lifshitz Theory |205j : 

• Standard Model Extension [T02] : 

• Eddington- inspired Born-Infeld gravity [43]; 

• New Massive Gravity [Ml EH] and Bi-Gravity Theories [3T31I3TT] : 

We will update this review with a description of these theories, once a detailed gravitational 
wave study for compact binaries or supernovae sources is carried out and the predictions for 
the gravitational waveform observables are made for any physical system plausibly detectable by 
current or near future gravitational wave experiments. 
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3 Detectors 

3.1 Gravitational Wave Interferometers 

Kilometer scale gravitational wave interferometers have been in operation for over a decade. This 
type of detectors use laser interferometry to monitor the locations of test masses at the ends of 
the arms with exquisite precision. Gravitational waves change the relative length of the optical 
cavities in the interferometer (or equivalently, the proper travel time of photons) resulting in a 
strain 

where AL ~ L^ ^ Ly is the path length difference between the two arms of the interferometer. 

Fractional changes in the difference in path lengths along the two arms can be monitored to 
better than 1 part in 10^'^. It is not hard to understand how such precision can be achieved. For a 
simple Michelson interferometer, a difference in path length of order the size of a fringe can easily 
be detected. For the typically-used, infrared lasers of wavelength A ^ 1 /im, and interferometer 
arms of length L = 4 km, the minimum detectable strain is 

This is still far off the 10^^° mark. In principle, however, changes in the length of the cavities 
corresponding to fractions of a single fringe can also be measured provided we have a sensitive 
photodiode at the dark port of the interferometer, and enough photons to perform the measure- 
ment. This way we can track changes in the amount of light incident on the photodiode as the 
lengths of the arms change and we move over a fringe. The rate at which photons arrive at the 
photodiode is a Poisson process and the fluctuations in the number of photons is ^ iV^/^, where 
N is the number of photons. Therefore we can track changes in the path length difference of order 

AX ^ 



iVl/2- 



The number of photons depends on the laser power P, and the amount of time available to 
perform the measurement. For a gravitational wave of frequency /, we can collect photons for a 
time t ^ 1//, so the number of photons is 

jhpV 

where hp is Planck's constant and u = c/X is the laser frequency. For a typical laser power P ^ 1 W, 
a gravitational wave frequency / = 100 Hz, and A ^ 1 /im the number of photons is 

iV- 10^^, 

so that the strain we are sensitive to becomes 

The sensitivity can be further improved by increasing the effective length of the arms. In the 
LIGO instruments, for example, each of the two arms forms a resonant Fabry-Perot cavity. For 
gravitational wave frequencies smaller than the inverse of the light storage time, the light in the 
cavities makes many back and forth trips in the arms while the wave is traversing the instrument. 
For gravitational waves of frequencies around 100 Hz and below, the light makes about a thousand 
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back and forth trips while the gravitational wave is traversing the interferometer, which results in 
a 3 order of magnitude improvement in sensitivity, 

/i- 10^21. 

For frequencies larger than 100 Hz the number of round trips the light makes in the Fabry-Perot 
cavities while the gravitational wave is traversing the instrument is reduced and the sensitivity is 
degraded. 

The proper light travel time of photons in interferometers is controlled by the metric pertur- 
bation, which can be expressed as a sum over polarization modes 

/i,,(t,f)=^;i^.(t,x), (48) 

A 

where A, labels the six possible polarization modes in metric theories of gravity. The metric 
perturbation for each mode can be written in terms of a plane wave expansion, 

ht,{t,x) = |_°° df j ^ dne^'^f^'-^-'yh^{f,(l)et^{Cl). (49) 

Here / is the frequency of the gravitational waves, k — 27r/J7 is the wave vector, f2 is a unit vector 
that points in the direction of propagation of the gravitational waves, e^ is the Ath polarization 
tensor, with i, j — x,y, z spatial indices. The metric perturbation due to mode A from the direction 
J7 can be written by integrating over all frequencies, 

/•OO 

hf^{t-Cl-x)= df e^^^f^'-''-'^h^{f, Cl)ef^{Cl). (50) 



By integrating Eq. ( 49 ) over all frequencies we have an expression for the metric perturbation from 
a particular direction fl, i.e. only a function oi t — fl ■ x. The full metric perturbation due to a 
gravitational wave from a direction fl can be written as a sum over all polarization modes 

%(i -n-x)=J2 h'^i^ - ^ ■ ^yfji^)- (51) 

A 

The response of an interferometer to gravitational waves is generally referred to as the antenna 
pattern response, and depends on the geometry of the detector and the direction and polarization 
of the gravitational wave. To derive the antenna pattern response of an interferometer for all six 
polarization modes we follow the discussion in |296j closely. For a gravitational wave propagating 
in the z direction, the polarization tensors are as follows 

4- = 




4 = I 1 I ,e^, = I I , (52) 

where the superscripts +, x, cc, y, 5, and £ correspond to the plus, cross, vector-x, vector-y, 
breathing, and longitudinal modes. 
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Figure 1: Detector coordinate system and gravitational wave coordinate system. 



Suppose that the coordinate system for the detector is i; = (1,0,0),?; = (0, l,0),i = (0,0,1), 
as in Fig.[T] Relative to the detector, the gravitational wave coordinate system is rotated by angles 
{0,(j)), x' = {cos 6 cos (j), cos 9 sin (j), — sin 9) , y' = (— sine/), cos 0, 0), and z' — (sin 6* cos 0, sin sin </>, cos ( 
We still have the freedom to perform a rotation about the gravitational wave propagation direction 
which introduces the polarization angle V', 

rh ^ x' cos '0 + y' sin tp , 
h — ~x' sin Tp + y' cos tp , 

n^z' . 

The coordinate systems {x,y,z) and {m,n,Cl) are also shown in Fig. Ill To generalize the polar- 
ization tensors in Eq. ( |52[ ) to a wave coming from a direction Ct, we use the unit vectors m, n, and 
Cl as follows 



e ' 


= m®m. — n®n, 


e^ 


— rh® fi + h® rh, 


e^ 


= m ® Vl + Cl ® m, 


e^ 


= n®Cl + (l®n, 


e" 


= m <S) m + n <® ji, 



= Q®n. 



(53) 



For LIGO and VIRGO the arms are perpendicular so that the antenna pattern response can be 
written as the difference of projection of the polarization tensor onto each of the interferometer 
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arms, 

This means that the strain measured by an interferometer due to a gravitational wave from direction 
fl and polarization angle tp takes the form 

h{t)=J2f^A{t-Cl-x)F'^{n,i;). (54) 

A 

Explicitly, the antenna pattern functions are, 

F^{9,(l),^) = -(1 + cos^0)cos20cos2'0 - cos0sin2(/)sin2V', 

F^{e,(f>,^) ^ --(1 + cos^ 6*) cos 20 sin 27/1- cos 6* sin 20 cos 2-0, 

F^{6,(f>,4>) — sin 6* (cos cos 20 cos — sin 20 sin 0;) , 
F^{6,(f>,^) — — sin0 (cos0cos20sin0j + sin20cos0'), 

F^e^cj)) = -isin2 6lcos20, 

i^^(6',0) = ^sin2 6lcos20. (55) 

The dependence on the polarization angles "0, reveals that the + and x polarizations are spin-2 
tensor modes, the x and y polarizations are spin-1 vector modes, and the b and £ polarizations 
are spin-0 scalar modes. Note that for interferometers the antenna pattern responses of the scalar 
modes are degenerate. Figure [2] shows the antenna patterns for the various polarizations given in 



Eq. (55 1 with = 0. The color indicates the strength of the response with red being the strongest 



and blue being the weakest. 

3.2 Pulsar Timing Arrays 

Neutron stars can emit powerful beams of radio waves from their magnetic poles. If the rotational 
and magnetic axes are not aligned, the beams sweep through space like the beacon on a lighthouse. 
If the line of sight is aligned with the magnetic axis at any point during the neutron star's rotation 
the star is observed as a source of periodic radio wave bursts. Such a neutron star is referred to as a 
pulsar. Due to their large moment of inertia pulsars are very stable rotators, and their radio pulses 
arrive on Earth with extraordinary regularity. Pulsar timing experiments exploit this regularity: 
gravitational waves are expected to cause fluctuations in the time of arrival of radio pulses from 
pulsars. 

The effect of a gravitational wave on the pulses propagating from a pulsar to Earth was first 
computed in the late 1970s by Sazhin and Detweiler [3411 1130| . Gravitational waves induce a 
redshift in the pulse train 

1 p^p^ 

2l + n.p' 



z{t,n) = - ^ ^^ ^ A%, (56) 



where p is a unit vector that points in the direction of the pulsar, Q is the unit vector of gravitational 
wave propagation, and 

Ahij = /i,y (ip, Cl) - h,j{t^, Q), (57) 

is the difference in the metric perturbation at the pulsar when the pulse was emitted and the metric 



perturbation on Earth when the pulse was received. The inner product in Eq. ( 56 ) is computed 
with the Euclidean metric. 
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(b) 




(c) 



(d) 




(e) 



Figure 2: Antenna pattern response functions of an interferometer (see Eqs. ( [55| )) for ip = 0. Panels 
(a) and (b) show the plus (|-F+|) and cross (|Fx |) modes, panels (c) and (d) the vector x and vector 
y modes {\Fx\ and \Fy\), and panel (e) shows the scalar modes (up to a sign, it is the same for 
both breathing and longitudinal). Color indicates the strength of the response with red being the 
strongest and blue being the weakest. The black lines near the center give the orientation of the 
interferometer arms. 



In pulsar timing experiments it is not the redshift, but rather the timing residual that is 
measured. The times of arrival of pulses are measured and the timing residual is produced by sub- 
tracting off a model that includes the rotational frequency of the pulsar, the spin-down (frequency 
derivative), binary parameters if the pulsar is in a binary, sky location and proper motion, etc. 
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The timing residual induced by a gravitational wave, R{t), is just the integral of the redshift 

R{t)= f dt' z{t'). (58) 



Times of arrival (TOAs) are measured a few times a year over the course of several years allowing 
for gravitational waves in the nano-Hertz band to be probed. Currently, the best timed pulsars 
have residual RMSs of a few 10s of ns over a few years. 

The equations above can be used to estimate the strain sensitivity of pulsar timing experiments. 
For gravitational waves of frequency / the expected induced residual is 

"^ 7' 

so that for pulsars with RMS residuals R ^ 100 ns, and gravitational waves of frequency / ^ 
10^® Hz, gravitational waves with strains 



Rf ^ 10 



-15 



would produce a measurable effect. 

To find the antenna pattern response of the pulsar-Earth system, we are free to place the pulsar 
on the z-axis. The response to gravitational waves of different polarizations can then be written 
as 



which allows us to express the Fourier transform of ( |56[ ) as 

z(/, n) - (e-2-/i(i+^p) - l) ^ HaU, (l)F^{h) (59) 

A 

where the sum is over all possible GW polarizations: A = +, x,x,y, 6, /, and L is the distance to 
the pulsar. 
Explicitly, 

f\ 

F+{e,ij) = sin2-cos2V' (60) 

n 

F''{e,ij) = - sin^ - sin 2i/; (61) 

F^9,^) = -^tIS^cos^, (62) 

^^(^'^) = ^ifS^^^^^^' (^^) 

F''{9) = sin^^ (64) 

fH9) = - ^°''^ . (65) 

^ ^ 2 1 + COS6' ^ ' 

Just like for the interferometer case, the dependence on the polarization angle ip, reveals that the 
+ and X polarizations are spin-2 tensor modes, the x and y polarizations are spin-1 vector modes, 
and the b and £ polarizations are spin-0 scalar modes. Unlike interferometers, the antenna pattern 
responses of the pulsar-Earth system do not depend on the azimuthal angle of the gravitational 
wave, and the scalar modes are not degenerate. 
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Figure 3: Antenna patterns for the pulsar-Earth system. The plus mode is shown in (a), breathing 
modes in (b), the vector-x mode in (c), and longitudinal modes in (d), as computed from Eq. (66 1. 
The cross mode and the vector-y mode are rotated versions of the plus mode and the vector-x mode, 
respectively, so we did not include them here. The GW propagates in the positive z-direction with 
the Earth at the origin, and the antenna pattern depends on the pulsar's location. The color 
indicates the strength of the response, red being the largest and blue the smallest. 



In the literature, it is common to write the antenna pattern response by fixing the GW direction 
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and changing the location of the pulsar. In this case the antenna pattern responses are |254[|1211[52] . 

F+iOp, (t)p) = sin^ ^ cos 2(j)p 

F''{ep,(f>p) = sm^^sm2(f>p 

i^xm J. \ 1 sin26'p 

F\ep) = sin^l 

~g _ 1 cos^ dp 

^^^^^ - 21 + cos9p' (^^) 

where 9p and (jjp are the polar and azimuthal angles of the vector pointing to the pulsar. Up to 



signs, these expressions are the same as Eq. (65 1 taking 9 ^ Op and ip ^ 'f'p- This is because fixing 
the GW propagation direction while allowing the pulsar location to change is analogous to fixing 
the pulsar position while allowing the direction of GW propagation to change-there is degeneracy 
in the GW polarization angle and the pulsar's azimuthal angle 0p. For example, changing the 
polarization angle of a gravitational wave traveling in the z-direction is the same as performing 
a rotation about the z-axis that changes the pulsar's azimuthal angle. Antenna patterns for the 
pulsar-Earth system using Eqs. (66) are shown in Fig. [SJ The color indicates the strength of the 



response, red being the largest and blue the smallest. 
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4 Testing Techniques 

4.1 Coalescence Analysis 

Gravitational waves emitted during the inspiral, merger and ringdown of compact binaries are the 
most studied in the context of data analysis and parameter estimation. In this section, we will 
review some of the main data analysis techniques employed in the context of parameter estimation 
and tests of General Relativity. We begin with a discussion of matched filtering and Fisher theory 
(for a detailed review, see |1551 15^ [1 18111561 1219| ). We then continue with a discussion of Bayesian 
parameter estimation and hypothesis testing (for a detailed review, see [3501 1182[ 11171 I^M] ). 

4.1.1 Matched Filtering and Fisher Analysis 

When the detector noise n{t) is Gaussian and stationary, and when the signal s(t) is known very 
well, the optimal detection strategy is matched filtering. For any given realization, such noise can 
be characterized by its power spectral density S'„(/), defined via 

{n{f)n*{n)^lsM)Hf-f'), (67) 

where recall that the overhead tilde stands for the Fourier transform, the overhead star stands for 
complex conjugation and the brackets for the expectation value. 

The detectability of a signal is determined by its signal-to-noise ratio or SNR, which is defined 
via 

P' = , ^''^^ , (68) 

where /i is a template with parameters A* and we have defined the inner product 

{A\B)^m\ -—d!. (69) 

Jo *5n 

If the templates do not exactly match the signal, then the signal-to-noise ratio is reduced by a 
factor of M, called the match: 

M^-id^L, (70) 

where 1 — M = MM is the mismatch. 

For the noise assumptions made here, the probability of measuring s(f) in the detector output, 
given a template /i, is given by 

pcxe-(^"''l^-'''/2, (71) 

and thus the waveform h that best fits the signal is that with best-fit parameters such that the 
argument of the exponential is minimized. For large signal-to-noise ratio, the best-fit parameters 
will have a multivariate Gaussian distribution centered on the true values of the signal A', and 
thus, the waveform parameters that best fit the signal minimize the argument of the exponential. 
The parameter errors (5 A' will be distributed according to 



where Yij is the Fisher matrix 



i{6X') ex e-3 


VijSysx^ 




dh\ 
dXi) 



(72) 
(73) 
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The root-mean-squared (Icr) error on a given parameter A' is then 

^{{5\^f)^^W\ (74) 



where S*^ = (Fy ) is the variance-covariance matrix and summation is not iniphcd in Eq. ( 74 ) 



(A' denotes a particular element of the vector A*). This root-mean-squared error is sometimes 



referred to as the statistical error in the measurement of A' . One can use Eq. ( 74 1 to estimate how 



well modified gravity parameters can be measured. Put another way, if a gravitational wave were 



detected and found consistent with General Relativity, Eq. ( 74 ) would provide an estimate of how 
close to zero these modified gravity parameters would have to be. 

The Fisher method to estimate projected constraints on modified gravity theory parameters 
is as follows. First, one constructs a waveform model in the particular modified gravity theory 
one wishes to constrain. Usually, this waveform will be similar to the General Relativity one, 
but it will contain an additional parameter, k, such that the templates parameters are now A* 
plus n. Let us assume that as k — >■ 0, the modified gravity waveform reduces to the General 
Relativity expectation. Then, the accuracy to which k, can be measured, or the accuracy to which 
we can say n is zero, is approximately (E'^")^/^, where the Fisher matrix associated with this 
variance-covariance matrix must be computed with the non-General Relativity model evaluated 
at the General Relativity limit (k — >■ 0). Such a method for estimating how well modified gravity 



theories can be constrained was pioneered by Will in |397[ I317J , and since then, it has been widely 
employed as a first-cut estimate of the accuracy to which different theories can be constrained. 

The Fisher method described above can dangerously lead to incorrect results if abused [3741 
1375] . One must understand that this method is suitable only if the noise is stationary and Gaussian 
and if the signal-to-noise ratio is sufficiently large. How large a signal-to-noise ratio is required for 
Fisher methods to work depends somewhat on the signals considered, but usually for applications 
concerning tests of General Relativity, one would be safe with p > 30 or so. In real data analysis, 
the first two conditions are almost never satisfied. Moreover, the first detections that will be made 
will probably be of low signal-to-noise ratio, ie. p ~ 8, for which again the Fisher method fails. In 
such cases, more sophisticated parameter estimation methods need to be employed. 

4.1.2 Bayesian Theory and Model Testing 

Bayesian theory is ideal for parameter estimation and model selection. Let us then assume that a 
detection has been made and that the gravitational wave signal in the data can be described by 
some model M, parameterized by the vector A*. Using Bayes' theorem, the posterior distribution 
function (PDF) or the probability density function for the model parameters, given data d and 
model M, is 

Obviously, the global maximum of the PDF in the parameter manifold gives the best fit parameters 
for that model. The prior probability density p(A*|M) represents our prior beliefs of the parameter 
range in model M. The marginalized likelihood or evidence^ is the normalization constant 

p{d\M) = / dX'p{d\X\M) p{X'\M) , (76) 



which clearly guarantees that the integral of Eq. (75 1 integrates to unity. The quantity p(d|A*, M) 



is the likelihood function, which is simply given by Eq. (71), with a given normalization. In 
that equation we used slightly different notation, with s being the data d and h the template 
associated with model M and parameterized by A*. The marginalized PDF, which represents the 
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probability density function for a given parameter A* (recall that A* is a particular element of A' 
after marginalizing over all other parameters, is given by 



p(A*|d,M)=/ dX'p{X'\M.)p{d\X\M) , (77) 

where the integration is not to be carried out over i. 

Let us now switch gears to model selection. In hypothesis testing, one wishes to determine 
whether the data is more consistent with hypothesis A (eg. that a General Relativity waveform 
correctly models the signal) or with hypothesis B (eg. that a non-General Relativity waveform 
correctly models the signal). Using Bayes' theorem, the PDF for model A given the data is 

As before, p{A) is the prior probability of hypothesis A, namely the strength of our prior belief 
that hypothesis A is correct. The normalization constant p{d) is given by 

p{d) = f dM p{d\M) p(M) . (79) 

where the integral is to be taken over all models. Thus, it is clear that this normalization constant 
does not depend on the model. Similar relations hold for hypothesis B by replacing A ^ B in 



Eq. (78) 



When hypothesis A and B refer to fundamental theories of nature we can take different view- 
points regarding the priors. If we argue that we know nothing about whether hypothesis A or B 
better describes nature, then we would assign equal priors to both hypotheses. If, on the other 
hand, we believe General Relativity is the correct theory of Nature, based on all previous experi- 
ments performed in the Solar System and with binary pulsars, then we would assign p{A) > p{B). 
This assigning of priors necessarily biases the inferences derived from the calculated posteriors, 
which is sometimes heavily debated when comparing Bayesian theory to a frequentist approach. 
However, this "biasing" is really unavoidable and merely a reflection of our state of knowledge of 
nature (for a more detailed discussion on such issues, please refer to [264J). 

The integral over all models in Eq. ( [79| can never be calculated in practice, simply because we 
do not know all models. Thus, one is forced to investigate relative probabilities between models, 
such that the normalization constant p{d) cancels out. The so-called odds-ratio is defined by 

p{A\d) p{A) 

where 'Ba.b = p{d\A) / p{d\B) is the Bayes Factor and the prefactor p{A) /p{B) is the prior odds. 
The odds-ratio is a convenient quantity to calculate because the evidence p{d), which is difficult 
to compute, actually cancels out. Recently, Vallisneri |376| has investigated the possibility of 
calculating the odds-ratio using only frequentist tools and without having to compute full evidences. 
The odds-ratio should be interpreted as the betting-odds of model A over model B. For example, 
an odds-ratio of unity means that both models are equally supported by the data, while an odds- 
ratio of 10^ means that there is a 100 to 1 odds that model A better describes the data than model 
B. 

The main difficulty in Bayesian inference (both in parameter estimation and model selection) is 
sampling the PDF sufficiently accurately. Several methods have been developed for this purpose, 
but currently the two main workhorses in gravitational wave data analysis are Markov chain Monte 
Carlo and Nested Sampling. In the former, one samples the likelihood through the Metropolis- 
Hastings algorithm [282. .193^ I116L 1331] . This is computationally expensive in high-dimensional 
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cases, and thus, there are several techniques to improve the efhciency of the method, e.g. parallel 
tempering |363j . Once the PDF has been sampled, one can then calculate the evidence integral, 
for example via thermodynamic integration |380[ 11521 1379] . In Nested Sampling, the evidence is 
calculated directly by laying out a fixed number of points in the prior volume, which are then 
allowed to move and coalesce toward regions of high posterior probability. With the evidence in 
hand, one can then infer the PDF. As in the previous case. Nested Sampling can be computationally 
expensive in high-dimensional cases. 

Del Pozzo et al |127j were the first to carry out a Bayesian implementation of model selection 
in the context of tests of General Relativity. Their analysis focused on tests of a particular massive 
graviton theory, using the gravitational wave signal from quasi-circular inspiral of non-spinning 
black holes. Cornish, et al |115[ I340J extended this analysis by considering mo del- independent 



deviations from General Relativity, using the ppE approach (Sec. 5.3.41 [4^. Recently, this was 
continued by Li et al |260[ I261| , who carried out a similar analysis to that of Cornish et al [1151 1340] 
on a large statistical sample of advanced LIGO detections using simulated data and a restricted 
ppE model. All of these studies suggest that Bayesian tests of General Relativity are possible, 
given sufficiently high signal-to-noise ratio events. Of course, whether deviations from General 
Relativity are observable will depend on the strong-field character and strength of the deviation, 
as well as the availability of sufficiently accurate General Relativity waveforms. 

4.1.3 Systematics in Model Selection 

The model selection techniques described above are affected by other systematics present in data 
analysis. In general, we can classify these into the following |377) : 

• Mismodeling Systematic, caused by inaccurate models of the gravitational wave template. 

• Instrumental Systematic, caused by inaccurate models of the gravitational wave response. 

• Astrophysical Systematic, caused by inaccurate models of the astrophysical environment. 

Mismodeling systematics are introduced due to the lack of an exact solution to the Einstein equa- 
tions from which to extract an exact template, given a particular astrophysical scenario. Inspiral 
templates, for example, are approximated through post-Newtonian theory and become increasingly 
less accurate as the binary components approach each other. Cutler and Vallisneri |120| were the 
first carry out a formal and practical investigation of such a systematic. 

Mismodeling systematics will prevent us from testing General Relativity effectively with signals 
that we do not understand sufficiently well. For example, when considering signals from black hole 
coalescences, if the the total mass of the binary is sufficiently high, the binary will merge in band. 
The higher the total mass, the fewer the inspiral cycles that will be in band, until eventually only 
the merger is in band. Since the merger phase is the least understood phase, it stands to reason 
that our ability to test General Relativity will deteriorate as the total mass increases. On the other 
hand, for neutron star binaries or very low-mass black hole binaries, the merger phase is expected 
to be essentially out of band for Ad. LIGO (above IkHz), and thus, the noise spectrum itself may 
shield us from our ignorance. 

Instrumental systematics are introduced by our ignorance of the transfer function, which con- 
nects the detector output to the incoming gravitational waves. Through sophisticated calibration 
studies with real data, one can approximate the transfer function very well [3j[l]. However, this 
function is not time- independent, because the noise in the instrument is not stationary or Gaussian. 
Thus, un-modeled drifts in the transfer function can introduce systematics in parameter estimation 
that are as large as 10% in the amplitude and the phase [3]. 

Instrumental systematics can affect tests of General Relativity, if these are performed with 
a single instrument. However, one expects multiple detectors to be online in the future and 
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for gravitational wave detections to be made in several of them simultaneously. Instrumental 
systematics should be present in all such detections, but since the noise will be mostly uncorrelated 
between different instruments, one should be able to ameliorate its effects through cross-correlating 
outputs from several instruments. 

Astrophysical systematics are induced by our lack of a priori knowledge of the gravitational 
wave source. As explained above, matched filtering requires knowledge of a waveform template 
with which to filter the data. Usually, we assume the sources are in a perfect vacuum and isolated. 
For example, when considering inspiral signals, we ignore any third bodies, electric or magnetic 
fields, neutron star hydrodynamics, the expansion of the Universe, etc. Fortunately, however, most 
of these effects are expected to be small: the probability of finding third bodies sufficiently close to 
a binary system is very small [421 ; for low redshift events, the expansion of the Universe induces 
an acceleration of the center of mass, which is also very small 1429] : electromagnetic fields and 
neutron star hydrodynamic effects may affect the inspiral of black holes and neutron stars, but 
not until the very last stages, when most signals will be out of band anyways. For example, tidal 
deformation effects enter a neutron star binary inspiral waveform at 5 post-Newtonian order, which 
therefore affects the signal outside of the most sensitive part of the Adv. LIGO sensitivity bucket. 

Perhaps the most dangerous source of astrophysical systematics is due to the assumptions made 
about the astrophysical systems we expect to observe. For example, when considering neutron 
star binary inspirals, one usually assumes the orbit will have circularized by the time it enters 
the sensitivity band. Moreover, one assumes that any residual spin angular momentum that the 
neutron stars may possess is very small and aligned or counter-aligned with the orbital angular 
momentum. These assumptions certainly simplify the construction of waveform templates, but if 
they happen to be wrong, they would introduce mismodeling systematics that could also affect 
parameter estimation and tests of General Relativity. 

4.2 Burst Analyses 

In alternative theories of gravity, gravitational- wave sources such as core collapse supernovae may 
result in the production of gravitational waves in more than just the plus- and cross-polarizations [3471 
13431 11911 13001 12991 I333J . Indeed, the near-spherical geometry of the collapse can be a source of 
scalar breathing-mode gravitational waves. The precise form of the waveform, however, is unknown 
because it is sensitive to the initial conditions. 

When searching for un-modeled bursts in alternative theories of gravity, a general approach in- 
volves optimized linear combinations of data streams from all available detectors to form maximum 
likelihood estimates of the waveforms in the various polarizations, and the use of null streams. In 
the context of ground-based detectors and General Relativity, these ideas were first explored by 
Giirsel and Tinto |188| and later by Chatterji et al. [M] with the aim of separating false-alarm 
events from real detections. The main idea was to construct a linear combination of data streams 
received by a network of detectors, so that the combination contained only noise. In General Rela- 
tivity, of course, one need only include /i+ and h^ polarizations, and thus a network of 3 detectors 
sufRces. This concept can be extended to develop null tests of General Relativity, as originally 
proposed by Chatziioannou, et al [HS] and recently implemented by Hayama, et al |199j . 

Let us consider a network of D > 6 detectors with uncorrelated noise and a detection by all 
D detectors. For a source that emits gravitational waves in the direction f2, a single data point 
(either in the time-domain, or a time-frequency pixel) from an array of D detectors (either pulsars 
or interferometers) can be written as 

d = Fh + n. (81) 
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Here 
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where n is a vector with the noise. The antenna pattern functions are given by the matrix, 
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(83) 



For simphcity we have suppressed the sky-location dependence of the antenna pattern functions. 



These can be either the interferometric antenna pattern functions in Eqs. (55), or the pulsar 



response functions in Eqs. (65 1. For interferometers, since the breathing and longitudinal antenna 



pattern response functions are degenerate, and even though i^ is a 6 x 13 matrix, there are only 5 
linearly independent vectors [731 [73 HH I199| . 

If we do not know the form of the signal present in our data, we can obtain maximum likelihood 
estimators for it. For simplicity, let us assume the data are Gaussian and of unit variance (the latter 
can be achieved by whitening the data). Just as we did in Eq. (71 1, we can write the probability 
of obtaining datum d, in the presence of a gravitational wave h as 



P{d\h) 



(27r)^/2 



exp 



\d~Fh\ 



(84) 



The logarithm of the likelihood ratio, ie. the logarithm of the ratio of the likelihood when a signal 
is present to that when a signal is absent, can then be written as 



L = ln 



P{d\h) 

Pirn 
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If we treat the waveform values for each datum as free parameters, we can maximize the likelihood 

ratio 

dL 







dh 



h — h-MAX 

and obtain maximum likelihood estimators for the gravitational wave, 

hMAx^iF^Py'F^d. 



(86) 



(87) 



We can further substitute this solution into the likelihood, to obtain the value of the likelihood at 
the maximum, 

^SL = 2L{hMAx) - d^P'^^d, (88) 

where 

The maximized likelihood can be thought of as the power in the signal, and can be used as a 
detection statistic. P^^ is a projection operator that projects the data into the subspace spanned 
by F. An orthogonal projector can also be constructed. 
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which projects the data onto a sub-space orthogonal to F. Thus one can construct a certain linear 
combination of data streams that has no component of a certain polarization by projecting them to 
a direction orthogonal to the direction defined by the beam pattern functions of this polarization 
mode 

This is called a null stream and, in the context of General Relativity, it was introduced as a means 
of separating false-alarm events due, say, to instrumental glitches from real detections [1881194) . 

With just 3 independent detectors, we can choose to eliminate the two tensor modes (the plus- 
and cross-polarizations) and construct a GR null stream: a linear combination of data streams 
that contains no signal consistent within General Relativity, but could contain a signal in another 
gravitational theory, as illustrated in Fig. [4] With such a GR null stream, one can carry out null 
tests of General Relativity and study whether such a stream contains any statistically significant 
deviations from noise. Notice that this approach does not require a template; if one were para- 
metrically constructed, such as in E5J, more powerful null tests could be applied. In the future, 
we expect several gravitational wave detectors to be online: the two Ad. LIGO ones in the United 
States, Ad. VIRGO in Italy, LIGO-India in India, and KAGRA in Japan. Given a gravitational 
wave observation that is detected by all 5 detectors, one can then construct 3 GR null streams, 
each with power in a signal direction. 

For pulsar timing experiments where one is dealing with the data streams of about a few tens 
of pulsars, waveform reconstruction for all polarization states, as well as numerous null streams, 
can be constructed. 




null -direct! on 



Figure 4: Schematic diagram of the projection of the data stream d orthogonal to the GR subspace 
spanned by F+ and F^ , along with a perpendicular subspace, for 3 detectors to build the GR null 
stream. 
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4.3 Stochastic Background Searches 

Much work has been done on the response of ground-based interferometers to non-tensorial polar- 
ization modes, stochastic background detection prospects, and data analysis techniques |268L 12901 
11711 1296| I114J . In the context of pulsar timing, the first work to deal with the detection of such 
backgrounds in the context of alternative theories of gravity is due to Lee et al |254j , who used a 
coherence statistic approach to the detection of non-Einsteinian polarizations. They studied the 
number of pulsars required to detect the various extra polarization modes, and found that pulsar 
timing arrays are especially sensitive to the longitudinal mode. Alves and Tinto |121j also found 
enhanced sensitivity to longitudinal and vector modes. Here we follow the work in [2961 E2] that 
deals with the LIGO and pulsar timing cases using the optimal statistic, a cross-correlation that 
maximizes the signal-to-noise ratio. 

In the context of the optimal statistic, the derivations of the effect of extra polarization states 
for ground-based instruments and pulsar timing are very similar. We begin with the metric per- 
turbation written in terms of a plane wave expansion, as in Eq. ( 49 1 . If we assume that the 
background is unpolarized, isotropic, and stationary, we have that 



Ch\{f, n)hA'{f, n')) = s''{n, n')SAA'S{f - nHAif), 



(92) 



where HA{f) is the gravitational wave power spectrum for polarization A. HA{f) is related to 
the energy density in gravitational waves per logarithmic frequency interval for that polarization 
through 

^Aif)^—^, (93) 
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where Pcrit — 'SHq/Stt is the closure density of the universe, and 
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expansion in Eq. (50), along with Eqs. (92) and (93 1 in Eq. (94), that 



is the energy density in gravitational waves for polarization A. It follows from the plane wave 

(95) 



HA{f) = ^,\f\-'nA(\f\), 



and therefore 
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(96) 



For both ground-based interferometers and pulsar timing experiments, an isotropic stochastic 
background of GWs appears in the data as correlated noise between measurements from different 
instruments. The data set from the a'^ instrument is of the form 



da{t) = Sa{t)+na{t), 



(97) 



where Sa{t) corresponds to the GW signal and na{t) to noise. The noise is assumed in this case to 
be stationary and Gaussian, and uncorrelated between different detectors. 
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for a y^ b. 
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Since the GW signal is correlated, we can use cross-correlations to search for it. The cross- 
correlation statistic is defined as 



i>ab 



T/2 i-T/2 

dt dt'da{t)db{t')Qab{t~t'), 

-T/2 J-T/2 



(100) 



where Qabit — t') is a filter function to be determined. Henceforth, no summation is implied on the 
detector indices {a,b, . . .). At this stage it is not clear why Qabit — t') depends on the pair of data 
sets being correlated. We will show how this comes about later. The optimal filter is determined 
by maximizing the expected signal-to-noise ratio 



SNR = 






(101) 



Here ^ab is the mean (Sab) and aab is the square root of the variance cr^^ = (Sab) ~ (5*06) ^■ 

The expressions for the mean and variance of the cross-correlation statistic, ^ab and cr^j, respec- 
tively, take the same form for both pulsar timing and ground-based instruments. In the frequency 



domain, Eq. ( 100 1 becomes 
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by the convolution theorem, and the mean fi is then 
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where St is the finite time approximation to the delta function, (5t(/) = sin7r/i/(7r/). With this 
in hand, the mean of the cross-correlation statistic is 



OU2 ^oo 

Ma. = Ye^^Ey jf\fr'QabimA{f)rtb{f), 



(104) 



and the variance in the weak signal limit is 

^ah = Wab) ^ (Sab) ~ (5'afc) 

rp />00 2 

df Pai\f\)Pb{\f\) Qabif) 

) 

where the one-sided power spectra of the noise are defined by 

«(/K(/')> = ^'5(/-/Va(|/|), 



(105) 
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in analogy to Eq. (67), where Pa plays here the role of Snif)- 



The mean and variance can be rewritten more compactly if we define a positive-definite inner 
product using the noise power spectra of the two data streams 
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again in analogy to the inner product in Eq. ( |69[ ), when considering inspirals. Using this definition 
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where we recall that the capital Latin indices {A,B, . . .) stand for the polarization content. From 
the definition of the signal-to-noise ratio and the Schwartz's inequality, it follows that the optimal 
filter is given by 



Qabif) = N^ 



(110) 



\f\' Pail f mil f\) ' 

where TV is an arbitrary normalization constant, normally chosen so that the mean of the statistic 
gives the amplitude of the stochastic background. 

The differences in the optimal filter between interferometers and pulsars arise only from dif- 
ferences in the so-called overlap reduction functions, rj^j(/). For ground based instruments, the 
signal data Sa are the strains given by Eq. ( 54 1 . The overlap reduction functions are then given by 



r-,(/) = / dnF^in)F,-in)e 
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where Xa and Xf, are the locations of the two interferometers. The integrals in this case all have 
solutions in terms of spherical Bessel functions 296 , which we do not summarize here for brevity. 



For pulsar timing arrays, the signal data Sa are the redshifts Za, given by Eq. (59 1. The overlap 
reduction functions are then given by 
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where La and Lj are the distances to the two pulsars. For all transverse modes pulsar timing 



experiments are in a regime where the exponential factors in Eq. (112) can be neglected [?7ll92) . 
and the overlap reduction functions effectively become frequency independent. For the + and x 
mode the overlap reduction function becomes 
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(113) 



where ^ab = cos~^(po ■ Pb) is the angle between the two pulsars. This quantity is proportional to 
the Hellings and Downs curve [202] . For the breathing mode, the overlap reduction function takes 
the closed form expression |254j : 



rab=4(3 + COSea6). 



(114) 



For the vector and and longitudinal modes the overlap reduction functions remain frequency de- 
pendent and there are no general analytic solutions. 

The result for the combination of cross-correlation pairs to form an optimal network statistic 
is also the same in both ground-based interferometer and pulsar timing cases: a sum of the cross- 
correlations of all detector pairs weighted by their variances. The detetctor network optimal 
statistic is. 
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where '^ab ^^ ^ sum over all detector pairs. 

In order to perform a search for a given polarization mode one first needs to compute the overlap 



reduction functions (using either Eq. ( 111 ) or ( 112 1) for that mode. With that in hand and a form 



for the stochastic background spectrum Oyi(/), one can construct optimal filters for all pairs in 
the detector network using Eq. (110), and perform the cross-correlations using either Eq. (100) (or 
equivalently Eq. (102)). Finally, we can calculate the overall network statistic Eq. (115), by first 
finding the variances using Eq. (105). 
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It is important to point out that the procedure outhned above is straightforward for ground- 
based interferometers. Pulsar timing data, however, are irregularly sampled, and have a pulsar 
timing model subtracted out. This needs to be accounted for, and generally, a time-domain ap- 
proach is more appropriate for these data sets. The procedure is similar to what we have outlined 
above, but power spectra and gravitational wave spectra in the frequency domain need to be re- 
placed by auto-covariance and cross-covariance matrices in the time domain that account for the 
model fitting (for an example of how to do this see ^ 147j ). 

Interestingly, Nishizawa et al. |296j show that with three spatially separated detectors the ten- 
sor, vector, and scalar contributions to the energy density in GWs can be measured independently. 
Lee et al. [254j and Alves and Tinto |121j showed that pulsar timing experiments are especially 
sensitive to the longitudinal mode, and to a lesser extent the vector modes. Chamberlin and 
Siemens [5^ showed that the sensitivity of the cross-correlation to the longitudinal mode using 
nearby pulsar pairs can be enhanced significantly compared to that of the transverse modes. For 
example, for the NANOGrav pulsar timing array, two pulsar pairs separated by 3° result in an 
enhancement of 4 orders of magnitude in sensitivity to the longitudinal mode relative to the trans- 
verse modes. The main contribution to this effect is due to GWs that are coming from roughly 
the same direction as the pulses from the pulsars. In this case, the induced redshift for any GW 
polarization mode is proportional to /L, the product of the GW frequency and the distance to the 
pulsar, which can be large. When the GWs and the pulse direction are exactly parallel, the redshift 
for the transverse and vector modes vanishes, but it is proportional to fL for the scalar-longitudinal 
mode. 

In [25511254] . Lee et al. studied the detectability of massive gravitons in pulsar timing arrays. 



They introduced a modification to Eq. ( 56 1 to account for graviton dispersion, and found the 



modified overlap reduction functions (i.e. modifications to the Hellings-Downs curves (113l) for 
various values of the graviton mass. They conclude that a large number of stable pulsars (> 60) 
are required to distinguish between the massive and massless cases. 
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5 Compact Binary Tests 

In this chapter, we discuss gravitational wave tests of General Relativity with signals emitted by 
compact binary systems. We begin by explaining the difference between direct and generic tests. 
We then proceed to describe the many direct or top-down tests and generic or bottom-up tests that 
have been proposed once gravitational waves are detected, including tests of the no-hair theorems. 
We concentrate here only on binaries composed of compact objects, such as neutron stars, black 
holes or other compact exotica. We will not discuss tests one could carry out with electromagnetic 
information from binary (or double binary) pulsars, as these are already described in |399| . We 
will also not review tests of General Relativity with accretion disk observations, for which we refer 
the interested reader to 



5.1 Direct and Generic Tests 

Gravitational wave tests of Einstein's theory can be classified in two distinct subgroups: direct tests 
and generic tests. Direct tests employ a top-down approach, where one starts from a particular 
modified gravity theory with a known action, derives the modified field equations and solves them 
for a particular gravitational wave-emitting system. On the other hand, generic tests adopt a 
bottom-up approach, where one takes a particular feature of General Relativity and asks what 
type of signature its absence would leave on the gravitational wave observable; one then asks 
whether the data presents a statistically significant anomaly pointing to that particular signature. 

Direct tests have by far been the traditional approach to testing General Relativity with grav- 
itational waves. The prototypical example here are tests of Jordan-Brans-Dicke-Fierz theory. As 
described in Sec. [2] one can solve the modified field equations for a binary system in the post- 
Newtonian approximation to find a prediction for the gravitational wave observable, as we will 
see in more detail later in this section. Other examples of direct tests include those concerning 
modified quadratic gravity models and non-commutative geometry theories. 

The main advantage of such direct tests is also its main disadvantage: one has to pick a 
particular modified gravity theory. Because of this, one has a well-defined set of field equations 
that one can solve, but at the same time, one can only make predictions about that particular 
modified gravity model. Unfortunately, we currently lack a particular modified gravity theory that 
is particularly compelling; many modified gravity theories exists, but none possess all the criteria 
described in Sec. [2] Lacking a clear alternative to General Relativity, it is not obvious which 
theory one should pick. Given that the full development (from the action to the gravitational 
wave observable) of any particular theory can be incredibly difficult, time and computationally 
consuming, carrying out direct tests of all possible modified gravity models once gravitational 
waves are detected is clearly unfeasible. 

Given this, one is led to generic tests of General Relativity, where one asks how the absence of 
specific features contained in General Relativity could impact the gravitational wave observable. 
For example, one can ask how such observable would be modified if the graviton had a mass, if the 
gravitational interaction were Lorentz or parity violating, or if there existed large extra dimensions. 
From these general considerations, one can then construct a "meta"-observable, ie. one that does 
not belong to a particular theory, but that interpolates over all known possibilities in a well-defined 
way. This model has become to be known as the parameterized post-Einsteinian framework, in 
analogy to the parameterized post-Newtonian scheme used to test General Relativity in the Solar 
System |399| . Given such a construction, one can then ask whether the data points to a statistically 
significant deviation from General Relativity. 

The main advantage of generic tests is precisely that one does not have to specify a particular 
model, but instead one lets the data select whether it contains any statistically significant deviations 
from our canonical beliefs. Such an approach is, of course, not new to physics, having most recently 
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been successfully employed by the WMAP team [55]. The intrinsic disadvantage of this method 
is that, if a deviation is found, there is no one-to-one mapping between it and a particular action, 
but instead one has to point to a class of possible models. Of course, such a disadvantage is not 
that limiting, since it would provide strong hints as to what type of symmetries or properties of 
General Relativity would have to be violated in an ultra-violet completion of Einstein's theory. 

5.2 Direct Tests 

5.2.1 Scalar- Tensor Theories 

Let us first concentrate on Jordan-Brans-Dicke-Fierz theory, where black holes and neutron stars 
have been shown to exist. In this theory, the gravitational mass depends on the value of the 
scalar field, as Newton's constant is effectively promoted to a function, thus leading to violations 
of the weak-equivalence principle f l45l I395[ I402J . The usual prescription for the modeling of binary 
systems in this theory is due to Eardley [145) . He showed that such a scalar-field effect can be 
captured by replacing the constant inertial mass by a function of the scalar field in the distributional 
stress-energy tensor and then Taylor expanding about the cosmological constant value of the scalar 
field at spatial infinity, ie. 
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These sensitivities encode how the gravitational mass changes due to a non-constant scalar field; 
one can think of them as measuring the gravitational binding energy per unit mass. In carrying 
out this expansion, one assumes that the scalar field takes on a constant value at spatial infinity 
(j) — >■ (/>o, disallowing any homogeneous solution to the scalar field evolution equation [Eq. (19)]. 



Moreover, one assumes that the internal gravitational field of each body leads to a non-trivial 
variation of the scalar field, which then leads to modifications to the gravitational binding energies 
of the bodies. 

With this at hand, one can solve the massless Jordan-Brans-Dicke-Fierz modified field equations 
[Eq. (|l9|)] for the non-dynamical, near-zone field of a compact objects to obtain 
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where we have defined the spatial field point distance r^ = |a;' — a;^|, the parameterized post- 
Newtonian quantity 7 = (1 +aJBD)(2 + u;bd)^^ and we have chosen units in which G = c = 1. This 
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solution is obtained in a post-Newtonian expansion i68], where the ehipses represent higher order 
terms in Vajc and ma/ra- From such an analysis, one can also show that compact objects follow 
geodesies of such a spacetime, to leading order in the post-Newtonian approximation |145| . except 
that Newton's constant is replaced by G — )■ S12 = 1 — (si -I- S2 — 2siS2)(2 -I- wbd)~^, in geometric 
units. 

As is clear from the above analysis, black hole and neutron star solutions in this theory generi- 
cally depend on the quantities wbd and s^. The former determines the strength of the correction, 
with the theory reducing to General Relativity in the wbd ~^ 00 limit. The latter depends on the 
compact object that is being studied. For neutron stars, this quantity can be computed as follows. 
First, using the Tolman-Oppenheimer-Volkoff equation, one notes that the mass m oc TV ex G~^/^, 
for a fixed equation of state and central density, with N the total baryon number. Thus, using 



Eq. (1171, one has that 

3 r d{\nma) 
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where the derivative is to be taken holding G fixed. In this way, given an equation of state 
and central density, one can compute the gravitational mass as a function of baryon number, 
and from this, obtain the neutron star sensitivities. Eardley jl45j . Will and Zaglauer j402j . and 
Zaglauer |435| have shown that these sensitivities are always in the range s^ g (0.19, 0.3) for a soft 
equation of state and Sa G (0.1,0.14) for a stiff one, in both cases monotonically increasing with 
mass in m^ e (1.1, 1.5) Af©. Recently, Gralla [179j has found a more general method to compute 
sensitivities is generic modified gravity theories. 

What is the sensitivity of black holes in generic scalar-tensor theories? Will and Zaglauer |435j 
have argued that the no-hair theorems require Sa = 1/2 for all black holes, no matter what their 
mass or spin is. As already explained in Sec. [21 stationary black holes that are the byproduct of 
gravitational collapse (ie. with matter that satisfies the energy conditions) in a general class of 
scalar-tensor theories are identical to their General Relativity counterparts [1961 13681 11441 [559p | 
This is because the scalar field satisfies a free wave equation in vacuum, which forces the scalar field 
to be constant in the exterior of a stationary, asymptotically fiat spacetime, provided one neglects 



a cosmological homogeneous solution. If the scalar field is to be constant, then by Eq. (118 1, 
Sa = 1/2 for a single black hole spacetime. 

Such an argument formally applies only to stationary scenarios, so one might wonder whether 
a similar argument holds for binary systems that are in a quasi-stationary arrangement. Will 
and Zaglauer |435j and Mirshekari and Will |283j extended this discussion to quasi-stationary 
spacetimes describing binary black holes to first and second post-Newtonian order respectively. 
They argued that the only possible deviations from -0 = are due to tidal deformations of the 
horizon due to the companion, which are known to arise at very high order in post-Newtonian 
theory, ip = 0[{ma/ra)^]. Recently, Yunes, et al. [426) extended this argument further by showing 
that to all orders in post-Newtonian theory, but in the extreme mass-ratio limit, black holes cannot 
have scalar hair in generic scalar tensor theories. Finally, Healy, et al. [201] have carried out a 
full numerical simulations of the non-linear field equations confirming this argument in the full 
non-linear regime. 

The activation of dynamics in the scalar field for a vacuum spacetime requires either a non- 
homogeneous distribution of initial scalar field (violating the constant cosmological scalar field 
condition at spatial infinity) or a pure geometrical source to the scalar field evolution equation. 
The latter would lead to the quadratic modified gravity theories discussed in Sec. |2.3.3[ As for 
the former, Horbatsch and Burgess \2U^ have argued that if, for example, one lets -0 = ^t, which 



^One should note in passing that more general black hole solutions in scalar-tensor theories have been found [2351 
1841 . These, however, usually violate the weak energy condition, and sometimes they require unreasonably small 
values of (.i^bd that have already been ruled out by observation. 
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clearly satisfies Dif; = in a Minkowski background Q then a Schwarzschild black hole will acquire 
modifications that are proportional to fi. Such cosmological hair, however, is likely to be suppressed 
by a long time-scale; in the example above fj, must have units of inverse time, and if it is to be 
associated with the expansion of the Universe, then it would be natural to assume /i = 0{H), where 
H is the Hubble parameter. Therefore, although such cosmological hair might have an effect on 
black holes in the early Universe, it should not affect black hole observations at moderate to low 
redshifts. 

Scalar field dynamics can be activated in non- vacuum spacetimes, even if initially the stars are 
not scalarized provided one considers a more general scalar tensor theory, like the one introduced 



by Damour and Esposito-Farese |122| 1123] , As discussed in Sec. 2.3.1 when the conformal factor 
takes on a particular functional form, non-linear effects induced when the gravitational energy 
exceeds a certain threshold can spontaneously scalarize merging neutron stars, as demonstrated 
recently by Barausse, et al |46j . Therefore, neutron stars in binaries are likely to have hair in 
generic scalar-tensor theories, even if they start their inspiral unscalarized. 

What do gravitational waves look like in Jordan-Brans-Dicke-Fierz theory? As described in 
Sec. |2.3.H both the scalar field perturbation ip and the new metric perturbation 9'^'^ satisfy a 
sourced wave equation [Eq. (19l], whose leading-order solution for a two-body inspiral is |397j 
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where R is the distance to the detector, n' is a unit vector pointing toward the detector, r is the 



magnitude of relative position vector x* 



- ^2, with x^ the trajectory of body a, fi — mim2/m 



is the reduced mass and m — rrii +m2 is the total mass, WJ2 
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We have here also introduced multi-index notation, such that A*-'-- = A^A^ . . .. Such a solution is 
derived using the Lorenz gauge condition 9^'^ ^i, = and in a post-Newtonian expansion, where we 
have left out subleading terms of relative order VI2 or m/r. 

Given the new metric perturbation 9^^ , one can reconstruct the gravitational wave h^^ metric 
perturbation, and from this, the response function, associated with the quasi-circular inspiral of 
compact binaries. After using Kepler's third law to simplify expressions [to = (312"^/'''^)^, where 
oj is the orbital angular frequency and m is the total mass and r is the orbital separation], one 
finds for a ground-based L-shaped detector [^ : 



hit) 



R 



(2^McF) 



-FFb sin^ 



2/3g-2i$ 



[^+(1 + 



;M, 



COS^ i] 



^1/3, 



2iFx COS t] 



■7 



,7i/5:_^(2,rMei^)^'-e 



i$ 



5(l-7)i^bSini 



^ {2nM^Ff' 1^ 



-Fb (F + 2A) , 



-S' 



(127) 



where 77 = ix/m is the symmetric mass ratio, Mc = rf'^m is the chirp mass, t is the inclination 
angle, and where we have used the beam-pattern functions in Eq. ( 55 ). In Eq. ( 127 ) and henceforth. 



'^The scalar field of Horbatsch and Burgess satisfies Di/) : 
gu = 0. 



: /^rj^j, and thus Di/" = for stationary spacetimes with 
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we linearize all expressions in 1 — 7 ^ 1. Jordan-Brans-Dicke-Fierz theory predicts the generic 
excitation of three polarizations: the usual plus and cross polarizations, and a breathing, scalar 
mode. We see that the latter contributes to the response at two, one and zero times the orbital 
frequency. One should note that all of these corrections arise during the generation of gravitational 
waves, and not due to a propagation effect. In fact, gravitational waves travel at the speed of light 
(and the graviton remains massless) in standard Jordan-Brans-Dicke-Fierz theory. 

The quantities $ and F are the orbital phase and frequency respectively, which are to be found 
by solving the differential equation 



^ = (1-7)^^^ 
at IT 



-M^2(27rMci^)3 
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where the ellipses stand for higher-order terms in the post-Newtonian approximation. The first 
term here corresponds to dipole radiation, which is activated by the scalar mode. That is, the scalar 
field carries energy away from the system modifying the energy balance law to [3971 13421 1401) 
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ebd = -gSL^ V (^ -^'^^ y ^'^^ (7) 
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where the ellipses stand again for higher-order terms in the post-Newtonian approximation. Solving 
the frequency evolution equation perturbative in I/cjbd ^ 1, one finds 
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where we have defined u = {2nM.cFy^^. In deriving these equations, we have neglected the last 



term in Eq. (128), as this is a constant that can be reabsorbed into the chirp mass. Notice that 



since the two definitions of chirp mass differ only by a term of O(ajgp), the first term of Eq. (128) 
is not modified. 

One of the main ingredients that goes into parameter estimation is the Fourier transform of 
the response function. This can be estimated in the so-called stationary-phase approximation, for 
a simple, non-spinning, quasi-circular inspiral. In this approximation, one assumes the phase is 
changing much more rapidly than the amplitude [SH 11181 11381 I418J . One finds [SSI 
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where the Brans-Dicke correction is kept only to leading order in Wg^ and v. In writing the Fourier 
response in this way, we had to redefine the phase of coalescence via 



(j)W = $^ — J^ 2 < arctan 



2 cos L Fx 
(1 + cos2 i) F+ 



(1 



>s'^i)FxFy 



x-fh 



^BD (1 + cos2 lY Fl + 4cos2 lFI 



(136) 



where dg^rn is the Kronecker delta and <i>c is the GR phase of coalescence (defined as an integration 
constant when the frequency diverges) . Of course, in this calculation we have neglected amplitude 
corrections that arise purely in General Relativity, if one were to carry out the post-Newtonian 
approximation to higher order. 

Many studies have been carried out to determine the level at which such corrections to the 
waveform could be measured or constrained once a gravitational wave has been detected. The first 
such study was carried out by Will |397) . who determined that given a LIGO detection at signal- 
to-noise ratio p = 10 of a (1.4,3)M0 black- hole/neutron-star non-spinning, quasi-circular inspiral, 
one could constrain wbd > 10^. Scharre and Will |342j carried out a similar analysis but for a LISA 
detection with p = 10 of a (1.4, 1O'^)M0 intermediate-mass black-hole/neutron-star, non-spinning, 
quasi-circular inspiral, and found that one could constrain wbd > 2.1 x lO'*. Such an analysis was 
then repeated by Will and Yunes |401j but as a function of the classic LISA instrument. They 
found that the bound is independent of the LISA arm length, but inversely proportional to the 
LISA position noise error, if the position error noise dominates over laser shot noise. All such 
studies considered an angle-averaged signal that neglected the spin of either body, assumptions 
that were relaxed by Berti, et al. [SSIIS^. They carried out Monte-Carlo simulations over all sky 
position of signals that included spin-orbit precession to find that the projected bound with LISA 
deteriorates to wbd > 0.7 x 10^ for the same system and signal-to-noise ratio. This was confirmed 
and extended by Yagi, et al. |411) , who in addition to spin-orbit precession allowed for non-circular 
(eccentric) inspirals. In fact, when eccentricity is included, the bound deteriorates even further 
to wbd > 0.5 X 10'*. The same authors also found that similar gravitational wave observations 
with the next-generation detector DECIGO could constrain wbd > 1.6 x 10^. Similarly, for a non- 
spinning neutron star/black hole binary, the future ground-based detector, the Einstein Telescope 
(ET) |325j . could place constraints about 5 times stronger than the Cassini bound, as shown in |33) . 

All such projected constraints are to be compared with the current Solar System bound of 
i^BD > 4 X lO'* placed through the tracking of the Cassini spacecraft |66|- Table [l] presents all such 
bounds for ease of comparisorj^ normalized to signal-to-noise ratio of 10. As it should be clear, it 
is unlikely that LIGO observations will be able to constrain ojbd better than current Solar System 
bounds. In fact, even LISA would probably not be able to do better than the Cassini bound. 
Table [l] also shows that the inclusion of more complexity in the waveform seems to dilute the level 
at which wbd can be constrained. This is because the inclusion of eccentricity and spin forces one 
to introduce more parameters in the waveform, without these modifications truly adding enough 
waveform complexity to break the induced degeneracies. One would then expect that the inclusion 
of amplitude modulation due to precession and higher harmonics should break such degeneracies, 
at least partially, as was found for massive binary black holes |249l I250J . Even then, however, 
it seems reasonable to expect that only third-generation detectors will be able to constrain wbd 
beyond Solar System levels. 

The main reason that Solar System constraints of Jordan-Brans-Dicke-Fierz theory cannot 
be beaten with gravitational wave observations is that the former are particularly well-suited to 
constrain weak-field deviations of General Relativity. One might have thought that scalar-tensor 
theories constitute strong-field tests of Einstein's theory, but this is not quite true, as argued in 
Sec. |2.3.1[ One can see this clearly by noting that scalar-tensor theory predicts dipolar radiation. 



*A11 LISA bounds refer to the classic LISA configuration. 
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Table 1: 


Comparison of 


Proposed Tests of Scalar- Tensor Theories. 


Reference 


Binary Mass 


^bd[10^] 


Properties 


m 


X 


4 


Solar System 


[397] 


(1.4,3)M0 


0.1 


LIGO, Fisher, Ang. Ave. 
circular, non-spinning 


|M51 


(1.4, lO^)A-f0 


24 


LISA, Fisher, Ang. Ave. 
circular, non-spinning 


m\ 


(1.4,10^)Mo 


20 


LISA study, Fisher, 
Ang. Ave., circular, non-spinning 


m 


(1.4, 1O3)M0 


0.7 


LISA, Fisher, Monte-Carlo, 
circular, w/spin-orbit 


mn 


(1.4, 1O^)M0 


0.5 


LISA, Fisher, Monte-Carlo 
eccentric, spin-orbit 


mn 


(1.4,1O)M0 


160 


DECIGO, Fisher, Monte-Carlo 
eccentric, spin-orbit 


[33] 


(1.4,1O)M0 


10 


FT, Fisher, Ang. Ave. 
circular, non-spinning 



which dominates at low velocities over the General Relativity prediction (precisely the opposite 
behavior that one would expect from a strong- field modification to Einstein's theory). 

One should note, however, that all the above analysis considered only the inspiral phase of 
coalescence, usually truncating their study at the innermost stable-circular orbit. The merger 
and ringdown phases, where most of the gravitational wave power resides, have so far been mostly 
neglected. One might expect that an increase in power will be accompanied by an increase in signal- 
to-noise ratio, thus allowing us to constrain wbd further, as this scales with 1/SNR [233] . Moreover, 
during merger and ringdown, dynamical strong-field gravity effects in scalar-tensor theories could 
affect neutron star parameters and their oscillations 356^ , as well as possibly inducing spontaneous 
scalarization [46] . All of these non-linear effect could easily lead to a strengthening of projected 
bounds. To date, however, no detailed analysis has attempted to determine how well one could 
constrain scalar-tensor theories using full information about the entire coalescence of a compact 
binary. 

The subclass of scalar-tensor models described by Jordan-Brans-Dicke-Fierz theory is not the 
only type of model that can be constrained with gravitational wave observations. In the extreme 
mass-ratio limit, for binaries consisting of a stellar-mass compact object spiraling into a supermas- 
sive black hole, Yunes, et al |426j have recently showed that generic scalar-tensor theories reduce 
to either massless or massive Jordan-Brans-Dicke-Fierz theory. The inclusion of a scalar field mass 
leads to an interesting possibility: floating orbits |87| . Such orbits arise when the small compact 
object experiences superradiance, leading to resonances in the scalar flux that can momentarily 
counteract the gravitational wave flux, leading to a temporarily stalled orbit that greatly modifies 
the orbital phase evolution. These authors showed that if an extreme mass-ratio inspiral is de- 
tected with a template consistent with General Relativity, this alone allows us to rule out a large 
region of (msjiOBi)) phase space, where nis is the mass of the scalar (see Fig. 1 in |426| ). This 
is because if such an inspiral had gone through a resonance, a General Relativity template would 
be grossly different from the signal. Such bounds are dramatically stronger than the current most 
stringent bound ojbd > 4 x 10^ and ms < 2.5 x lO"^*' eV obtained from Cassini measurements of 
the Shapiro time-delay in the Solar System [THj. Even if resonances are not hit, Berti et al [SS] have 
estimated that second-generation ground-based interferometers could constrain the combination 
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'^s/('^bd)^ ^ lO^^'^ eV with the observation of gravitational waves from neutron-star/binary 
inspirals at signal-to-noise ratio of 10. These bounds can also be stronger than current constraints, 
especially for large scalar mass. 

Lastly one should mention possible gravitational wave constraints on other types of scalar tensor 
theories. Let us first consider Brans-Dicke type scalar-tensor theories, where the coupling constant 
is allowed to vary. Will [397] has argued that the constraints described in Table fl| go through, with 
the change 



2Sl,2 291.2 

■ > — 



^BD ^ + ^BD 



2^3 j3 



(3 + 2cjbd)2 



(137) 



where w^q = dujBB/dtj). In the wbd ^ 1 limit, this implies the replacement wbd ^ ^^^BDil + 
'^bd/(2'^bd)]"^- *-'f course, this assumes that there is neither a potential nor a geometric source 
driving the evolution of the scalar field. 

Another interesting scalar-tensor theory to consider is that studied by Damour and Esposito- 



Farese [1221 I123J . As explained in Sec. 2.3.1 this theory is defined by the action of Eq. (14) 
with the conformal factor A{tp) = e^"^' . In standard Brans-Dicke theory, binary neutron stars 
generate gravitational waves that are barely modified from the General Relativity expectation. 
This is because the Brans-Dicke correction is proportional to the difference in the neutron star 



sensitivities (see Eqs. ( 125 ) and ( 135 1) , which must be close to zero for neutron stars of similar mass. 
Similarly, binary black holes generate gravitational waves that are identical to those of General 
Relativity, since the difference in black hole sensitivities is exactly zero. In the theory considered 
by Damour and Esposito-Farese, however, when the gravitational energy is large enough, as in the 
very late inspiral, non-linear effects can lead to drastic modifications from the General Relativity 
expectation, such as spontaneously scalarization |46j . Unfortunately, most of this happens at 
rather high frequency, and thus, it is not clear whether such effects are observable with current 
ground-based detectors. 

5.2.2 Modified Quadratic Gravity 

Black holes exist in the classes of modified quadratic gravity that have been so far considered. In 



non-dynamical theories (when /3 = and the scalar-fields are constant, refer to Eq. (25)), Stein 
and Yunes [434j have shown that all metrics that are Ricci tensor flat are also solutions to the 
modified field equations (see also [324J). This is not so for dynamical theories, since then the d field 
is sourced by curvature, leading to corrections to the field equations proportional to the Riemann 
tensor and its dual. 

In dynamical Chern-Simons gravity, stationary and spherically symmetric spacetimes are still 
described by General Relativity solutions, but stationary and axisymmetric spacetimes are not. 
Instead, they are represented by |427[ I241J 



^r4 1^ "'' y v^ Toy 



dsls = dsh,, + 7^^1 + V— + T7^— ) ^^^' ^^^'^^ + 0(aVM2) , (138) 



with the scalar field 



where dsxorr i^ ^^^ ^i^^ element of the Kerr metric and we recall that acs = — 4a4 in the notation 
of Sec. |2.3.3[ These expressions are obtained in Boyer-Lindquist coordinates and in the small- 
rotation/small-couphng limit to 0{a/M) in [i?fl[MT] and to 0{a'^/M^) in gTB]. The linear- in- 
spin corrections modify the frame-dragging effect and they are of 3.5 post-Newtonian order. The 
quadratic-in-spin corrections modify the quadrupole moment, which induces 2 post-Newtonian 
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order corrections to the binding energy. The stabihty of these black holes, however, has not yet 
been demonstrated. 

In Einstein-Dilaton-Gauss-Bonnet gravity, stationary and spherically symmetric spacetimes are 
described, in the small-coupling approximation, by the line element |434j 

dsluGB - -/schw (1 + h) df + /g-J^^ (1 + k) dr^ + r^dn^ , (140) 

in Schwarzschild coordinates, where dO? is the line element on the two-sphere, /schw = 1 ^ 2M/r 
is the Schwarzschild factor and we have defined 

/^=^r77j:77— (1 + 26- + -— + -— -80— ) , (141) 
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while the corresponding scalar field is 
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This solution is not restricted just to Einstein-Dilaton-Gauss-Bonnet gravity, but it is also the most 
general, stationary and spherically symmetric solution in quadratic gravity. Linear slow-rotation 
corrections to this solution have been found in |310j . Although the stability of these black holes has 
not yet been demonstrated, other dilatonic black hole solutions obtained numerically (equivalent 
to those in Einstein-Dilaton-Gauss-Bonnet theory in the limit of small fields) [228] have been found 
to be stable under axial perturbations |2291 13701 1308| . 

Neutron stars also exist in quadratic modified gravity. In dynamical Chern-Simons gravity, 
the mass-radius relation remains unmodified to first-order in the slow-rotation expansion, but the 
moment of inertia changes to this order |430[ [T7] , while the quadrupole moment and the mass 
measured at spatial infinity change to quadratic order in spin [409] . This is because the mass- 
radius relation, to first order in slow-rotation, depends on the spherically-symmetric part of the 
metric, which is unmodified in dynamical Chern-Simons gravity. In Einstein-Dilaton-Gauss-Bonnet 
gravity, the mass-radius relation is modified [307]. As in General Relativity, these functions must 
be solved for numerically and they depend on the equation of state. 

Gravitational waves are also modified in quadratic modified gravity. In dynamical Chern- 
Simons gravity, Garfinklc et al [170 J have shown that the propagation of such waves on a Minkowski 
background remains unaltered, and thus, all modifications arise during the generation stage. In 
Einstein-Dilaton-Gauss-Bonnet theory, no such analysis of the propagation of gravitational waves 
has yet been carried out. Yagi et al |408j studied the generation mechanism in both theories 
during the quasi-circular inspiral of comparable-mass, spinning black holes in the post-Newtonian 
and small-coupling approximations. They found that a standard post-Newtonian analysis fails for 
such theories because the assumption that black holes can be described by a distributional stress- 
energy tensor without any further structure fails. They also found that since black holes acquire 
scalar hair in these theories, and this scalar field is anchored to the curvature profiles, as black 
holes move, the scalar fields must follow the singularities, leading to dipole scalar-field emission. 

During a quasi-circular inspiral of spinning black holes in dynamical Chern-Simons gravity, 
the total gravitational wave energy flux carried out to spatial infinity (equal to minus the rate of 
change of a binary's binding energy by the balance law) is modified from the General Relativity 
expectation to leading order by |408j 
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S, 



GR 



^^^" lie [^ + ' ^(^ ■ '-)>] + S^ (^^^^ ('*-'- - '"-"-) > r'^ • (^'') 
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due to scalar field radiation and corrections to the metric perturbation that are of magnetic- 
type, quadrupole form. In this equation, E'en = (32/5)?7^i;^2 is the leading order General Rel- 
ativity prediction for the total energy flux, (,4 = a\/{j3Krn^) is the dimensionless Chern-Simons 
coupling parameter, W12 is the magnitude of the relative velocity with unit vector v\2, A* = 
{m2/mi){ai/m)S\ — {mi/m2){a2/'rn)S\, where a^ is the Kerr spin parameter of the Ath black 
hole and S\ is the unit vector in the direction of the spin angular momentum, the unit vector n\2 
points from body one to two, and the angle brackets stand for an average over several gravitational 
wave wavelengths. If the black holes are not spinning, then the correction to the scalar energy flux 
is greatly suppressed |408) 

ApCS r, 

"-^no-spin ■^ r2 /- 14 /I/ICA 

— s = o VC4W12 , (145) 

where we have defined the reduced mass difference 5m = (jtii — ni2)/ni. Notice that this is a 7 



post-Newtonian order correction, instead of a 2 post-Newtonian correction as in Eq. (144). In 
the non-spinning limit, the dynamical Chern-Simons correction to the metric tensor induces a 6 
post-Newtonian order correction to the gravitational energy flux |408j , which is consistent with the 
numerical results of |309) . 

On the other hand, in Einstein-Dilaton-Gauss-Bonnet gravity, the corrections to the energy 
flux are [408] 

rpEDGB r 

no— spin — 4 r2 /■ —2 lTAe\ 

^^^ = 96^ '5„C3.i2, (146) 

which is a —1 post-Newtonian correction. This is because the scalar field i?edgb behaves like a 



monopole (see Eq. (143)), and when such a scalar monopole is dragged by the black hole, it emits 
electric-type, dipole scalar radiation. Any hairy black hole with monopole hair will thus emit 
dipolar radiation, leading to —1 post-Newtonian corrections in the energy flux carried to spatial 
infinity. 

Such modifications to the energy fiux modify the rate of change of the binary's binding energy 
through the balance law, E = —Ei,, which in turn modify the rate of change of the gravitational 
wave frequency and phase, F ~ —E (dEi^/dF)^^ . For dynamical Chern-Simons gravity (when 
the spins are aligned with the orbital angular momentum) and for Einstein-Dilaton-Gauss-Bonnet 
theory (in the non-spinning case) , the Fourier transform of the gravitational wave response function 
in the stationary phase approximation becomes [4081 I415J 

/ldCS,EDGB = /iGRe''3^c=.™<^B"''^''"°''" , (147) 

where hcR is the Fourier transform of the response in General Relativity, u = (irMcfY^^ with / 
the gravitational wave frequency and [4081 I415J 
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where we have defined the symmetric and antisymmetric spin combinations Xs,a = (cii/mi ± 02/7722) /2. 
We have here neglected any possible amplitude correction, but we have included both deforma- 
tions to the binding energy and Kepler's third law, in addition to changes in the energy flux, when 
computing the phase correction. In Einstein-Dilaton-Gauss-Bonnet theory, however, the binding 
energy is modifled at higher post-Newtonian order, and thus, corrections to the energy flux control 
the modifications to the gravitational wave response function. 
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From the above analysis, it should be clear that the corrections to the gravitational wave 
observable in quadratic modified gravity are always proportional to the quantity ^3^4 = ^3.4/™* = 
a§ 4/(/3Krn'*). Thus, any measurement that is consistent with General Relativity will allow a 
constraint of the form ^3^4 < NS, where TV is a number of order unity, and 6 is the accuracy of the 
measurement. Solving for the coupling constants of the theory, such a measurement would lead 
to ^34 < (NSy^'^m [353 ' . Therefore, constraints on quadratic modified gravity will weaken for 
systems with larger characteristic mass. This can be understood by noticing that the corrections 
to the action scale with positive powers of the Riemann tensor, while this scales inversely with the 
mass of the object, ie. the smaller a compact object is, the larger its curvature. Such an analysis 
then automatically predicts that LIGO will be able to place stronger constraints than LISA-like 
missions on such theories, because LIGO operate in the 100 Hz frequency band, allowing for the 
detection of stellar-mass inspirals, while LISA-like missions operate in the mHz band, and is limited 
to supermassive black holes inspirals. 

How well can these modifications be measured with gravitational wave observations? Yagi, et 
al |408j predicted, based on the results of Cornish et al [115j . that a sky-averaged LIGO gravita- 
tional wave observation with signal-to-noise ratio of 10 of the quasi-circular inspiral of non spinning 

1 /4 

black holes with masses (6, 12) Mq would allow a constraint of ^3 < 20 km, where we recall that 
£.3 — Oii/iPi^)- A similar sky-averaged, eLISA observation of a quasi-circular, spin-aligned black 
hole inspiral with masses (10^,3 x 10^) Mq would constrain 1^3 < 10^ km '408' . The loss in 
constraining power comes from the fact that the constraint on ^3 will scale with the total mass of 
the binary, which is six orders of magnitude larger for space-borne sources. This constraints are not 
stronger than current bounds from the existence of compact objects |307j (^3 < 26 km) and from 
the change in the orbital period of the low-mass x-ray binary A0620 — 00 (^3 < 1.9 km) |406j . but 
they are independent of the nature of the object and sample the theory in a different energy scale. 
In dynamical Chern-Simons gravity, one expects similar projected gravitational wave constraints 
on ^4, namely ^4' < 0(M), where M is the total mass of the binary system in kilometers. There- 
fore, for binaries detectable with ground-based interferometers, one expects constraints of order 
^4 < 10 km. In this case, such a constraint would be roughly six orders of magnitude stronger 
than current LAGEOS bounds |17j . Dynamical Chern-Simons gravity cannot be constraint with 
binary pulsar observations, since the theory's corrections to the post-Keplerian observables are too 
high post-Newtonian order, given the current observational uncertainties j409| . The gravitational 
wave constraint, however, is more difficult to achieve in the dynamical Chern-Simons case, because 
the correction to the gravitational wave phase is degenerate with spin. Yagi, et al. plS , however, 
argued that precession should break this degeneracy, and if a signal with sufficiently high signal- 
to-noise ratio is observed, such bounds would be possible. One must be careful, of course, to check 
that the small-coupling approximation is still satisfied when saturating such a constraint |415) . 

5.2.3 Non-Commutative Geometry 

Black holes exist in non-commutative geometry theories, as discussed in Sec. |2.3.5[ What is more, 
the usual Schwarzschild and Kerr solutions of General Relativity persist in these theories. This 
is not because such solutions have vanishing Weyl tensor, but because the quantity y^^C^aup 
happens to vanish for such metrics. Similarly, one would expect that the two-body, post-Newtonian 
metric that describes a binary black hole system should also satisfy the non-commutative geometry 
field equations, although this has not been proven explicitly. Similarly, although neutron star 
spacetimes have not yet been considered in non-commutative geometries, it is likely that if such 
spacetimes are stationary and satisfy the Einstein equations, they will also satisfy the modified 
field equations. Much more work on this is still needed to establish all of these concepts on a firmer 
basis. 

Gravitational waves exist in non-commutative gravity. Their generation for a compact binary 
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system in a circular orbit was analyzed by Nelson et al, in |293l I292J . They began by showing 
that a transverse-traceless gauge exists in this theory, although the transverse-traceless operator is 
slightly different from that in General Relativity. They then proceeded to solve the modified field 



equations for the metric perturbation [Eq. (42)] via a Green's function approach: 

/Afl 
,,^ ,,,, , J \i)UN(f~v?~V?) , (150) 

where recall that /3^ = (— 327rao)~^ acts like a mass term, the integral is taken over the entire past 
light cone, di{-) is the Bessel function of the first kind, \r\ is the distance from the source to the 
observer and the quadrupole moment is defined as usual: 
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where T*'^ is the time-time component of the matter stress-energy tensor. Of course, this is only 
the first term in an infinite multipole expansion. 



Although the integral in Eq. (|150| has not yet been solved in the post-Newtonian approximation 

2n 



Nelson et al |2931 1292] did solve for its time derivative to find 

n 
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where fi — 2'nF is the orbital angular frequency and we have defined 
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and one has assumed that the binary is in the x-y plane and the observer is on the z-axis. If 
one expands these expressions about /3 = oo however, one recovers the General Relativity solution 
to leading order, plus corrections that decay faster than l/r. This then automatically implies 
that such modifications to the generation mechanism will be difficult to observe for sources at 
astronomical distances. 

Given such a solution, one can compute the fiux of energy carried by gravitational waves to 
spatial infinity. Stein and Yunes [361j have shown that in quadratic gravity theories, this flux is 
still given by 

E ^ '^ f dnr^ (h^.h^'^) , (155) 

where /i^^ is the trace-reversed metric perturbation, the integral is taken over a 2-sphere at spatial 
infinity, and we recall that the angle brackets stand for an average over several wavelengths. Given 



the solution in Eq. (152 1, one finds that the energy flux is 
9 
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The asymptotic expansion of the term in between square brackets about /? = cx) is 
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which then leads to an energy flux identical to that in General Relativity, as any subdoniinant 
term goes to zero when the 2-sphere of integration is taken to spatial infinity. In that case, there 
are no modifications to the rate of change of the orbital frequency. Of course, if one were not to 
expand about /3 = oo, then the energy flux would lead to certain resonances at /3 = 2i7, but the 
energy flux is only well-defined at future null infinity. 

The above analysis was used by Nelson, et al J 2931 1292] to compute the rate of change of the 
orbital period of binary pulsars, in the hopes of using this to constrain /3. Using data from the 
double binary pulsar, they stipulated an order-of-magnitude constraint of /3 > 10~^^ m^^. Such 
an analysis, however, could be revisited to relax a few assumptions used in [29311292] . First, binary 
pulsar constraints on modified gravity theories require the use of at least three observables. These 
observables can be, for example, the rate of change of the period P, the line of nodes (l and the 
perihelion shift w. Any one observable depends on the parameters {rrii, 7712) in General Relativity 
or (?7ii,m2,/3) in non-commutative geometries, where mi. 2 are the component masses. Therefore, 
each observable corresponds to a surface of co-dimension one, ie. a two-dimensional surface or sheet 
in the three-dimensional space (mi,TO2,/3). If the binary pulsar observations are consistent with 
Einstein's theory, then all sheets will intersect at some point, within a certain uncertainty volume 
given by the observational error. The simultaneous fitting of all these observables is what allows one 
to place a bound on /3. The analysis of [2931 1292] assumed that all binary pulsar observables were 
known, except for /3, but degeneracies between (7711,7712,/?) could potentially dilute constraints on 
these quantities. Moreover, this analysis should be generalized to eccentric and inclined binaries, 
since binary pulsars are known to not be on exactly circular orbits. 

But perhaps the most important modification that ought to be made has to do with the 



calculation of the energy fiux itself. The expression for E in Eq. ( 155 1 in terms of derivatives of the 
metric perturbation derives from the effective gravitational wave stress-energy tensor, obtained by 
perturbatively expanding the action or the field equations and averaging over several wavelengths 
(the so-called Isaacson procedure |2121 I213J ). In modified gravity theories, the definition of the 
effective stress-energy tensor in terms of the metric perturbation is usually modified, as found for 
example in 1 3611 . In the case of non-commutative geometries. Stein and Yunes ^361^ showed that 



Eq. (155) still holds, provided one considers fiuxes at spatial infinity. The analysis of |2931 1292] . 
however, evaluated this energy flux at a fixed distance, instead of taking the r — ^ cxd limit. 

The balance law relates the rate of change of a binary's binding energy with the gravitational 
wave flux emitted by the binary, but for it to hold, one must require the following: (i) that the 
binary be isolated and it possesses a well-defined binding energy; (ii) the total stress-energy of the 
spacetime satisfies a local covariant conservation law. If (ii) holds, one can use this conservation law 
to relate the rate of change of the volume integral of the energy density, ie. the energy fiux, to the 
volume integral of the current density, which can be rewritten as an integral over the boundary of 
the volume through Stokes' theorem. Since in principle one can choose any integration volume, any 
physically meaningful result should be independent of the surface of that volume. This is indeed 
the case in General Relativity, provided one takes the integration 2-sphere to spatial infinity. 
Presumably, if one included all the relevant terms in E, without taking the limit to i*^, one would 
still find a result that is independent of the surface of this two-sphere. This, however, has not yet 
been verified. Therefore, the analysis of |293[I292 ] should be taken as an interesting first step toward 
understanding possible changes in the gravitational wave metric perturbation in non-commutative 
geometries. 

Not much beyond this has been done regarding non-commutative geometries and gravitational 
waves. In particular, one lacks a study of whether gravitational wave propagation is modified, what 
the final response function is, which of course depends on the time-evolution of all propagating 
gravitational wave degrees of freedom, and whether there are only the two usual dynamical degrees 
of freedom in the metric perturbation. 
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5.3 Generic Tests 

5.3.1 Massive Graviton Theories and Lorentz Violation 

Several massive graviton theories have been proposed to later be discarded due to ghosts, non- 
linear or radiative instabilities. Thus, little work has gone into studying whether black holes and 
neutron stars in these theories persist and are stable, and how the generation of gravitational waves 
is modified. Such questions will depend on the specific massive gravity model considered, and of 
course, if a Vainshtein mechanism is employed, then there will not be any modifications. 

A few generic properties of such theories, however, can still be stated. One of them is that the 
non-dynamical (near-zone) gravitational field will be corrected, leading to Yukawa-like modifica- 
tions to the gravitational potential |398) 

VuGir) = — e-'-/^« , or VuGir) = — (l + juGe^^^^A , (158) 

r r \ / 

where r is the distance from the source to a field point. The latter parameterization, for example, 
arises in gravitational theories with compactified extra dimensions |232j . Such corrections lead to a 
so-called fifth force, which then in turn allows us to place constraints on irig through Solar System 
observations |365j . Nobody has yet considered how such modifications to the near-zone metric 
could affect the binding energy of compact binaries and their associated gravitational waves. 

Another generic consequence of a graviton mass is the appearance of additional propagating 
degrees of freedom in the gravitational wave metric perturbation. In particular, one expects scalar, 
longitudinal modes to be excited (see eg. |133j ). This is, for example, the case if the action is of 
Pauli-Fierz type [1541 1133] . Such longitudinal modes arise due to the non- vanishing of the ^2 and 
^^3 Newman-Penrose scalars, and can be associated with the presence of spin-0 particles, if the 
theory is of Type N in the E{2) classification [555] . The specific form of the scalar mode will 
depend on the structure of the modified field equations, and thus, it is not possible to generically 
predict its associated contribution to the response function. 

A robust prediction of massive graviton theories relates to how the propagation of gravitational 
waves is affected. If the graviton has a mass, its velocity of propagation will differ from the speed of 



light, as given for example in Eq. (23). Will |398j showed that such a modification in the dispersion 
relation leads to a correction in the relation between the difference in time of emission Ate and 
arrival Ata of two gravitons: 



At a = {1 + Z) 



At, 



D f I 1 



(159) 



where z is here the redshift, Xg is the graviton's Compton wavelength, /e and f, are the emission 
frequencies of the two gravitons and D is the distance measure 

1 + z /"^ dz' 

"^ffcTio (1 + z')'[^m(1 + ;^')' + ^a]^/'' 

where Hq is the present value of the Hubble parameter, i7j\/ is the matter energy density and J7a 
is the vacuum energy density (for a zero spatial-curvature Universe). 

Even if the gravitational wave at the source is unmodified, the graviton time delay will leave an 
imprint on the Fourier transform of the response function by the time it reaches the detector |398) . 
This is because the Fourier phase is proportional to 

* (X 2^ / [t{f) - t,]df' , (161) 

■Jf. 
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where now t is not a constant but a function of frequency as given by Eq. ( 159 1 . Carrying out the 



integration, one finds that the Fourier transform of the response function becomes 

huG = /iGRe^^"'^"''"'' , (162) 

where /igr is the Fourier transform of the response function in General Relativity, we recall that 
u — {nMcfY^^ S'^^d "^s have defined 

A2(l + z) 

Such a correction is of 1 post-Newtonian order relative to the leading-order, Newtonian term in 
the Fourier phase. Notice also that there are no modifications to the amplitude at all. 

Numerous studies have considered possible bounds on Xg. The most stringent Solar System 
constraint is Xg > 2.8 x 10^^ km and it comes from observations of Kepler's third law (mainly Mars' 



orbit), which if the graviton had a mass would be modified by the Yukawa factor in Eq. (158 1. 
Observations of the rate of decay of the period in binary pulsars [1561 H5] can also be used to place 
the more stringent constraint A > 1.5 x 10"'^^ km. Gravitational wave observations of binary systems 
could also be used to constrain the mass of the graviton once gravitational waves are detected. One 
possible test is to compare the times of arrival of coincident gravitational wave and electromagnetic 
signals, for example in white-dwarf binary systems. Larson and Hiscock |251j and Gutler et al |119j 
estimated that one could constrain Ag > 3 x lO^'^ km with classic LISA. Will |398j was the first to 
consider constraints on Xg from gravitational wave observations only. He considered sky-averaged, 
quasi-circular inspirals and found that LIGO observations of lOAf© equal-mass black holes would 
lead to a constraint of Ag > 6 x 10^^ km with a Fisher analysis. Such constraints are improved 
to Ag > 6.9 X 10^^ km with classic LISA observations of IO'^Mq, equal-mass black holes. This 
increase comes about because the massive graviton correction accumulates with distance traveled 



(see Eq. ( 162 )). Since classic LISA would have been able to observe sources at Gpc scales with high 
signal-to-noise ratio, its constraints on Ag would have been similarly stronger than what one would 
achieve with LIGO observations. Will's study was later generalized by Will and Yunes [401], who 
considered how the detector characteristics affected the possible bounds on Ag. They found that 
this bound scales with the square-root of the LISA arm length and inversely with the square root of 
the LISA acceleration noise. The initial study of Will was then expanded by Berti et al [SBj, Yagi 
and Tanaka [411j, Arun and Will [Mj; Stavridis and Will ^360j and Berti, et al [S^ to allow for non 
sky-averaged responses, spin-orbit and spin-spin coupling, higher-harmonics in the gravitational 
wave amplitude, eccentricity and multiple detections. Although the bound deteriorates on average 
for sources that are not optimally oriented relative to the detector, the bound improves when one 
includes spin couplings, higher harmonics, eccentricity, and multiple detections as the additional 
information and power encoded in the waveform increases, helping to break parameter degeneracies. 
All of these studies, however, neglected the merger and ringdown phases of the coalescence, an 
assumption that was relaxed by Keppel and Ajith |233j . leading to the strongest projected bounds 
Ag > 4 X 10^^ km. Moreover, all studies until then had computed bounds using a Fisher analysis 
prescription, an assumption relaxed by del Pozzo et al |127j . who found that a Bayesian analysis 
with priors consistent with Solar System experiments leads to bounds stronger than Fisher ones by 
roughly a factor of two. All of these results are summarized in Table [21 normalizing everything to 
signal-to-noise ratio of 10. In summary, projected constraints on Ag are generically stronger than 
current Solar System or binary pulsar constraints by several orders of magnitude, given a LISA 
observation of massive black hole mergers. Even an Adv. LIGO observation would do better than 
current Solar System constraints by a factor between a few [127 to an order of magnitude |233| , 
depending on the source. 

Before proceeding, we should note that the correction to the propagation of gravitational waves 
due to a non-zero graviton mass are not exclusive to binary systems. In fact, any gravitational wave 
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Table 2: Comparison of Proposed Tests of Massive Graviton Theories. 



Reference 


Binary Mass 


Ag[10i5km] 


Properties 


i365^ 


X 


0.0028 


Solar System Dynamics 


[398 


(1O,1O)M0 


0.006 


LIGO, Fisher, Ang. Ave. 
circular, non-spinning 


PH] 


{W,W)Mq 


69 


LISA, Fisher, Ang. Ave. 
circular, non-spinning 


USB EH] 


(O.5,O.5)M0 


0.03 


LISA, WD- WD, coincident 
with electromagnetic signal 


mn 


{10',10'')Mq 


50 


LISA, Fisher, Ang. Ave. 
circular, non-spinning 


m 


{10^,10^)Mq 


10 


LISA, Fisher, Monte-Carlo, 
circular, w/spin-orbit 


m 


(10^10'^)Mo 


10 


LISA, Fisher, Ang. Ave 
higher-harmonics, circular, non-spinning 


mn 


(1O^1O''')M0 


22 


LISA, Fisher, Monte-Carlo 
eccentric, spin-orbit 


ma 


(lo^looAf0 


2.4 


DECIGO, Fisher, Monte-Carlo 
eccentric, spin-orbit 


|360| 


(10^106)Mo 


50 


LISA, Fisher, Monte-Carlo 
circular, w/spin modulations 


1233) 


{10',10'')Mq 


400 


LISA, Fisher, ang-averaged 
circular, non-spinning, w/merger 


\m\ 


(13,3)Af0 


0.006-0.014 


LIGO, Bayesian, ang-averaged 
circular, non-spinning 


[64 


(13,3)M0 


30 


cLISA, Fisher, Monte-Carlo 
multiple detections, circular, non-spinning 
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that propagates a significant distance from the source will suffer from the time delays described in 
this section. Binary inspirals are particularly useful as probes of this effect because one knows the 
functional form of the waveform, and thus, one can employ matched filtering to obtain a strong 
constraint. But in principle, one could use gravitational wave bursts from supernovae or other 
sources. 

We have so far concentrated on massive graviton theories, but, as discussed in Sec. |2.3.2[ 
there is a strong connection between such theories and Lorentz-violation. Modifications to the 
dispersion relation are usually a result of a modification of the Lorentz group or its action in real 
or momentum space. For this reason, it is interesting to consider generic Lorentz- violating-inspired. 



modified dispersion relations of the form of Eq. ( 24 1 , or more precisely |284j 

JL 

r,2 



^1 = 1- AE°'^^-'^ , (164) 

c^ 



where q;lv controls the structure of the modification and A its amplitude. When aLv = and 



A — rn^c^ one recovers the standard modified dispersion relation of Eq. ([23|). Equation (|164|) 
introduces a generalized time delay between subsequent gravitons of the form 



Ma = {l + Z) 



Atp 



Da.^ 



2\2 — ctLV y j'2— qlv fl2—a-i_ 



(165) 



where we have defined \a = hpA^^^""^^ '^\ with hp Planck's constant, and the generalized distance 
measure |284| 

Such a modification then leads to the following correction to the Fourier transform of the response 
function [284] 

hm = /iGRe^^-"'"" , (167) 

where /igr is the Fourier transform of the response function in General Relativity and we have 
defined [^ 

fl"^v#i^_I ±l£^l^_l!k^ , 6"{;^^i = 3(aLV - 1) . (168) 

The case q;lv = 1 is special leading to the Fourier phase correction |284j 

5*„,v=i-^ln7.. (169) 

The reason for this is that when qlv ~ 1 the Fourier phase is proportional to the integral of 1//, 
which then leads to a natural logarithm. 

Different q;lv limits deserve further discussion here. When aLV — 0, one of course recovers 
the standard massive graviton result with the mapping \^^ — > \^^ + X. . When aLv = 2, the 



dispersion relation is identical to that in Eq. ( 23 ) , but with a redefinition of the speed of light 



and should thus be unobservable. Indeed, in this limit the correction to the Fourier phase in 



Eq. ( 167 1 becomes linear in frequency, and this is 100% degenerate with the time of coalescence 
parameter in the standard General Relativity Fourier phase. Finally, relative to the standard 
General Relativity terms that arise in the post-Newtonian expansion of the Fourier phase, the 
new corrections are of (1 -I- 3aLv/2) post-Newtonian order. Then, if LIGO gravitational wave 
observations were incapable of discerning between a 4 post-Newtonian and a 5 post-Newtonian 
waveform, then such observations would not be able to see the modified dispersion effect if qlv > 2. 
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Mirshekari et al |284) confirmed this expectation with a Fisher analysis of non-spinning, comparable 
mass quasi-circular inspirals. They found that for qlv = 3, one can place very weak bounds on 
A^, namely A < 10~^ eV~^ with a LIGO observation of a (1.4, 1A)Mq neutron star inspiral, 
A < 0.2 eV~^ with an enhanced-LISA or NGO observation of a (10^, 10^)Mq black hole inspiral, 
assuming a signal-to-noise ratio of 10 and 100 respectively. A word of caution is due here, though, 
as these analysis neglect any Lorentz- violating correction to the generation of gravitational waves, 
including the excitation of additional polarization modes. One would expect that the inclusion of 
such effects would only strengthen the bounds one could place on Lorentz-violating theories, but 
these must be done on a theory by theory basis. 

5.3.2 Variable G Theories and Large Extra-Dimensions 

The lack of particular Lagrangian associated with variable G theories or extra dimensions makes 
it difficult to ascertain whether black hole or neutron star binaries exist in such theories. Whether 
this is so will depend on the particular variable G model considered. In spite of this, if such binaries 
do exist, the gravitational waves emitted by such systems will carry some generic modifications 
relative to the General Relativity expectation. 

Most current tests of the variability of Newton's gravitational constant rely on electromagnetic 
observations of massive bodies, such as neutron stars. As discussed in Sec. |2.3.4[ scalar-tensor 
theories can be interpreted as variable-G theories, where the variability of G is really a varia- 
tion in the coupling between gravity and matter. Newton's constant, however, serves the more 
fundamental role of defining the relationship between geometry or length and energy, and such a 
relationship is not altered in most scalar-tensor theories, unless the scalar fields are allowed to vary 
on a cosmological scale (background, homogeneous scalar solution). 

For this reason, one might wish to consider a possible temporal variation of Newton's constant 
in pure vacuum spacetimes, such as in binary black hole inspirals. Such temporal variation would 
encode G/G at the time and location of the merger event. Thus, once a sufficiently large number 
of gravitational wave events has been observed and found consistent with General Relativity, one 
could reconstruct a constraint map that bounds G/G along our past light-cone (as a function of 
redshift and sky position) . Since our past-light cone with gravitational waves would have extended 
to roughly redshift 10 with classic LISA (limited by the existence of merger events at such high 
redshifts) , such a constraint map would have been much more complete than what one can achieve 
with current tests at redshift almost zero. Big Bang nucleosynthesis constraints also allow us to 
bound a linear drift in G/G from z ^ 10"^ to zero, but these become degenerate with limits on 
the number of relativistic species. Moreover, these bounds exploit the huge lever-arm provided 
by integrating over cosmic time, but they are insensitive to local, oscillatory variations of G with 
periods much less than the cosmic observation time. Thus, gravitational-wave constraint maps 
would test one of the pillars of General Relativity: local position invariance. This principle (encoded 
in the equivalence principle) states that the laws of physics (and thus the fundamental constants 
of nature) are the same everywhere in the Universe. 

Let us then promote G to a function of time of the form |429j 

G(i,x,2/,z)«Gc + Gc(tc-i) , (170) 

where Gc — G{tc,Xc,yc,Zc) and Gc — {dG/dt){tc,Xc,yc, Zc) are constants, and the sub-index 
c means that these quantities are evaluated at coalescence. Clearly, this is a Taylor expansion 
to first order in time and position about the coalescence event (tcx].)^ which is valid provided 
the spatial variation of G is much smaller than its temporal variation, ie. |V*G| <C G, and the 
characteristic period of the temporal variation is longer than the observation window (at most, 
^^Dbs < 3 years for classic LISA), so that GcTohs <C Gc- Similar parameterization of G(i) have been 
used to study deviations from Newton's second law in the Solar System |134[ 13911 [388. ,371j . One 
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can thus think of this modification as the consequence of some effective theory that could represent 
the predictions of several different alternative theories. 

The promotion of Newton's constant to a function of time changes the rate of change of the 
orbital frequency, which then directly impacts the gravitational wave phase evolution. To leading 
order, Yunes, et al |429) find 

F = Fgr + ^U-^xS'^'iGMc) , (171) 

where -Fgr is the rate of change of the orbital frequency in General Relativity, due to the emission 
of gravitational waves and x — (27rMF)^/'^. Such a modification to the orbital frequency evolution 
leads to the following modification |429j to the Fourier transform of the response function in the 
stationary-phase approximation [5TJ 11181 11381 1418| 

h = hcR (1 + a^u'^o) e''9G"'<* , (172) 

where we recall again that u = (ttMc/)^/"^ and have defined the constant parameters |429j 

«G = -^^(GcM,), /3^ = ___^(G,M,), a = -8, 6 = -13, (173) 

to leading-order in the post-Newtonian approximation. We note that this corresponds to a cor- 
rection of —4 post-Newtonian order in the phase, relative to the leading-order term, and that the 
corrections are independent of the symmetric mass ratio, scaling only with the redshifted chirp 
mass M^. Due to this, one expects the strongest effects to be seen in low- frequency gravitational 
waves, such as those one could detect with LISA or DECIGO/BBO. 

Given such corrections to the gravitational wave response function, one can investigate the 
level to which a gravitational wave observation consistent with General Relativity would allow us 
to constrain G^. Yunes, et al [429] carried out such a study and found that for comparable-mass 
black hole inspirals of total redshifted mass rUz = 10^ Mq with LISA, one could constrain Gc/Gc ^ 
10~^ yr""'^ or better to redshift 10 (assuming a signal-to-noise ratio of 10'^). Similar constraints are 
possible with observations of extreme mass-ratio inspirals. The constraint is strengthened when 
one considers intermediate mass black hole inspirals, where one would be able to achieve a bound 
of Gc/Gc ^ 10""'^^ yr~^. Although this is not as stringent as the strongest constraints from other 



observations (see Sec. 2.3.4), we recall that gravitational waves constraints would measure local 
variations at the source, as opposed to local variations at zero redshift or integrated variations 
from the very early Universe. 

The effect of promoting Newton's constant to a function of time is degenerate with several 
different effects. One such effect is a temporal variability of the black hole masses, ie. if m 7^ 0. 
Such time-variation could be induced by gravitational leakage into the bulk in certain brane-world 



scenarios |226j . as we explained in Sec. 2.3.4 For a black hole of mass M, the rate of black hole 
evaporation is given by 



dM 
~dF 



■«'<'°-'(^)^(jo^)'«---'- ("*) 



where £ is the size of the large extra dimension. As expected, such a modification to a binary 
black hole inspiral will lead to a correction to the Fourier transform of the response function that 



is identical in structure to that of Eq. (172), but the parameters {/S^jb^) — > (/?ed,&ed) change 

to [HD] 

/3ed = -8.378 X 10-8 fJ-\ fi-^r] + 34r^A , 5ed = -13 . (175) 
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A similar expression is found for a neutron star/black-hole inspiral, except that the ?7-dependent 
factor in between parenthesis is corrected. 

Given a gravitational wave detection consistent with General Relativity, one could then in 
principle place an upper bound on i. Yagi et al 14101 carried out a Fisher analysis and found 
that a 1-year LISA detection would constrain i < 10^ /im with a (10, 1O^)M0 binary inspiral at 
signal-to-noise ratio of 100. This constraint is roughly two orders of magnitude weaker than current 
table-top experiment constraints ^. Moreover, the constraint weakens somewhat for more generic 
inspirals, due to degeneracies between £ and eccentricity and spin. A similar observation with the 
third generation detector DECIGO/BBO, however, should be able to beat current constraints by 
roughly one order of magnitude. Such a constraint could be strengthened by roughly one order 
of magnitude further, if one included the statistical enhancement in parameter estimation due to 
detection of order 10^ sources by DECIGO/BBO. 

Another way to place a constraint on £ is to consider the possibility of floating orbits |276) . 
Such orbits occur if the rate of change of the semi-major axis due to gravitational wave emission is 
balanced by the counter-acceleration induced by £. The latter arises because the mass loss corrects 
the rate of change of the semi-major axis o by d = —{M/M)a, due to conservation of specific 
orbital angular momentum. Therefore, there is a critical semi-major axis a^ at which the two 
forces balance and the binary stagnates. McWilliams [276] has argued that if a gravitational wave 
is detected at a given semi-major axis, then £ must be such that Oc is smaller than the observed 
semi-major axis. This would in principle place a constraint on £ of about 0(20 /im). Of course, 
such stagnate orbits would be rather strange, as if one takes this concept to the extreme and 
allows £ to be such that Uc is outside the innermost stable circular orbit, one would end up with a 
collection of black holes all sitting in a ring at a fixed semi-major axis. There are no astrophysical 
observations to support this idea, and given gravitational wave observations, one could use this to 
rule out a certain region of £ space. Yagi, et al f^lO improved and corrected this analysis to find 
that such a constraint is weaker than what one could achieve via matched filtering with a waveform 



of the form in Eq. (172), using the DECIGO detector. 

The G correction to the gravitational wave phase evolution is also degenerate with cosmological 
acceleration. That is, if a gravitational wave is generated at high-redshift, its phase will be affected 
by the acceleration of the Universe. To zeroth-order, the correction is a simple redshift of all 
physical scales. However, if one allows the redshift to be a function of time 

z^Zc + z,{t - Q ^z, + Ho \{1 + zcf - (1 + z,f'^ ^^]^f 1 (i - t,) , (176) 



then the observed waveform at the detector becomes structurally identical to Eq. (172) but with 

the parameters 

25 
/3i = ^^^icM,, 6i = -13. (177) 

Using the measured values of the cosmological parameters from the WMAP analysis |2401 1141) , 
however, one finds that this effect is roughly 10"'^ times smaller than that of a possible G correction 
at the level of the possible bounds quoted above |429j . Of course, if one could in the future constrain 
G better by 3 orders of magnitude, possible degeneracies with i would become an issue. 

A final possible degeneracy arises with the effect of a third body |421j , accretion disk migra- 
tion |237l I424J and the interaction of a binary with a circumbinary accretion disk [200j . All of these 
effects introduce corrections to the gravitational wave phase of negative PN order, just like the 
effect of a variable gravitational constant. Degeneracies of this type, however, are only expected to 
affect a small subset of binary black hole observations, namely those with a third body sufficiently 
close to the binary, or a sufficiently massive accretion disk. 
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5.3.3 Parity Violation 



As discussed in Sec. 2.3.6 the simplest action to model parity violation in the gravitational inter- 
action is given in Eq. ( 44 ) . Black holes and neutron stars exist in this theory, albeit non-rotating. 
A generic feature of this theory is that parity violation imprints onto the propagation of gravita- 
tional waves, an effect that has been dubbed amplitude birefringence. Such birefringence is not 
to be confused with optical or electromagnetic birefringence, in which the gauge boson interacts 
with a medium and is doubly-refracted into two separate rays. In amplitude birefringence, right- 
(left)-circularly polarized gravitational waves are enhanced or suppressed (suppressed or enhanced) 
relative to the General Relativity expectation as they propagate |216|. I265| [9l 14221 IHl I425J . 

One can understand amplitude birefringence in gravitational wave propagation due to a possible 
non-commutativity of the parity operator and the Hamiltonian. The Hamiltonian is the generator 
of time evolution, and thus, one can write |425j 



where / is the gravitational wave angular frequency, t is time, and h^^x,k are the gravitational wave 
Fourier components with wavenumber k. The quantity Uc models possible background curvature 
effect, with Uc = 1 for propagation on a Minkowski metric, and v proportional to redshift for 
propagation on a Friedman- Robertson- Walker metric j247j . The quantity v models possible parity- 
violating effects, with u = in General Relativity. One can rewrite the above equation in terms of 
right- and left-circular polarizations, /ir.l = (ft-+ ± i/ix)/v2 to find 



HkW \ ^ -^ft { U,+V 0\ f ;iR.k(0) 

/lL,k(i) J \ Uc-V J \ /lL,k(0) 



(179) 



Amplitude birefringence has the effect of modifying the eigenvalues of the diagonal propagator 
matrix for right and left-polarized waves, with right modes amplified or suppressed and left modes 
suppressed or enhanced relative to General Relativity, depending on the sign of w. In addition 
to this parity-violating propagation effects, parity violation should also leave an imprint in the 
generation of gravitational waves. Such effects, however, need to be analyzed on a theory by 
theory basis. Moreover, the propagation distance-independent nature of generation effects should 
make them easily distinguishable from the propagation effects we consider here. 

The degree of parity violation, v, can be expressed entirely in terms of the waveform observable 
viapS] 



1 / hn ft.L 



2V/JGR ;,GR;-i(^'^L-<5</.R,), (180) 

where h^^ is the General Relativity expectation for a right- or left-polarized gravitational wave. 
In the last equality we have also introduced the notation d(f) = (p — (p'^^, where t/)^^ is the General 
Relativity gravitational wave phase and 
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/^o.R.Le-^t^^"'-""^'] , (181) 



where /iq.r.l is a constant factor, k is the conformal wave number and (?7, x*) are conformal 
coordinates for propagation in a Friedmann-Robertson- Walker Universe. The precise form of v 
will depend on the particular theory under consideration. For example, in non-dynamical Chern- 
Simons gravity with a field ■& = i?(t), and in an expansion about 2; <C 1, one finds |425j 

V = %fz Uo - ^ ) = ^^fD {Ho^o - ^0) , (182) 
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where 'do is the Chern-Simons scalar field at the detector, with a the Chern-Simons coupling 



constant [see e.g. Eq. (44|], z is redshift, D is the comoving distance and Hq is the value of the 



Hubble parameter today and / is the observed gravitational wave frequency. When considering 
propagation on a Minkowski background, one obtains the above equation in the limit as a — > 0, so 
the second term dominates, where a is the scale-factor. To leading-order in a curvature expansion, 
the parity-violating coefficient v will always be linear in frequency, as shown in Eq. 



( |182| ). For 

more general parity violation and flat-spacetime propagation, v will be proportional to {fD)f°'a, 
where a is a coupling constant of the theory (or a certain derivative of a coupling field) with units 
of [Length]" (in the previous case, a = 0, so the correction was simply proportional to fDa, where 
a oc -d). 

How does such parity violation affect the waveform? By using Eq. (179) one can easily show 



that the Fourier transform of the response function becomes [HI I422[ I425f 

hpv = {F+ - i V Fy,)h+ + {Fy, + i V F+)hy 



(183) 



One can of course rewrite this in terms of a real amplitude correction and a real phase correction. 
Expanding in i; <C 1 to leading order, we find [425] 
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where /igr is the Fourier transform of the response function in General Relativity and we have 
defined 



Qgr = Vn (1 + cos2 + 4cos2 iFi 
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(185) 
(186) 

(187) 



We see then that amplitude birefringence modifies both the amplitude and the phase of the response 



Eq. ([T84| as [425] 

where we have defined the coefficients 



function. Using the non-dynamical Chern-Simons expression for v in Eq. (182), we can rewrite 

(188) 
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6pv = 6 , 
(190) 



where we recall that u = {irJAcfY^^ ■ The phase correction corresponds to a term of 5.5 post- 
Newtonian order relative to the Newtonian contribution, and, it scales quadratically with the 
Chern-Simons coupling field d, which is why it was left out in |425j . The amplitude correction, on 
the other hand, is of 1.5 post-Newtonian order relative to the Newtonian contribution. Since both 
of these appear as positive-order, post-Newtonian corrections, there is a possibility of degeneracy 
between them and standard waveform template parameters. 

Given such a modification to the response function, one can ask whether such parity violation 
is observable with current detectors. Alexander et al 1422) argued that a gravitational wave 
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observation with LISA would be able to constrain an integrated measure of v, because LISA 
can observe massive black hole mergers to cosmological distances, while amplitude birefringence 
accumulates with distance traveled. For such an analysis, one cannot Taylor expand {) about its 
present value, and instead, one finds that 



e^'^/^^^) , (191) 



where we have defined 



C(,) = ^ fd^il + zf^ 



^^^ ^'^^ (192) 
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We can solve the above equation to find 



1 + e2'r/C(z) 



2 dz dz 



TTfCiz), (193) 



where in the second equality we have linearized about v ^ I and /C ^ 1. Alexander et al [HI 
I422J realized that this induces a time-dependent change in the inclination angle (ie. the apparent 
orientation of the binary's orbital angular momentum with respect to the observer's line-of-sight), 
since the latter can be defined by the ratio /ir/Zil- They then carried out a simplified Fisher analysis 
and found that a LISA observation of the inspiral of two massive black holes with components 
masses 10^Mo(l-|-z)~-'^ at redshift z — 15 would allow us to constrain the integrated dimensionless 
measure C < 10"^^ to la. One might worry that such an effect would be degenerate with other 
standard General Relativity processes that induce similar time-dependencies, such as spin-orbit 
coupling. However, this time-dependence is very different from that of the parity-violating effect, 
and thus, Alexander et al [HI 1422) argued that these effects would be weakly correlated. 

Another test of parity violation was proposed by Yunes et al |425j . who considered the co- 
incident detection of a gravitational wave and a gamma-ray burst with the SWIFT |173| and 
GLAST/Fermi [SU] gamma-ray satellites, and the ground-based LIGO [5] and Virgo dj gravita- 
tional wave detectors. If the progenitor of the gamma-ray burst is a neutron-star/neutron-star 
merger, the gamma-ray jet is expected to be highly coUimated. Therefore, an electromagnetic 
observation of such an event implies that the binary's orbital angular momentum at merger must 
be pointing along the line of sight to Earth, leading to a strongly circularly-polarized gravitational 
wave signal and to maximal parity violation. If the gamma-ray burst observation were to provide 
an accurate sky location, one would be able to obtain an accurate distance measurement from the 
gravitational wave signal alone. Moreover, since GLAST/Fermi observations of gamma-ray bursts 
occur at low redshift, one would also possess a purely electromagnetic measurement of the distance 
to the source. Amplitude birefringence would manifest itself as a discrepancy between these two 
distance measurements. Therefore, if no discrepancy is found, the error ellipse on the distance 
measurement would allow us to place an upper limit on any possible gravitational parity violation. 
Because of the nature of such a test, one is here constraining generic parity violation over distances 
of hundreds of Mpc, along the light-cone on which the gravitational waves propagate. 

The coincident gamma-ray burst/gravitational wave test compares favorably to the pure LISA 
test, with the sensitivity to parity violation being about 2-3 orders of magnitude better in the 
former case. This is because, although the fractional error in the gravitational wave distance 
measurement is much smaller for LISA than for LIGO, since it is inversely proportional to the 
signal-to-noise ratio, the parity violating effect also depends on the gravitational wave frequency, 
which is much larger for neutron-star inspirals than massive black hole coalescences. Mathemati- 
cally, the simplest models of gravitational parity violation will lead to a signature in the response 
function that is proportional to the gravitational wave wavelengttjj Aqw oc Df. Although the 

^Even if it is not linear, the effect should scale with positive powers of Aqw It is difficult to think of any parity 
violating theory that would lead to an inversely proportional relation. 
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coincident test requires small distances and low signal-to-noise ratios (by roughly 1-2 orders of 
magnitude), the frequency is also larger by a factor of 5-6 orders of magnitude for the LIGO- Virgo 
network. 

The coincident gamma-ray burst/gravitational wave test also compares favorably to current 
Solar System constraints. Using the motion of the LAGEOS satellites, Smith et al. J351I have 
placed the la bound do < 2000 km assuming ??o = 0- A similar assumption leads to a 2ct bound 
of 1^0 < 200 km with a coincident gamma-ray burst/gravitational wave observation. Moreover, the 
latter test also allows us to constrain the second time-derivative of the scalar field. Finally, a LISA 
observation would constrain the integrated history of d along the past light-cone on which the 
gravitational wave propagated. These tests, however, are not as stringent as the recently proposed 
test by Dyda et al |143j . do < 10^^ km, assuming the effective theory cut-off scale is less than 
10 eV and obtained by demanding that the energy density in photons created by vacuum decay 
over the lifetime of the Universe not violate observational bounds. 

The coincident test is somewhat idealistic in that there are certain astrophysical uncertainties 
that could hamper the degree to which we could constrain parity violation. One of the most 
important uncertainties relates to our knowledge of the inclination angle, as gamma-ray burst jets 
are not necessarily perfectly aligned with the line of sight. If the inclination angle is not known 
a priori, it will become degenerate with the distance in the waveform template, decreasing the 
accuracy to which the luminosity could be extracted from a pure gravitational wave observation 
by at least a factor of two. Even after taking such uncertainties into account, Yunes et al |425j 
found that ^o could be constrained much better with gravitational waves than with current Solar 
System observations. 

5.3.4 Parameterized Post-Einsteinian Framework 

One of the biggest disadvantages of a top-down or direct approach toward testing General Relativity 
is that one must pick a particular theory from the beginning of the analysis. However, given the 
large number of possible modifications to Einstein's theory and the lack of a particularly compelling 
alternative, it is entirely possible that none of these will represent the correct gravitational theory 
in the strong field. Thus, if one carries out a top-down approach, one will be forced to make the 
assumption that we, as theorists, know which modifications of gravity are possible and which 
are not |428j . The parameterized post-Einsteinian (ppE) approach is a framework developed 
specifically to alleviate such a bias by allowing the data to select the correct theory of Nature 
through the systematic study of statistically significant anomalies. 

For detection purposes, one usually expects to use match filters that are consistent with General 
Relativity. But if General Relativity happened to be wrong in the strong field, it is possible that a 
General Relativity template will still extract the signal, but with the wrong parameters. That is, 
the best fit parameters obtained from a matched filtering analysis with General Relativity templates 
will be biased by the assumption that General Relativity is sufficiently accurate to model the entire 
coalescence. This fundamental bias could lead to a highly distorted image of the gravitational wave 
Universe. In fact, recent work by Vallisneri and Yunes |377) indicates that such fundamental bias 
could indeed be present in observations of neutron star inspirals, if General Relativity is not quite 
the right theory in the strong-field. 

One of the primary motivations for the development of the ppE scheme was to alleviate funda- 
mental bias, and one of its most dangerous incarnations: stealth-bias |115j . If General Relativity 
is not the right theory of nature, yet all our future detections are of low signal-to-noise ratio, we 
may estimate the wrong parameters from a matched-filtering analysis, yet without being able to 
identify that there is an non-General Relativity anomaly in the data. Thus, stealth bias is noth- 
ing but fundamental bias hidden by our limited signal-to-noise ratio observations. Vallisneri and 
Yunes |377j have found that such stealth-bias is indeed possible in a certain sector of parameter 



70 



space, inducing errors in parameter estimation that could be larger than statistical ones, without 
us being able to identify the presence of a non-General Relativity anomaly. 

5.3.4.1 Historical Development 

The ppE scheme was designed in close analogy with the parameterized post-Newtonian (ppN) 
framework, developed in the 1970s to test General Relativity with Solar System observations (see 
e.g. [399j for a review). In the Solar System, all direct observables depend on a single quantity, 
the metric, which can be obtained by a small- velocity/weak-field post-Newtonian expansion of the 
field equations of whatever theory one is considering. Thus, Will and Nordtvedt |298l 13931 14001 
12971 I394J proposed the generalization of the Solar System metric into a meta-metric that could 
effectively interpolate between the predictions of many different alternative theories. This meta- 
metric depends on the product of certain Green function potentials and ppN parameters. For 
example, the spatial-spatial components of the meta-metric take the form 

g,j = 5,, (1 + 27C/ + . . . ) (194) 

where 5ij is the Kronecker delta, U is the Newtonian potential and 7 is one of the ppN parameters, 
which acquires different values in different theories: 7 = 1 in General Relativity, 7 = (1-|-wbd)(2 + 
wbd)^^ ^ l^'^BD '^^ Jordan-Brans-Dicke-Fierz theory, etc. Therefore, any Solar System observable 
could then be written in terms of system parameters, such as the masses of the planets, and the 
ppN parameters. An observation consistent with General Relativity allows for a bound on these 
parameters, thus simultaneously constraining a large class of modified gravity theories. 

The idea behind the ppE framework was to develop a formalism that allowed for similar generic 
tests but with gravitational waves instead of Solar System observations. The first such attempt 
was by Arun, et al [23 I285| who considered the quasi-circular inspiral of compact objects. They 
suggested the waveform template family 
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This waveform depends on the standard system parameters that are always present in General 
Relativity waveforms, plus one theory parameter /3pnt that is to be constrained. The quantity 6pN 
is a number chosen by the data analyst and is restricted to be equal to one of the post-Newtonian 
predictions for the phase frequency exponents, ie. 6pN = (—5, —3, —2, —1, . . .). 



The template family in Eq. ( 195 1 allows for the so-called post-Newtonian tests of General 
Relativity, ie. consistency checks of the signal with the post-Newtonian expansion. For example, 
let us imagine that a gravitational wave has been detected with sufficient signal-to-noise ratio that 
the chirp mass and mass ratio have been measured from the Newtonian and 1 post-Newtonian 
terms in the waveform phase. One can then ask whether the 1.5 post-Newtonian term in the phase 
is consistent with these values of chirp mass and mass ratio. Put another way, each term in the 
phase can be thought of as a curve in (Mc,/;) space. If General Relativity is correct, all these 
curves should intersect inside some uncertainty box, just like when one tests General Relativity 
with binary pulsar data. From that standpoint, these tests can be thought of as null-tests of 
General Relativity and one can ask: given an event, is the data consistent with the hypothesis 
AppE = for the restricted set of frequency exponents fopN? 

A Fisher and a Bayesian data analysis study of how well /3pnt could be constrained given a cer- 
tain fepN was carried out in [2851 I21H I260J . Mishra et al [285[ considered the quasi-circular inspiral 
of non-spinning compact objects and showed that Ad. LIGO observations would allow one to con- 
strain /3pivix to 6% up to the 1.5 post-Newtonian order correction (5pN — ^2). Third-generation 
detectors, such as ET, should allow for better constraints on all post-Newtonian coefficients to 
roughly 2%. Clearly, the higher the value of 6pN, the worse the bound on /3pnt because the power 
contained in higher frequency exponent terms decreases, ie. the number of useful additional cycles 
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induced by the /^pnt"^^'^ term decreases as &pn increases. Huwyler et al [211] repeated this analysis 
but for LISA observations of the quasi-circular inspiral of black hole binaries with spin precession. 
They found that the inclusion of precessing spins forces one to introduce more parameters in the 
waveform, which dilutes information and weakens constraints on /3pnt by as much as a factor of 
5. Li et al |260j carried out a Bayesian analysis of the odds-ratio between General Relativity and 
restricted ppE templates given a non-spinning, quasi-circular compact binary inspiral observation 
with Ad. LIGO and Virgo. They calculated the odds ratio for each value of bp^ listed above and 
then combined all of this into a single probability measure that allows one to quantify how likely 
the data is to be consistent with General Relativity. 

5.3.4.2 The Simplest ppE Model 



One of the main disadvantage of the post-Newtonian template family in Eq. ( 195 1 is that 
it is not rooted on a theoretical understanding of modified gravity theories. To alleviate this 
problem, Yunes and Pretorius |428| re-considered the quasi-circular inspiral of compact objects. 
They proposed a more general ppE template family through generic deformations of the £ ^ 2 
harmonic of the response function in Fourier space : 
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where now (appE, appE,/3ppE, &ppe) are all free parameters to be fitted by the data, in addition 
to the usual system parameters. This waveform family reproduces all predictions from known 
modified gravity theories: when (appE,/3ppE) = (0,0), the waveform reduces exactly to General 
Relativity, while for other parameters one reproduces the modified gravity predictions of Table |3] 
In this Table, recall that 5* is the difference in the square of the sensitivities and wbd is the Brans- 



Table 3: Parameters that define the deformation of the response function in a variety of modified 
gravity theories. The notation • means that a value for this parameter is irrelevant, as its amplitude 
is zero. 



Theory 


appE 


flppE 


/3ppE 


foppE 


Jordan-Brans- 
Dicke-Fierz 


5 S^ 2/5 
96 wbd ' 


-2 


5 S^ 2/5 
3584 WBD ' 


-7 


Dissipative 
Einstein-Dilaton- 
Gauss-Bonnet 
gravity 







5 ^,„-18/5x2 


-7 


Massive Graviton 







Ai(l-K2) 


-3 


Lorentz Violation 







(l-7Lv) Af-^LV (H-z)i-"'LV 


-37LV-3 


G(t) Theory 


-^GM, 


-8 


25 /y TV4- 

65536 "-^cJVLc 


-13 


Extra Dimensions 






-25riml^^-n3-26ry + 2V) 


-13 


Non-Dynamical 

Chern-Simons 

Gravity 


QfpV 


3 


/3pv 


6 


Dynamical 

Chern-Simons 

Gravity 







^dCS 


-1 



Dicke coupling parameter (see Sec. 5.2.1 and we have here neglected the scalar mode), Cs is the 
coupling parameter in Einstein-Dilaton-Gauss-Bonnet theory (see Sec. 5.2.21, where we have here 
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included both the dissipative and the conservative corrections, D is a certain distance measure 
and Xg is the Conipton wavelength of the graviton (see Sec. 5.3.1), Alv is a distance scale at 



which Lorentz-violation becomes important and 7t,v i s the graviton momentum exponent in the 



deformation of the dispersion relation (see Sec. 5.3.1), Gc is the value of the time derivative of 
Newton's constant at coalescence and dM/dt is the mass loss due to enhanced Hawking radiation 
in extra-dimensional scenarios (see Sec. 5.3.2[), /3dcs is given in Eq. (148) and (apv;/^pv) are given 



in Eqs. (189) and (190) of Sec. 5.3.3 



Although there are only a few modified gravity theories where the leading-order post-Newtonian 
correction to the Fourier transform of the response function can be parameterized by post-Newtonian 



waveforms of Eq. ( 195 1, all such predictions can be modeled with the ppE templates of Eq. ( 196 1 



In fact, only massive graviton theories, certain classes of Lorentz- violating theories and dynamical 



Chern-Simons gravity lead to waveform corrections that can be parameterized via Eq. (195). For 



example, the lack of amplitude corrections in Eq. ( 195 ) does not allow for tests of gravitational 



parity violation or non-dynamical Chern-Simons gravity. 

This, however, does not imply that Eq. ( 196 ) can parameterize all possible deformations of 
General Relativity. First, Eq. (196) can be understood as a single-parameter deformation away 
from Einstein's theory. If the correct theory of Nature happens to be a deformation of General 
Relativity with several parameters (eg. several coupling constants, mass terms, potentials, etc. ), 
then Eq. ( 196 ) will only be able to parameterize the one that leads to the most useful-cycles. This 
was recently verified by Sampson et al [340]. Second, Eq. (196) assumes that the modification can 
be represented as a power series in velocity, with possibly non-integer values. Such an assumption 
does not allow for possible logarithmic terms, which are known to arise due to non-linear memory 
interactions at sufficiently high-post-Newtonian order. It also does not allow for interactions that 
are screened, eg. in theories with massive degrees of freedom. Nonetheless, the parameterization 



in Eq. (196) will still be able to signal that the detection is not a pure Einstein event, at the cost 



of biasing their true value. 



The inspiral ppE model of Eq. ( 196 ) is motivated not only from examples of modified gravity 



predictions, but from generic modifications to the physical quantities that drive the inspiral: the 
binding energy or Hamiltonian and the radiation-reaction force or the fluxes of the constants of the 
motion. Yunes and Pretorius |428l and Chatziioannou, et al. [95 considered generic modifications 
of the form 
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where (p, g) G Z, since otherwise one would lose analyticity in the limit of zero- velocities for 
circular inspirals, and where {A, B) are parameters that depend on the modified gravity theory 
and, in principle, could depend on dimensionless quantities like the symmetric mass ratio. Such 
modifications lead to the following corrections to the SPA Fourier transform of the t = 2 time- 
domain response function for a quasi-circular binary inspiral template (to leading order in the 
deformations and in post-Newtonian theory) 
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Of course, usually one of these two modifications dominates over the other, depending on whether 
q > p or p < q. In Jordan-Brans-Dicke-Fierz theory, for example, the radiation-reaction correction 
dominates as q < p. If, in addition to this modifications in the generation of gravitational waves, 
one also allows for modifications in the propagation, one is then lead to the following template 
family [95] 

^pptU,2 = ^ (ttM,/)-^/^ e'^*«- [l + c/3ppE (^M J)''--/^+^/3] e^'^pp-"'^^" g-PPB^''-- ^ (200) 

Here, (6ppE, /3ppE) and (fcppE, KppE) are ppE parameters induced by modifications to the generation 
and propagation of gravitational waves respectively, where still (6ppE,fcppE) G 2, while c is fully 
determined by the former set via 

16 (3-6)(42b + 61 + 5b^) 
'^™"''" 15 5&2+466 + 81 ' ^ ' 

if the modifications to the binding energy dominate, 

16 

Cdiss = --^(3-6)6, (202) 

15 

if the modifications to the energy flux dominate, or 

32 b{3-b){Ub + 71 + 5b^) 
'''°*^ - " 15 562+466 + 81 ■ ^^^^> 

if both corrections enter at the same post-Newtonian order. Noticing again that if only a single 
term in the phase correction dominates in the post-Newtonian approximation (or both will enter 



at the same post- Newtonian order), one can map Eq. (198) to Eq. (1961 by a suitable redefinition 
of constants. 

5.3.4.3 More Complex ppE Models 

Of course, one can introduce more ppE parameters to increase the complexity of the waveform 



family, and thus, Eq. (196) should be thought of as a minimal choice. In fact, one expects any 
modified theory of gravity to introduce not just a single parametric modification to the amplitude 
and the phase of the signal, but two new functional degrees of freedom: 

appE^""^ ^ SAppEiX"" , e° ; u) , /3ppEii'"^ ^ ^*ppE ( A" , 0" ; u) , (204) 

where these functions will depend on the frequency u, as well as on system parameters A° and 
theory parameters 0°. In a post-Newtonian expansion, one expects these functions to reduce to 



leading-order to the left-hand sides of Eq. (204), but also to acquire post-Newtonian corrections of 
the form 

MppE(A^ r ; u) = appE(A^ ^ )w"-- ^ a„^ppE(A^ ^ )u" , (205) 

n 

'5*ppE(A^ r ; u) = /3ppE(A^ r )u''pp- Y. f^n,ppE{X\ e'^)u- , (206) 

n 

where here the structure of the series is assumed to be of the form u" with m > 0. Such a model, 
also suggested by Yuncs and Prctorius [428J , would introduce too many new parameters that would 
dilute the information content of the waveform model. Recently, Sampson et al |340j demonstrated 



that the simplest ppE model of Eq. ( 196 ) suffices to signal a deviation from General Relativity, 



even if the injection contains three terms in the phase. In fact, this is precisely one of the most 
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important differences between the ppE and ppN frameworks: in the latter, it does not matter 
how many ppN parameters are introduced, but in the former, one wishes to minimize the number 
of theory parameters. One must then find a balance between the number of additional theory 
parameters to introduce and the amount of bias contained in the templates. 

At this junction, one must emphasize that frequency exponents in the amplitude and phase 
correction were above assumed to be integers, ie. (flppE, 6ppE, n) € Z. This must be the case if these 
corrections arise due to modifications that can be represented as integer powers of the momenta 
or velocity. We are not aware of any theory that predicts corrections proportional to fractional 
powers of the velocity for circular inspirals. Moreover, one can show that theories that introduce 
non-integer powers of the velocity into the equations of motion will lead to issues with analyticity 
at zero velocity and a breakdown of uniqueness of solutions [95]. In spite of this, modified theories 
can introduce logarithmic terms, that for example enter at high post-Newtonian order in General 
Relativity due to non-linear propagation effects (see eg. [SBj and references therein). Moreover, 
certain modified gravity theories introduce screened modifications that become "active" only above 
a certain frequency. Such effects would be modeled through a Heaviside function, for example 
needed when dealing with massive Brans-Dicke gravity ir32( i87 l fTS l I426J . Even this non-polynomial 
injections, however, would be detectable with the simplest ppE model. In essence, one finds similar 
results as if one were trying to fit a 3-parameter injection with the simplest 1-parameter ppE 
model [MP] . 

One can of course also generalize the inspiral ppE waveform families to more general orbits, 
for example through the inclusion of spins aligned or counter-aligned with the orbital angular 



momentum. More general inspirals would still lead to waveform families of the form of Eq. ( 196 ) 



or (2001, but where the parameters (appEj/3ppE) would now depend on the mass ratio, mass 
difference, and the spin parameters of the black holes. With a single detection, one cannot break 
the degeneracy in the ppE parameters and separately fit for its system parameter dependencies. 
Given multiple detections, however, one should be able to break such a degeneracy, at least to a 
certain degree |115j . Such breaking of degeneracies begins to become possible when the number of 
detections exceeds the number of additional parameters required to capture the physical parameter 
dependencies of (appE, /^pps)- 

PpE waveforms can be extended to account for the merger and ringdown phases of coalescence. 
Yunes and Pretorius have suggested the following template family to account also for this |428j 

ppE / < /iM, 

.tfuii = { 7«^e*(^+™) /iM < / < /mrd, (207) 

C 1+4^2^2 J(y_yj^„)d / > /mrd, 



'•ppEjfuU 



where the subscript IM and MRD stand for inspiral-merger and merger-ringdown. The merger 
phase (/iM < / < /mrd) is here modeled as an interpolating region between the inspiral and 
ringdown, where the merger parameters (7, S) are set by continuity and differentiability and the 
ppE merger parameters (c, e) should be fit for. In the ringdown phase (/ > /mrd), the response 
function is modeled as a single-mode generalized Lorentzian, with real and imaginary dominant 
frequencies /rd and t, ringdown parameter ( also set by continuity and differentiability, and the 
ppE ringdown parameters (k, d) are to be fit for. The transition frequencies (/im, /mrd) can either 
be treated as ppE parameters or set via some physical criteria, such as at light-ring frequency and 
the fundamental ringdown frequency, respectively. 

Recently, there has been effort to generalize the ppE templates to allow for the excitation of non- 
GR gravitational wave polarizations. Modifications to only the two General Relativity polarizations 
map to corrections to terms in the time-domain Fourier transform that are proportional to the £ = 2 
harmonic of the orbital phase. Arun suggested, however, that if additional polarizations are present, 
other terms proportional to the £ = and £ = 1 harmonic will also arise ,34j . Chatziioannou, Yunes 
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and Cornish [95j have found that the presence of such harmonics can be captured through the more 
complete single-detector template family 
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where we have defined Ui = (27rMc//^)^^'^. The ppE theory parameters are now 9 — (6ppE, /?ppE, fcppE, ^ppE, 7ppE, *&c )• 
Of course, one may ignore (fcppE, KppE) altogether, if one wishes to ignore propagation effects. Such 
a parameterization recovers the predictions of Jordan-Brans-Dicke-Fierz theory for a single-detector 
response function [33], as well as Arun's analysis for generic dipole radiation [33]. 

One might worry that the corrections introduced by the £ = 1 harmonic, ie. terms proportional 



to 7ppE in Eq. ( 208 1 , will be degenerate with post-Newtonian corrections to the amplitude of the 
i = 2 mode (not displayed in Eq. (208)). This, however, is clearly not the case, as the latter 
scale as (ttMc/)"^/^"''"'/'^ with n an integer greater than 0, while the i = 1 mode is proportional 
to (ttMc/)""^'^, which would correspond to a (—0.5) post-Newtonian order correction, ie. n = —1. 
On the other hand, the ppE amplitude corrections to the i = 2 mode, ie. terms proportional to 



/3ppE in the amplitude of Eq. (208), can be degenerate with such post-Newtonian corrections when 
6ppE is an integer greater than —4. 

5.3.4.4 Applications of the ppE Formalism 



The two models in Eq. ( 196 ) and ( 200 ) answer different questions. The latter contains a stronger 
prior (that ppE frequency exponents be integers), and thus, it is ideal for fitting a particular set 
of theoretical models. On the other hand, Eq. ( 196 ) with continuous ppE frequency exponents 



allows one to search for generic deviations that are statistically significant, without imposing such 



theoretical priors. That is, if a deviation from General Relativity is present, then Eq. ( 196 ) is more 
likely to be able to fit it, than Eq. (200 1. If one prioritizes the introduction of the least number 
of new parameters, Eq. (1961 with (appE,&ppE) G ^ can still recover deviations from General 



Relativity, even if the latter cannot be represented as a correction proportional to an integer power 
of velocity. 

Given these ppE waveforms, how should they be used in a data analysis pipeline? The main 
idea behind the ppE framework is to match filter or perform Bayesian statistics with ppE enhanced 
template banks to allow the data to select the best-fit values of 9'^. As discussed in |428|,I115] and 
then later in |260j , one might wish to first run detection searches with General Relativity template 
banks, and then, once a signal has been found, do a Bayesian model selection analysis with ppE 
templates. The first such Bayesian analysis was carried out by Cornish et al |115j . who concluded 
that an Ad. LIGO detection at signal-to-noise ratio of 20 for a quasi-circular, non-spinning black 
hole inspiral would allow us to constrain ctppE and /3ppE much better than existent constraints for 
sufficiently strong-field corrections, eg. 5ppE > —5. This is because for lower values of the frequency 
exponents, the corrections to the waveform are weak-field and better constrained with binary pulsar 
observations |423| . The large statistical study of Li, et al 1260' uses a reduced set of ppE waveforms 
and investigates our ability to detect deviations of GR when considering a catalogue of Advanced 
LIGO/Virgo detections. Of course, the disadvantage of such a pipeline is that it requires a first 
detection, and if the gravitational interaction is too different from General Relativity's prediction, 
it is possible that a search with General Relativity templates might miss the signal all together; 
we dim this possibility to be less likely. 

A built-in problem with the ppE and the ppN formalisms is that if a non-zero ppE or ppN 
parameter is detected, then one cannot necessarily map it back to a particular modified gravity 
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action. On the contrary, as suggested in Table [31 there can be more than one theory that predicts 
structurally similar corrections to the Fourier transform of the response function. For example, 
both Jordan-Brans-Dicke-Fierz theory and the dissipative sector of Einstein-Dilaton-Gauss-Bonnet 
theory predict the same type of leading-order correction to the waveform phase. However, if a given 
ppE parameter is measured to be non-zero, this could provide very useful information as to the 
type of correction that should be investigated further at the level of the action. The information 
that could be extracted is presented in Table |4j which is derived from knowledge of the type of 
corrections that lead to Table [3j Moreover, if a follow-up search is done with the ppE model in 

Table 4: Interpretation of non-zero ppE parameters. 



flppE 


&ppE 


Interpretation 


1 




Parity Violation 


-8 


-13 


Anomalous Acceleration, Extra 
Dimensions, Violation of Posi- 
tion Invariance 




-7 


Dipole Gravitational Radiation, 
Electric Dipole Scalar Radiation 




-3 


Massive Graviton Propagation 


a spin 


-1 


Magnetic Dipole Scalar Radia- 
tion, Quadrupole Moment Cor- 
rection, Scalar Dipole Force 



Eq. ( 200 1 , one could infer whether the correction is one due to modifications to the generation or 



the propagation of gravitational waves. In this way, a non-zero ppE detection could inform theorist 
of what type of General Relativity modification is preferred by Nature. 

5.3.4.5 Degeneracies 

Much care must be taken, however, to avoid confusing a ppE theory modification with some 
other systematic, such as an astrophysical, a mismodeling or an instrumental effect. Instrumental 
effects can be easily remedied by requiring that several instruments, with presumably unrelated 
instrumental systematics, independently derive a posterior probability for (ctppE, /3ppE) that peaks 
away from zero. Astrophysical uncertainties can also be alleviated by requiring that different 
events leads to the same posteriors for ppE parameters (after breaking degeneracies with system 
parameters). Astrophysically, however, there are a limited number of scenarios that could lead to 
corrections in the waveforms that are large enough to interfere with these tests. For comparable 
mass ratio inspirals, this is usually not a problem as the inertia of each binary component is too large 
for any astrophysical environment to affect the orbital trajectory [200J. Magnetohydrodynamic 
effects could affect the merger of neutron-star binaries, but this usually occurs outside of the 
sensitivity band of ground-based interferometers. In extreme mass-ratio inspirals, however, the 
small compact object can be easily nudged away by astrophysical effects, such as the presence 
of an accretion disk |424[ 1237) or a third supermassive black hole [421] . These astrophysical 
effects, however, present the interesting feature that they correct the waveform in a form similar 
to Eq. (196) but with 6ppE < —5. This is because the larger the orbital separation, the stronger 



the perturbations of the astrophysical environment, either because the compact object gets closer 
to the third body or because it leaves the inner edge of the accretion disk and the disk density 
increases with separation. Such effects, however, are not likely to be present in all sources observed, 
as few extreme mass-ratio inspirals are expected to be embedded in an accretion disk or sufficiently 
close to a third body (< 0.1 pc) for the latter to have an effect on the waveform. 

Perhaps the most dangerous systematic is mismodeling, which is due to the use of approxima- 
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tion schemes when constructing waveform templates. For example, in the inspiral one uses the 
post-Newtonian approximation, series expanding and truncating the waveform at a given power of 
orbital velocity. Moreover, neutron stars are usually modeled as test-particles (with a Dirac distri- 
butional density profile), when in reality they have a finite radius, which will depend on its equation 
of state. Such finite-size effects enter at 5 post-Newtonian order, but with a post-Newtonian coef- 
ficient that can be rather large. Ignorance of the post- Newtonian series beyond 3 post-Newtonian 
order can lead to systematics in the determination of physical parameters and possibly also to 
confusion when carrying out ppE-like tests. Much more work is needed to determine the systems 
and signal-to-noise ratios for which such systematics are truly a problem. 

5.3.5 Searching for Non-Tensorial Gravitational Wave Polarizations 

Another way to search for generic deviations from General Relativity is to ask whether any grav- 
itational wave signal detected contains more than the two traditional polarizations expected in 
General Relativity. A general approach to answer this question is through null streams, as dis- 



cussed in Sec. 4.3 This concept was first studied by Giirsel and Tinto |188j and later by Chatter ji 
et al. [94] with the aim to separate false-alarm events from real detections. Chatziioannou, et al [95] 
proposed the extension of the idea of null streams to develop null tests of General Relativity, which 
has recently been implemented by Hayama, et al |199] . 

Given a gravitational wave detection, one can ask whether the data is consistent with two 
polarizations by constructing a null stream through the combination of data streams from 3 or 



more detectors. As explained in Sec. 4.3 such a null stream should be consistent with noise in 
General Relativity, while it would present a systematic deviation from noise if the gravitational 
wave metric perturbation possessed more than 2 polarizations. Notice that such a test would not 
require a template; if one were parametrically constructed, such as in [SS], more powerful null tests 
could be applied to such a null steam. In the future, we expect several gravitational wave detectors 
to be online: the two Ad. LIGO ones in the United States, Ad. VIRGO in Italy, LIGO-India in 
India, and KAGRA in Japan. Given a gravitational wave observation that is detected by all 5 
detectors, one can then construct 3 enhanced General Relativity-null streams, each with power in 
a signal null direction. 

5.3.6 I-Love-Q Tests 

Neutron stars in the slow-rotation limit can be characterized by their mass and radius (to zcroth- 
order in spin), by their moment of inertia (to first-order in spin), and by their quadrupole moment 
and Love numbers (to second-order in spin) . One may expect these quantities to be quite sensitive 
to the neutron star's internal structure, which can be parameterized by its equation of state, ie. the 
relation between its internal pressure and its internal energy density. Since the equation of state 
cannot be well-constrained at super-nuclear densities in the laboratory, one is left with a variety 
of possibilities that predicts different neutron star mass-radius relations. 

Recently, however, Yagi and Yunes [41311414] have demonstrated that there is a relation between 
the moment of inertia (/), the Love numbers (A), and the quadrupole moment (Q), the I-Love- 
Q relations that is essentially insensitive to the equation of state. Figure [Sj shows two of these 
relations (the normalized I-Love and Q-Love relations - see caption) for a variety of equations 
of state, including APR .8j, SLy [135) 1348] . Lattimer-Swesty with nuclear incomprcssibility of 
220MeV (LS220) [MIEQJ, Shen [3451 |346l [30l] , the latter two with temperature of O.OlMeV and 
an electron fraction of 30%, and polytropic equations of state with indices of n = 0.6, 0.8 and 1. 
The bottom panels show the difference between the numerical results and the analytical, fitting 
curve. Observe that all equations of state lead to the same I-Love and Q-Love relations, with 
discrepancies smaller than 1% for realistic neutron-star masses. These results have been recently 
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verified in [272] through the post-Newtonian- Affine approach |1531I273] . which proves the I-Love-Q 
relations hold not even during the inspiral, but also close to plunge and merger. 



M, (APR) |M 




Figure 5: (Top) Fitting curves (solid curve) and numerical results (points) of the universal I- 
Love (left) and Q-Love (right) relations for various equations of state, normalized as / = I/M^g, 
^(tid) ^ AC'^i^/MSg and Q = -Q("'V[^ns('5'A^ns)^]' ^^ns is the neutron-star mass, A^'^^) is the 
tidal Love number, Q^'^°^> is the rotation-induced quadrupole moment, and S is magnitude of the 
neutron star spin angular momentum. The neutron star central density is the parameter varied 
along each curve, or equivalently the NS compactness. The top axis shows the neutron star mass 
for the APR equation of state, with the vertical dashed line showing Mns — IMq. (Bottom) 
Relative fractional errors between the fitting curve and the numerical results. Observe that these 
relations are essentially independent of the equation of state, with loss of universality at the 1% 
level. 

Given the independent measurement of any two members of the TLove-Q trio, one could 
carry out a (null) model-independent and equation-of-state-independent test of General Rela- 
tivity |413l 1414] . For example, assume that electromagnetic observations of the double binary 
pulsar J0737-3039 have measured the moment of inertia to 10% accuracy |253L 12441 I243J . The 
slow-rotation approximation is perfectly valid for this double binary pulsar, due to its relatively 
long spin period. Assume further that a gravitational wave observation of a binary neutron star 
inspiral, with individual masses similar to that of the primary in J0737-3039, manages to measure 
the neutron star tidal Love number to 60% accuracy |413l 1414] . These observations then lead to 
an error box in the LLove plane, which must contain the curve in the left-panel of Fig. [5] 

A similar test could be carried out by using data from only binary pulsar observations or only 
gravitational wave detections. In the case of the latter, one would have to simultaneously measure 
or constrain the value of the quadrupole moment and the Love number, since the moment of inertia 
is not measurable with gravitational wave observations. In the case of the former, one would have 
to extract the moment of inertia and the quadrupole moment, the latter of which will be difficult 
to measure. Therefore, the combination of electromagnetic and gravitational wave observations 
would be the ideal way to carry out such tests. 

Such a test of General Relativity, of course, is powerful only as long as modified gravity theories 
predict I-Love-Q relations that are not degenerated with the General Relativistic ones. Yagi and 
Yunes |413l I414| investigated such a relation in dynamical Chern-Simons gravity to find that such 
degeneracy is only present in the limit Ccs — ^ 0. That is, for any finite value of ^4, the dynamical 
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Chern-Simons I-Love-Q relation differs from that of General Relativity, with the distance to the 
General Relativity expectation increasing for larger Ccs- Yagi and Yunes |4131I414| predicted that a 
test similar to the one described above could constrain dynamical Chern-Simons gravity to roughly 

Ccs* < IOMns ~ 15 km. 

The test described above, of course, only holds provided the I-Love-Q relations are valid, which 
in turn depends on the assumptions made in deriving them. In particular, Yagi and Yunes [4131 
1414] assumed that the neutron stars are uniformly and slowly-rotating, as well as only slightly 
tidally deformed by their rotational velocity or companion. These assumptions would not be 
valid for newly-born neutron stars, which are probably differentially rotating and doing so fast. 
The gravitational waves emitted by neutron star inspirals, however, are expected to have binary 
components that are old and not rapidly spinning, by the time they enter the detector sensitivity 
band [67] . Some short-period, millisecond pulsars may spin at a non-negligible rate, for which the 
normalized moment of inertia, quadrupole moment and Love number would not be independent 
of the rotational angular velocity. If then, however, the above tests should still be possible, since 
binary pulsar observations would also automatically determine the rotational angular velocity, for 
which a unique I-Love-Q relation should exist in General Relativity. 

5.4 Tests of the No-Hair Theorems 

Another important class of generic tests of General Relativity are those that concern the so-called 
no-hair theorems. Since much work has been done on this area, we have decided to separate this 
topic from the main generic tests subsection. In what follows, we describe what these theorems are 
and the possible tests one could carry out with gravitational wave observations emitted by black 
hole binary systems. 

5.4.1 The No-Hair Theorems 

The no-hair theorems state that the only stationary, vacuum solution to the Einstein equations 
that is non-singular outside the event horizon is completely characterized by three quantities: its 
mass M, its spin S and its charge Q. This conclusion is arrived at by combining several different 
theorems. First, Hawking [19511194] proved that a stationary black hole must have an event horizon 
with a spherical topology and that it must be either static or axially symmetric. Israel [2141 1215] 
then proved that the exterior gravitational field of such static black holes is uniquely determined by 
M and Q and it must be given by the Schwarzschild or the Reissner-Nordstrom metrics. Carter ^T\ 
constructed a similar proof for uncharged stationary, axially-symmetric black holes, where this 
time black holes fall into disjoint families, not deformable into each other and with an exterior 
gravitational field uniquely determined by M and S. Robinson |327j and Mazur |275| later proved 
that such black holes must be described by either the Kerr or the Kerr-Newman metric. See 
also [2861 I315J for more details. 

The no-hair theorems apply under a restrictive set of conditions. First, the theorems only apply 
in stationary situations. Black hole horizons can be tidally deformed in dynamical situations, and 
if so, Hawking's theorems [1951 1194| about spherical horizon topologies do not apply. This then 
implies that all other theorems described above also do not apply, and thus, dynamical black holes 
will generically have hair. Second, the theorems only apply in vacuum. Consider, for example, 
an axially symmetric black hole in the presence of a non-symmetrical matter distribution outside 
the event horizon. One might naively think that this would tidally distort the event horizon, 
leading to a rotating, stationary black hole that is not axisymmetric. Hawking and Hartle |198) . 
however, showed that in such a case the matter distribution torques the black hole forcing it to 
spin down, thus leading to a non-stationary scenario. If the black hole is non-stationary, then again 
the no-hair theorems do not apply by the arguments described at the beginning of this paragraph. 
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and thus non-isolated black holes can have hair. Third, the theorems only apply within General 
Relativity, ie. through the use of the Einstein equations. Therefore, it is plausible that black holes 
in modified gravity theories or in General Relativity with singularities outside any event horizons 
(naked singularities) will have hair. 

The no-hair theorems imply that the exterior gravitational field of isolated, stationary, un- 
charged and vacuum black holes (in General Relativity and provided the spacetimc is regular 
outside all event horizons) can be written as an infinite sum of mass and current multipole mo- 
ments, where only two of them are independent: the mass monopole moment M and the current 
dipole moment S. If the no-hair theorems hold, all other multipole moments can be determined 
from pTfiHTMlfTOO] 

Me + iSg = M{iaY , (210) 

where Mi and Se are the £th mass and current multipole moments. Even if the black hole progenitor 
was not stationary or axisymmetric, the no-hair theorems guarantee that any excess multipole 
moments will be shed-off during gravitational collapse [3201 1321| . Eventually, after the black hole 
has settled down and reached an equilibrium configuration, it will be described purely in terms of 
Mq — M and 51 = 5* = Ma? ^ where a is the Kerr spin parameter. 

An astrophysical observation of a hairy black hole would not imply that the no- hair theorems are 
wrong, but rather than one of the assumptions made in deriving these theorems is not appropriate 
to describe Nature. As described above, the three main assumptions are stationarity, vacuum 
and that General Relativity and the regularity condition hold. Astrophysical black holes will 
generically be hairy due to a violation of the first two assumptions, since they will neither be 
perfectly stationary, nor exist in a perfect vacuum. Astrophysical black holes will always suffer 
small perturbations by other stars, electromagnetic fields, other forms of matter, like dust, plasma 



or dark matter, etc, which will induce non-zero deviations from Eq. (210) and thus evade the no- 
hair theorems. In all cases of interest, however, such perturbations are expected to be too small 
to be observable, which is why one argues that even astrophysical black holes should obey the 
no-hair theorems if General Relativity holds. Put another way, an observation of the violation 
of the no-hair theorems would be more likely to indicate a failure of General Relativity in the 
strong-field, than an unreasonably large amount of astrophysical hair. 

Tests of the no-hair theorems come in two flavors: through electromagnetic observations [2211 
[2221, V2M, "225] and through gravitational wave observations [3341 |335l [Ml [HSl HQI HSl [2591 ESS 
[432, 383, 3821 1162|, 13811 I328J . The former rely on radiation emitted by accelerating particles in an 
accretion disk around black holes. Such tests, however, are not clean as they require the modeling 
of complicated astrophysics, with matter and electromagnetic fields. Gravitational wave tests are 
clean in that respect, but unlike electromagnetic tests, they cannot be carried out yet due to 
lack of data. Other electromagnetic tests of the no-hair theorems exist, for example through the 
observation of close stellar orbits around Sgr A* [280) 12811 13371 and pulsar-black hole binaries [392j . 
but these cannot yet probe the near-horizon, strong-field regime, since electromagnetic observations 
cannot yet resolve horizon scales. See [323| for reviews on this topic. 

5.4.2 Extreme Mass-Ratio Tests of The No-Hair Theorem 

Gravitational wave tests of the no-hair theorems require the detection of either extreme mass- 
ratio inspirals or the ringdown of comparable-mass black hole mergers with future space-borne 
gravitational wave detectors [201 [21] • Extreme mass-ratio inspirals consist of a stellar- mass compact 
object spiraling into a supermassive black hole in a generic orbit within astronomical units from 
the event horizon of the supermassive object (32]. These events outlive the observation time of 
future detectors, emitting millions of gravitational wave cycles, with the stellar-mass compact 
object essentially acting as a tracer of the supermassive black hole spacetime [358j . Ringdown 
gravitational waves are always emitted after black holes merge and the remnant settles down 
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into its final configuration. During the ringdown, tfie higlily-distorted remnant radiates all excess 
degrees of freedom and this radiation carries a signature of whether the no-hair theorems hold in 
its quasi- normal mode spectrum (see e.g. j62j for a recent review). 

Both electromagnetic and gravitational wave tests need a metric with which to model accretion 
disks, quasi-periodic oscillations, or extreme mass-ratio inspirals. One can classify these metrics 
as direct or generic, paralleling the discussion in Sec. |5.2[ Direct metrics are exact solutions to a 
specific set of field equations, with which one can derive observables. Examples of such metrics are 
the Manko-Novikov metric [271] and the slowly-spinning black hole metric in dynamical Chern- 
Simons gravity |427| . When computing observables with these metrics, one usually assumes that 
all radiative and dynamical process (e.g. the radiation-reaction force) are as predicted in General 
Relativity. Generic metrics are those that parametrically modify the Kerr spacetime, such that 
for certain parameter choices one recovers identically the Kerr metric, while for others, one has a 
deformation of Kerr. Generic metrics can be further classified into two subclasses, Ricci-fiat versus 
non-Ricci-flat, depending on whether they satisfy i?^i/ = 0. 

Let us first consider direct metric tests of the no-hair theorem. The most studied direct metric 
is the Manko-Novikov one, which although an exact, stationary and axisymmetric solution to the 
vacuum Einstein equations, does not represent a black hole, as the event horizon is broken along the 
equator by a ring singularity |271| . Just like the Kerr metric, the Manko-Novikov metric possesses 
an ergoregion, but unlike the former, it also possesses regions of closed time-like curves that overlap 
the ergoregion. Nonetheless, an appealing property of this metric is that it deviates continuously 
from the Kerr metric through certain parameters that characterize the higher multiple moments 
of the solution. 

The first geodesic study of Manko-Novikov spacetimes was carried out by Gair, et al |163) . 
They found that there are two ring-like regions of bound orbits; an outer one where orbits look 
regular and integrable, as there exists four isolating integrals of the motion; and an inner one where 
orbits are chaotic and thus ergodic. Gair, et al [163) suggested that orbits that transition from the 
integrable to the chaotic region would leave a clear observable signature in the frequency spectrum 
of the emitted gravitational waves. However, they also noted that chaotic regions exist only very 
close to the central body and are probably not astrophysically accessible. The study of Gair, et 
al [163] was recently confirmed and followed up by Contopoulos, et al |109| . They studied a wide 
range of geodesies and found that, in addition to an inner chaotic region and an outer regular 
region, there are also certain Birkhoff islands of stability. When an extreme mass-ratio inspiral 
traverses such a region, the ratio of resonant fundamental frequencies would remain constant in 
time, instead of increasing monotonically. Such a feature would impact the gravitational waves 
emitted by such a system, and it would signal that the orbit equations are non-integrable and the 
central object is not a Kerr black hole. 

The study of chaotic motion in geodesies of non-Kerr spacetimes is by no means new. Chaos has 
also been found in geodesies of Zipoy-Voorhees-Weyl and Curzon spacetimes |354ll355j . of perturbed 
Schwarzschild spacetimes |257| , of Schwarzschild spacetimes with a dipolar halo |256l 12581 I187J of 
Erez- Rosen spacetimes |186j . and of deformed generalizations of the Tomimatsy-Sato spacetime 
|139j . One might worry that such chaotic orbits will depend on the particular spacetime considered, 
but recently Apostolatos, et al [28] and Lukes- Gerakopoulos, et al }266j have argued that the 
Birkhoff islands of stability are a general feature. Although the Kolmogorov, Arnold, and Moser 
theorem [2391 [5^ 1288) states that phase orbit tori of an integrable system are only deformed if 
the Hamiltonian is perturbed, the Poincare-Birkhoff theorem [262" states that resonant tori of 
integrable systems actually disintegrate, leaving behind a chain of Birkhoff islands. These islands 
are only characterized by the ratio of winding frequencies that equals a rational number, and thus, 
they constitute a distinct and generic feature of non-integrable systems [28( I266J . Given an extreme 
mass-ratio gravitational wave detection, one can monitor the ratio of fundamental frequencies and 
search for plateaus in their evolution, which would signal non-integrability. Of course, whether 
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detectors can resolve such plateaus depends on the initial conditions of the orbits and the physical 
system under consideration (these determine the thickness of the islands) , as well as the mass ratio 
(this determines the radiation-reaction timescale) and the distance and mass of the central black 
hole (this determines the signal-to-noise ratio). 

Another example of a direct metric test of the no-hair theorem is through the use of the slowly- 
rotating dynamical Chern-Simons black hole metric |427j . Unlike the Manko-Novikov metric, the 
dynamical Chern-Simons one does represent a black hole, ie. it possesses an event horizon, but 
it evades the no-hair theorems because it is not a solution to the Einstein equations. Sopuerta 
and Yunes |353j carried out the first extreme mass-ratio inspiral analysis when the background 
supermassive black hole object is taken to be such a Chern-Simons black hole. They used a semi- 
relativistic model |332| to evolve extreme mass-ratio inspirals and found that the leading-order 
modification comes from a modification to the geodesic trajectories, induced by the non-Kerr 
modifications of the background. Because the latter correspond to a strong-field modification to 
General Relativity, modifications in the trajectories are most prominent for zoom-whirl orbits, as 
the small compact object zooms around the supermassive black hole in a region of unstable orbits, 
close to the event horizon. These modifications were then found to propagate into the gravitational 
waves emitted, leading to a dephasing that could be observed or ruled out with future gravitational 
wave observations to roughly the horizon scale of the supermassive black hole, as has been recently 
confirmed by Canizares et al |86j . These studies, however, may be underestimates, given that they 
treat the black hole background in dynamical Chern-Simons gravity only to first-order in spin. 

A final example of a direct metric test of the no-hair theorems is by considering black holes that 
are not in vacuum. Barausse, et al [47 studied extreme mass-ratio inspirals in a Kerr black hole 
background that is perturbed by a self-gravitating, homogeneous torus that is compact, massive 
and close to the Kerr black hole. They found that the presence of this torus impacts the gravita- 
tional waves emitted during such inspirals, but only weakly, making it difficult to distinguish the 
presence of matter. Yunes, et al |424j and Kocsis, et al j237j carried out a similar study, where this 
time they considered a small compact object inspiraling completely within a geometrically thin, 
radiation-pressure dominated accretion disk. They found that disk-induced migration can modify 
the radiation-reaction force sufficiently so as to leave observable signatures in the waveform, pro- 
vided the accretion disk is sufficiently dense in the radiation-dominated regime and a gap opens 
up. These tests of the no-hair theorem, however, will be rather difficult as most extreme mass-ratio 
inspirals are not expected to be in an accretion disk. 

Let us now consider generic metric tests of the no-hair theorem. Generic Ricci-flat deformed 
metrics will lead to Laplace-type equations for the deformation functions in the far-field since they 
must satisfy i?^^ == to linear order in the perturbations. The solution to such an equation can 
be expanded in a sum of mass and current multipole moments, when expressed in asymptotically 
Cartesian and mass-centered coordinates [369 . These multipoles can be expressed via [105^ i383l 

Mi + iSi^MiiaY + SMi + idSe, (211) 

where SMg and SSe are mass and current multipole deformations. Ryan [33411335] showed that the 
measurement of three or more multipole moments would allow for a test of the no-hair theorem. 
Generic non-Ricci flat metrics, on the other hand, will not necessarily lead to Laplace-type equa- 



tions for the deformation functions in the far field, and thus, the far-field solution and Eq. (211) 
will depend on a sum of i and m multipole moments. 

The first attempt to construct a generic, Ricci-flat metric was by Collins and Hughes |105j : the 
so-called bumpy black hole metric. In this approach, the metric is assumed to be of the form 

g^.^g'-^.''"'>+^Sg,,., (212) 

where e ^ 1 is a book keeping parameter that enforces that 6gfj_^ is a perturbation of the Kerr 
background. This metric is then required to satisfy the Einstein equations linearized in e, which 
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then leads to differential equations for the metric deformation. Collins and Hughes |105j assumed 
a non-spinning, stationary spacetime, and thus Sg^^ only possessed two degrees of freedom, both of 
which were functions of radius only: '0i(''): which must be a harmonic function and which changes 
the Newtonian part of the gravitational field at spatial infinity; and 71(7") which is completely 
determined through the linearized Einstein equations once tpi is specified. One then has the 
freedom to choose how to prescribe V'l and Collins and Hughes investigated [105j two choices 
that corresponded physically to point-like and ring-like naked singularities, thus violating cosmic 
censorship |312j . Vigeland and Hughes |383| and Vigeland |382j then extended this analysis to 
stationary, axisymmetric spacetimes via the Newman- Janis method |294l I136J , showed how such 



metric deformations modify Eq. (211), and computed how these bumps imprint themselves onto 
the orbital frequencies and thus the gravitational waves emitted during an extreme mass-ratio 
inspiral. 

That the bumps represent unphysical matter should not be a surprise, since by the no-hair 
theorems, if the bumps are to satisfy the vacuum Einstein equations they must either break sta- 
tionarity or violate the regularity condition. Naked singularities are an example of the latter. A 
Lorentz-violating massive field coupled to the Einstein tensor is another example 11401 . Gravita- 
tional wave tests with bumpy black holes must then be understood as null tests: one assumes the 
default hypothesis that General Relativity is correct and then sets out to test whether the data 
rejects or fails to reject this hypothesis (a null hypothesis can never be proven). Unfortunately, 
however, bumpy black hole metrics cannot parameterize spacetimes in modified gravity theories 
that lead to corrections in the field equations that are not proportional to the Ricci tensor, such 
as for example in dynamical Chern-Simons or in Einstein-Dilaton-Gauss-Bonnet modified gravity. 

Other bumpy black hole metrics have also been recently proposed. Glampedakis and Babak |176| 
proposed a different type of stationary and axisymmetric bumpy black hole through the Hartle- 
Thorne metric |192j . with modifications to the quadrupole moment. They then constructed a 
"kludge" extreme mass-ratio inspiral waveform and estimated how well the quadrupole deforma- 
tion could be measured ^]. This metric, however, is valid only when the supermassive black hole is 
slowly-rotating, as it derives from the Hartle-Thorne ansatz. Recently, Johansen and Psaltis [223] 
proposed yet another metric to represent bumpy stationary and spherically symmetric spacetimes. 
This metric introduces one new degree of freedom, which is a function of radius only and as- 
sumed to be a series in M/r. Johansen and Psaltis then rotated this nietic via the Newman-Janis 
method |294| I136J to obtain a new bumpy metric for axially symmetric spacetimes. Such a metric, 
however, possesses a naked ring singularity on the equator, and naked singularities on the poles. 
As before, none of these bumpy metrics can be mapped to known modified gravity black hole 
solutions, in the Glampedakis and Babak case [176} because the Einstein equations are assumed 
to hold to leading order in the spin, while in the Johansen and Psaltis case |223) because a single 
degree of freedom is not sufficient to model the three degrees of freedom contained in stationary 
and axisymmetric spacetimes |362l 1381) . 

The only generic non-Ricci-flat bumpy black hole metrics so far is that of Vigeland, Yunes 
and Stein [381j. They allowed generic deformations in the metric tensor, only requiring that the 
new metric perturbatively retained the Killing symmetries of the Kerr spacetime: the existence 
of two Killing vectors associated with stationarity and axisymmetry, as well as the perturbative 
existence of a Killing tensor (and thus a Carter-like constant), at least to leading order in the 
metric deformation. Such requirements imply that the geodesic equations in this new background 
are fully integrable, at least perturbatively in the metric deformation, which then allows one to 
solve for the orbital motion of extreme mass-ratio inspirals by adapting previously existing tools. 
Brink [TTl [78l [79l [8OI [81] studied the existence of such a second-order Killing tensor in generic, 
vacuum, stationary and axisymmetric spacetimes in Einstein's theory and found that these are 
difficult to construct exactly. By relaxing this exact requirement, Vigeland, Yunes and Stein [381] 
found that the existence of a perturbative Killing tensor poses simple differential conditions on 
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the metric perturbation that can be analytically solved. Moreover, they also showed how this new 
bumpy metric can reproduce all known modified gravity black hole solutions in the appropriate 
limits, provided these have an at least approximate Killing tensor; thus, these metrics are still 
vacuum solutions even though i? 7^ 0, since they satisfy a set of modified field equations. Although 
unclear at this junction, it seems that the imposition that the spacetime retains the Kerr Killing 
symmetries leads to a bumpy metric that is well-behaved everywhere outside the event horizon 
(no singularities, no closed-time-like curves, no loss of Lorentz signature) |384| . Recently, Gair and 
Yunes |162) studied how the geodesic equations are modified for a test-particle in a generic orbit in 
such a spacetime and showed that the bumps are indeed encoded in the orbital motion, and thus, 
in the gravitational waves emitted during an extreme mass-ratio inspiral. 

One might be concerned that such no-hair tests of General Relativity cannot constrain modified 
gravity theories, because Kerr black holes can also be solutions in the latter [324) . This is indeed 
true provided the modified field equations depend only on the Ricci tensor or scalar. In Einstein- 
Dilaton-Gauss-Bonnet or dynamical Chern-Simons gravity, the modified field equations depend on 



the Riemann tensor, and thus, Ricci- flat metric need not solve these modified set |434j . Moreover, 
just because the metric background is identically Kerr does not imply that inspiral gravitational 
waves will be identical to those predicted in General Relativity. All studies carried out to date, 
be it direct metric tests or generic metric tests, assume that the only quantity that is modified is 
the metric tensor, or equivalently, the Hamiltonian or binding energy. Inspiral motion, of course, 
does not depend just on this quantity, but also on the radiation-reaction force that pushes the 
small object from geodesic to geodesic. Moreover, the gravitational waves generated during such 
an inspiral depend on the field equations of the theory considered. Therefore, all metric tests 
discussed above should be considered as partial tests. In general, strong-field modified gravity 
theories will modify the Hamiltonian, the radiation-reaction force and the wave generation. 

5.4.3 Ringdown Tests of The No-Hair Theorem 

Let us now consider tests of the no-hair theorems with gravitational waves emitted by comparable- 
mass binaries during the ringdown phase. Gravitational waves emitted during ringdown can be 
described by a superposition of exponentially damped sinusoids |63| : 



h+\ 
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(213) 
where r is the distance from the source to the detector, the asterisk stands for complex conjugation, 
the real mode amplitudes Ai^m,n and -A^ ^ „ and the real phases (j>nem and (j>nem depend on the 
initial conditions, Simn are spheroidal functions evaluated at the complex quasinormal ringdown 
frequencies ujnim = '^'^fnim + i/Tnim, and the real physical frequency fnim and the real damping 
times Tnim are both functions of the mass M and the Kerr spin parameter a only, provided the 
no-hair theorems hold. These frequencies and damping times can be computed numerically or 
semi-analytically, given a particular black hole metric (see |62| for a recent review). The Fourier 
transform of a given (i, m, n) mode is p3j 



^V'''"'"^(^) = -A™„ [e'<""5,„™6+(c.) + e-^-"^;_6_(c.)] , (214) 
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where we have defined ■AJ^m^^'^''"" = ■^imne^'^'"^" ±A'e *'^'?mi as well as the Lorentzian bmctions 

^±M = -. r-^ ^ (216) 
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Ringdown gravitational waves will all be of the form of Eq. (2131 provided that the characteristic 
nature of the differential equation that controls the evolution of ringdown modes is not modified, 
ie. provided that one only modifies the potential in the Teukolsky equation or other subdominant 
terms, which in turn depend on the modified field equations. 

Tests of the no-hair theorems through the observation of black hole ringdown date back to 
Detweiler |131| . and it was recently worked out in detail by Dreyer, et al. |I3Z]- Let us first 
imagine that a single complex mode is detected uJij^mini and one measures separately its real 
and imaginary parts. Of course, from such a measurement, one cannot extract the measured 
harmonic triplet {ii,mi,ni), but instead one only measures the complex frequency w^imini- This 
information is not sufficient to extract the mass and spin angular momentum of the black hole 
because different quintuplets {M,a,£,m,n) can lead to the same complex frequency uji-^mini- The 
best way to think of this is graphically: a given observation of w^ ^ ^ traces a line in the complex 

^iimini = Mu!f l^ „ plane; a given {£, m, n) triplet defines a complex frequency uoimn that also 
traces a curve in the complex ^imn plane; each intersection of the measured line ^i-^^mini with ^imn 
defines a possible doublet (Af, a); since different {(., m, n) triplets lead to different uJimn curves and 
thus different intersections, one ends up with a set of doublets 5*1, out of which only one represents 
the correct black hole parameters. We thus conclude that a single mode observation of ringdown 
gravitational waves is not sufficient to test the no- hair theorem |1371 163] . 

Let us then imagine that one has detected two complex modes, w^^mini and cu^^^^n^- Each 
detection leads to a separate line Qi-^mini and i7^2m2n2 in the complex plane. As before, each 
{n,i, m) triplet leads to separate curves ^e,mn which will intersect with both ^i-^^mini and ^i^rn2n2 
in the complex plane. Each intersection between ^tmn and rif^m^nj^ leads to a set of doublets 5*1, 
while each intersection between ri^mn and ^i^rn2n2 leads to another set of doublets 5*2. If the no- 
hair theorems hold, however, sets 5*1 and 5*2 must have at least one element in common. Therefore, 
a two- mode detection allows for tests of the no- hair theorem [1371 [55] . When dealing with a quasi- 
circular binary black hole inspiral within General Relativity, however, one knows that the dominant 
mode is the I = m — 2 one. In such a case, the observation of this complex mode by itself allows 
one to extract the mass and spin angular momentum of the black hole. Then, the detection of the 
real frequency in an additional mode can be used to test the no-hair theorem ^63t i60j . 

Although the logic behind these tests is clear, one must study them carefully to determine 
whether all systematic and statistical errors are sufficiently under control so that they are feasible. 
Bcrti et al 63j .60J investigated such tests carefully through a frequentist approach. First, they 
found that a matched-filtering type analysis with two-mode ringdown templates would increase 
the volume of the template manifold by roughly three orders of magnitude. A better strategy 
then is perhaps to carry out a Bayesian analysis, like that of Gossan et al. [2271 1178] : through 
such a study one can determine whether a given detection is consistent with a two-mode or a 
one-mode hypothesis. Berti et al [531 ISO] also calculated that a signal-to-noise ratio of O(IO^) 
would be sufficient to detect the presence of two modes in the ringdown signal and to resolve 
their frequencies, so that no-hair tests would be possible. Strong signals are necessary because 
one must be able to distinguish at least two modes in the signal. Unfortunately, however, whether 
the ringdown leads to such strong signal-to-noise ratios and whether the sub-dominant ringdown 
modes are of a sufficiently large amplitude depends on a plethora of conditions: the location of 
the source in the sky, the mass of the final black hole, which depends on the rest mass fraction 
that is converted into ringdown gravitational waves (the ringdown efficiency), the mass ratio of the 
progenitor, the magnitude and direction of the spin angular momentum of the final remanent and 
probably also of the progenitor and the initial conditions that lead to ringdown. Thus, although 
such tests are possible, one would have to be quite fortunate to detect a signal with the right 
properties so that a two-mode extraction and a test of the no-hair theorems is feasible. 
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5.4.4 The Hairy Search for Exotica 

Another way to test General Relativity is to modify the matter-sector of the theory through the 
introduction of matter corrections to the Einstein- Hilbert action that violate the assumptions made 
in the no-hair theorems. More precisely, one can study whether gravitational waves emitted by 
binaries composed of strange stars, like quark stars, or horizonless objects, such as boson stars or 
gravastars, are different from waves emitted by more traditional neutron star or black hole binaries. 
In what follows, we will describe such hairy tests of the existence of compact exotica. 

Boson stars are a classic example of a compact object that is essentially indistinguishable from a 
black hole in the weak field, but which differs drastically from one in the strong field due to its lack 
of an event horizon. A boson star is a coherent scalar-field configuration supported against gravi- 
tational collapse by its self-interaction. One can construct several Lagrangian densities that would 
allow for the existence of such an object, including mini-boson stars |1581I159] , axially-symmetric 
solitons j336] . and nonsolitonic stars supported by a non-canonical scalar potential energy |106) . 
Boson stars are well-motivated from fundamental theory, since they are the gravitationally-coupled 
limit of q-balls [1011 I246J , a coherent scalar condensate that can be described classically as a non- 
topological soliton and that arises unavoidably in viable supersymmetric extensions of the standard 
model }245, . In all studies carried out to date, boson stars have been studied within General Rel- 
ativity, but they are also allowed in scalar-tensor theories |41) . 

At this junction, one should point out that the choice of a boson star is by no means special; 
the key point here is to select a straw-man to determine whether gravitational waves emitted 
during the coalescence of compact binaries are sensitive to the presence of an event horizon or the 
evasion of the no-hair theorems induced by a non- vacuum spacetime. Of course, depending on the 
specific model chosen, it is possible that the exotic object will be unstable to binary evolution or 
even to its own rotation. For example, in the case of an extreme mass-ratio inspiral, one could 
imagine that as the small compact object enters the boson star's surface, it will accrete scalar 
field, forcing the boson star to collapse into a black hole. Alternatively, one can imagine that 
as two supermassive boson stars merge, the remnant might collapse into a black hole, emitting 
the usual General Relativity quasinormal modes. What is worse, even when such objects are in 
isolation, they are unstable under small perturbations if their angular momentum is large, possibly 
leading to gravitational collapse into a black hole or possibly a scalar explosion [88l [89] . Since most 
astrophysical black hole candidates are believed to have high spins, such instabilities somewhat 
limit the interest of horizonless objects. Even then, however, the existence of slowly spinning or 
non spinning horizonless compact objects cannot be currently ruled out by observation. 

Boson stars evade the no-hair theorems within General Relativity because they arc not vacuum 
spacetimes, and thus, their metric and quasinormal mode spectrum cannot be described by just 
their mass and spin angular momentum; one also requires other quantities intrinsic to the scalar- 
field energy momentum tensor, scalar hair. Therefore, as before, two types of gravitational wave 
tests for scalar hair have been proposed: extreme mass-ratio inspiral tests and ringdown tests. As 
for the former, several studies have been carried out that considered a supermassive boson star 
background. Kesden, et al [234| showed that stable circular orbits exist both outside and inside 
of the surface of the boson star, provided the small compact object interacts with the background 
only gravitationally. This is because the effective potential for geodesic motion in such a boson 
star background lacks the Schwarzschild-like singular behavior at small radius, instead turning over 
and allowing for a new minimum. Gravitational waves emitted in such a system would then stably 
continue beyond what one would expect if the background had been a supermassive black hole; 
in the latter case the small compact object would simply disappear into the horizon. Kesden et 
al |234j found that orbits inside the boson star exhibit strong precession, exciting high frequency 
harmonics in the waveform, and thus allowing one to easily distinguish between such boson stars 
from black hole backgrounds. 
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Just like the inspiral phase is modified by the presence of a boson star, the merger phase is 
also greatly altered, but this must be treated fully numerically. A few studies have found that the 
merger of boson stars leads to a spinning bar configuration that either fragments or collapses into a 
Kerr black hole |304l 1303) . Of course, the gravitational waves emitted during such a merger will be 
drastically different from those produced when black holes merge. Unfortunately, the complexity 
of such simulations makes predictions difficult for any one given example, and the generalization 
to other more complicated scenarios, such as theories with modified field equations, is currently 
not feasible. 

Recently, Pani et al [3051 I306J revisited this problem, but instead of considering a supermassive 
boson star, they considered a gravastar. This object consists of a Schwarzschild exterior and a de 
Sitter interior, separated by an infinitely thin shell with finite tension |274| l93]. Pani et al |306j 
calculated the gravitational waves emitted by a stellar-mass compact object in a quasi-circular 
orbit around such a gravastar background. In addition to considering a different background, 
Pani et al used a radiative-adiabatic waveform generation model to describe the gravitational 
waves [3l6l[2Ml[2l0llll9l|4l7l|420], instead of the kludge scheme used by Kesden et al [44l l40l liTZ] . 
Pani el al |306) concluded that the waves emitted during such inspirals are sufficiently different 
that they could be used to discern between a Kerr black hole and a gravastar. 

On the ringdown side of no-hair tests, several studies have been carried out. Berti and Car- 
doso [HI] calculated the quasi-normal mode spectrum of boson stars. Chirenti and Rezzolla ^93] 
studied the non-radial, axial perturbations of gravastars, and Pani, et al |305j the non-radial, axial 
and polar oscillations of gravastars. Medved, et al [27711278] considered the quasinormal ringdown 
spectrum of skyrmion black holes J349J. In all cases, it was found that the quasi- normal mode 
spectrum of such objects could be used to discern between them and Kerr black holes. Of course, 
such tests still require the detection of ringdown gravitational waves with the right properties, such 
that more than one mode can be discerned and extracted from the signal (see Sec. 5.4.3). 



6 Musings About the Future 

Gravitational waves hold the key to testing Einstein's theory of General Relativity to new exiting 
levels in the previously unexplored strong-field regime. Depending on the type of wave that is 
detected, eg. compact binary inspirals, mergers, ringdowns, continuous sources, supcrnovae, etc, 
different tests will be possible. Irrespective of the type of wave detected, two research trends seem 
currently to be arising: direct tests and generic tests. These trends aim at answering different ques- 
tions. With direct tests, one wishes to determine whether a certain modified theory is consistent 
with the data. Generic tests, on the other hand, ask whether the data is statistically consistent 
with our canonical believes. Or put another way: are there any statistically significant deviations 
present in the data from what we expected to observe? This approach is currently used in cosnio- 
logical observations, for example by the WMAP team, and it is particularly well-suited when one 
tries to remain agnostic as to which is the correct theory of Nature. Given that we currently have 
no data in the strong-field, it might be too restrictive to assume General Relativity is correct prior 
to verifying that this is the case. 

What one would like to believe is that gravitational waves will be detected by the end of 
this decade, either through ground-based detectors or through pulsar timing arrays. Given this, 
there is a concrete effort to develop the proper formalism and implementation pipelines to test 
Einstein's theory once data becomes available. Currently, the research groups separate into two 
distinct classes: theory and implementation. The theory part of the research load is being carried 
out at a variety of institutions without a given focal point. The implementation part is being 
done mostly within the LIGO Scientific Collaboration and the pulsar timing consortia. Cross- 
communication between the theory and implementation groups has recently flourished and one 
expects more interdisciplinary work in the future. 

Many accomplishments have been made in the past 50 years that it is almost impossible to 
list them here. From the implementation side, perhaps one of the most important ones is the 
actual construction and operation of the initial LIGO instruments at design sensitivity in all 
of its frequency domain. This is a tremendously important engineering and physics challenge. 
Similarly, the constructing of impressive pulsar timing arrays, and the timing of these pulses to 
nano-second precision is an instrumental and data analysis feat to be admired. Without these 
observatories no waves would be detectable in the future, and of course, no tests of Einstein's 
theory would be feasible. On the theory side, perhaps the most important accomplishment has 
been the understanding of the inspiral phase to really high post-Newtonian order and the merger 
phase with numerical simulations. The latter, in particular, had been an unsolved problem for over 
50 years, until very recently. It is these accomplishments that then allow us to postulate modifled 
inspiral template families, since we understand what the General Relativity expectation is. This is 
particularly true if one is considering small deformations away from Einstein's theory, as it would 
be impossible to perturb about an unknown solution. 

The main questions that are currently at the forefront are the following. On the theory side of 
things, one would wish to understand the inspiral to high post-Newtonian order in certain strong- 
field modifications to General Relativity, like dynamical Chern-Simons gravity or Einstein-Dilaton- 
Gauss-Bonnet theory. One would also like to investigate theories with preferred frames, such as 
Einstein- Aether theory or Horava-Lifshitz gravity, which will lead to Lorentz violating observables. 
Understanding these theories to high post-Newtonian order is particularly important for those 
that predict dipolar gravitational emission, such as Einstein-Dilaton-Gauss-Bonnet theory. Such 
corrections dominate over Einstein's quadrupole emission at sufficiently low velocities. 

A full inspiral- merger-ringdown template is of course not complete unless we also understand the 
merger. This would require full numerical simulations, which are of course very taxing even within 
General Relativity. Once one modifies the Einstein field equations, the characteristic structure of 
the evolution equations will also likely change, and it is unclear whether the standard evolution 
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methods will continue to work. Moreover, when dealing with the merger phase, one is usually 
forced to treat the modified theory as exact, instead of as an effective theory. Without the latter, 
it is likely that certain modified theories will not have a well-posed initial value problem, which 
would force any numerical evolution to fail. Of course, one could order-reduce these equations and 
then use these to evolve black hole spacetimes. Much work still remains to be done to understand 
whether this is feasible. 

On the implementation side of things, there is also much work that remains to be done. Cur- 
rently, efforts are only beginning on the implementation of Bayesian frameworks for hypothesis 
testing. This seems today like one of the most promising approaches to testing Einstein's the- 
ory with gravitational waves. Current studies concentrate mostly on single-detectors, but by the 
beginning of the next decade we expect 4 or 5 detectors to be online, and thus, one would like 
to see these implementations extended. The use of multiple detectors also opens the door to the 
extraction of new information, such as multiple polarization modes, a precise location of the source 
in the sky, etc. Moreover, the evidence for a given model increases dramatically if the event is 
observed in several detectors. One therefore expects that the strongest tests of General Relativity 
will come from leveraging the data from all detectors in a multiply coincident event. 

Ultimately, research is moving toward the construction of robust techniques to test Einstein's 
theory. A general push is currently observed toward the testing of general principles that serve as 
foundations of General Relativity. This allows one to answer general questions like for example: 
does the graviton have a mass? are compact objects represented by the Kerr metric and the no- hair 
theorems satisfied? does the propagating metric perturbation possess only two transverse-traceless 
polarization modes? what is the rate of change of a binary's binding energy? do naked singularities 
exist in Nature and are orbits chaotic? Is Lorentz-violation present in the propagation of gravitons? 
These are examples of questions that can be answered once gravitational waves are detected. The 
more questions of this type that are generated and the more robust the methods to answer them 
are, the more stringent the test of Einstein's theories and the more information we will obtain 
about the gravitational interaction in a previously unexplored regime. 
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